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Microwave Apparatus for the Measurement of the Refraction, 
Dispersion and Absorption of Gases at Relatively High Pressure. 


A. BATTAGLIA, F. BRUIN (*) and A. GozzINI 


Istituto di Fisica dell’ Università - Pisa 


(ricevuto il 23 Giugno 1957) 


Summary. — An apparatus of high sensitivity is described using a single 
frequency modulated klystron oscillator. The resonance curve of a sample 
cavity is compared to the one of a reference cavity by means of an elec- 
tronic device using microsecond pulse technique. Without special pre- 
cautions the sensitivity for measurements of the power absorption co- 
efficient is 10-> m- and 10-$ for the refractive index. 


1. General outline of apparatus. 


The apparatus described below was designed to search for still unknown 
microwave absorption bands of gases and vapours at relatively high pressures, 
when resonance lines are too broad to apply the conventional waveguide 
techniques of microwave spectroscopy. The application of the latter technique, 
which is extremely accurate and sensitive but also rather time absorbing, 
could then be restricted to those frequency regions selected in advance by 
our apparatus. Apart from this the set-up turned out to be very useful as 
well to investigate pressure effects of known absorption bands. 

For gases at relatively high pressure over a small band in the spectrum 
the absorption and refraction are independent of frequency so that these 


(*) Now at the Natuurkundig Laboratorium. Universiteit van Amsterdam. 
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quantities may be measured with resonance cavities. The small space re- 
quired (as compared to a waveguide) has the practical advantage that only 
relatively little gas is needed and that the device may easily be temperature 
controlled. Absorption measurements of this kind were first carried out by 
BLEANEY and PENROSE (!*). Their method was further developed by BIRN- 
BAUM, KRYDER and LYONS (*), by ZIEMAN (*), by HUGHES and ARMSTRONG (5) 
and by JASINSKY and BERRY (6). 

One of the main difficulties when using cavities for the measurement of 
small absorptions or refractions is the drift and the low frequency fluctuation 
of the frequency and the amplitude of the wave generated by the klystron 
oscillator: The change in frequency during a comparison of an empty cavity 
with a cavity containing the sample limits the sensitivity. Of course the 
sensitivity can be increased when the effect of these fluctuations on the measu- 
rement can be diminished. Usually this is done by stabilizing one or two 
klystron oscillators in frequency against the applied cavity resonators. In the 
present apparatus, however, high sensitivity is achieved more simply, namely 
by comparing both cavities in a time interval of only some microseconds. 

In this way low frequency 


a fluctuations and drift of the 
ema 
; ee Al rate] — klystron are of no importance 
He I mplhirier . +]: * 
=; “3 2 and no electronic stabilization 
Snr of the klystron frequency is 
e |APAffer 204. necessary. Early construct- 


Modulator 


| C REPAS ES ions of the present apparatus 

X-Axis may be found in short notes. 
by GOZZINI, PoLAcco (7) and 

La Puiser Cas BATTAGLIA (8). We will now 


give a more detailed and up 


cas to date description. 


Synchroscope A block diagram of a set- 
‘up is shown in Fig. 1. A klys- 


Fig. 1. — Block diagram of the microwave and È ; 
electronic part of an apparatus used to measure tron K is modulated in fre- 
small absorptions and refractions in gases. quency on the reflector by a 


B. BLEANEY and R. P. PENROSE: Proc. Roy. Soc., 189, 358 (1947). 
B. BLEANEY and R. P. PENROSE: Proc. Phys. Soc., 60, 83 (1948). 
G. BIRNBAUM, S. J. KRYDER and H. Lyons: Journ. Appl. Phys., 22, 95 (1951). 
C. M. Zieman: Journ. Appl. Phys., 23, 154, (1952). 
J. V. Hueues and: H. L. ARMSTRONG: Journ. Appl. Phys., 23, 501 (1952). 
W. JASINSKY and J. A. Berry: Proc. I.E.E., 101, Pt. HI, 337 (1954). 

(7) A. Gozzini and E. Poracco: Oompt. Rend. Acad. Sci. (Paris), 237, 1497. 
1652 (1953). 
(8) A. BATTAGLIA and A. Gozzini: Onde Elect., 35, 500 (1955). 
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triangular voltage of small amplitude and of low frequeney. A transmission 
cavity €, is tuned to the mean frequency », and we suppose a second trans- 
mission cavity 03, coupled via a crystal multiplier, to be tuned to a fre- 
quency 2»,. The resonance curves detected by the crystals +, and x, are am- 
plified considerably. The resulting resonance curves may be viewed on two 
small monitor oscilloscopes M, and M,: By means of the strongly amplified 
signals two discriminators both fire a pulse at the moment when the reso- 
nance curve of interest passes a certain fraction of its maximum height. 
These pulses are shaped in the pulsers to have a steep rise and a length of 2 us. 
The pulse from channel 1 is used to trigger a synchroscope, the high speed time 
base of which is calibrated in microseconds by intensity modulation. The 
pulse from channel 2 is shown on the screen once for increasing and once 
for decreasing klystron modulation voltage. The pulses on the screen may 
be distinguished by applying an additional square wave voltage to the y-axis 
of the synchroscope as is indicated in the block diagram. The pulses may 
be made coincident by tuning either C, or C;. They may be shifted together 
along the x-axis by the variation of a delay in channel 1. 

When a little gas is introduced into ©, or C; the pulses will be considerably 
displaced. Using well known formulae for cavity resonators (1) it will be shown 
that a shift of their « barycentre » is proportional to the power absorption 2a” 
of the gas, whereas a shift of one pulse relative to the other is proportional 
to the real part of the refractive index n’. Gases could be investigated at two 
different frequencies, one being twice the other, simply by putting them in 
either €, or C:. 

If sufficient microwave power is available it would be possible to use three 
cavities instead of C, and C,, two of which could be tuned to the same fre- 
quency and one to double this frequency. By means of a switch one could 
then select the two cavities which would be suited best for a particular expe- 
riment. For this reason in Fig. 1 the combination C,, C,, C3; is drawn. 


2. — Analysis of operation. 
Consider two high Q transmission cavities ©, and C, having resonance 


frequencies v, and »,. For C, the relative microwave power transmission is 
given by (°) 


a paci eo 


(9) C. G. Montgomery, R. H. Diske and E. M. PurcELL: M.I.T. Lab. Series, 
Vol. 8. 
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the same holding for C, replacing Q, and Q, and », by »,. If we assume that 
both cavities are followed by a matched detector, not necessarily quadratic, 
and a linear video amplifier, and if we apply linear frequency modulation to 
the microwave generator, then the output voltages V of the amplifiers will 
be functions which are symmetryc with respect to the frequencies », or », 
at which they reach maximum values. If the frequency modulation law is 
triangular the oscillator frequency may be represented by: 


(2a) y=y+ut when O<t< T/2, 


(2b) ÿ =» hott) when Tp<t< 7, 


v =|dy/dt| being the modulation velocity and 7 the repetition time. 

If we assume that a pulse is fired each time the amplifier outputs reach 
preset fractions of their maximum output, for ©; being f, = V/V(v,), and only 
after such a maximum is passed, then for ©, these pulses are fired at the fre- 
quencies 


(3a) vi =v;(1+k,/20,), when dv/dt>0, 
(3b) vy, =v(1— k,/2Q,), when dr/di<0, 


which are symmetric with respect to the resonance frequency »,. The constant 

k, depends on f, and the detector characteristics. 

Ni Similar relations hold for the second channel (C:), 
—V>0 


us with as a rule kj ~k, and Q, <Q. 


di res As explained in Sect. 1 the pulses from chan- 
ae foe nel 2 trigger the sweep of a synchroscope. The 
0 — — — — +———v<0 velocity of the spot on the screen, which is very 
——"- — high, we call w. The pulses from channel 1, delayed 
: à by and adjustable time 7, produce deflections along 
Fig. 2. — Picture on the È 
the y-axis. One therefore observes on the screen 
screen of the synchroscope : : 
showing the various distan- two pulses of which the distances from the start- 
ces occurring in theformulae. ing point of the sweep (0 in Fig. 2) will be 


w | kor. kyw a dy 

(4a) d+ = wt a 0, a a (v, — vi) when el 
w | Kav kw a dy 

(4b) d- = wer DI Fe | ” (vv) when "E #0. 


‘ so that the distance between the observed pulses 


Sg ae | 
(5) d,=|dt—d = 
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is independent of the constants 7, k, and k, and the Q-values of the cavities. 
On the other hand the distance from the starting point 0 of the sweep to the 
« barycentre » of the observed pulses is 


(dde i w 
(6) (sn z| 


kos pal 
O Gt 


If a gas is let into one or into both cavities, by means of formula (5) one may 
obtain all quantities related to the refraction of this gas, whereas from (6) 
one may deduce the absorption. In order to elucidate this we will discuss 
some examples of practical interest. 


2°1. Refraction and absorption. — The empty cavities are tuned to the same 
frequency v = v, = »,, then pulses appear superimposed on the CRT screen. 
If a gas is admitted into C, the distances d, and d, according to (5) and (6) 
will change by the amounts 


2w 


(7) op = a Oto 
è a) ei! 
(8) dò, = op ko È a. ; 


in which dy is the frequency shift due to the admission of gas and Q is the 
quality factor of the gas-filled cavity C,: Now dy is related to the real part 
of the refractive index n' by (°) 


il 
(9) DEN (1 — =) ww vo(n' —1), 


whereas the power absorption coefficient is given by (°) 


2709 | 1 I 
(10) 2a" = Fe z| | 
It therefore follows that 


AMI E 
(11) n Su 


\* TO ee iy, ‘sl TARA A TIEE 
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se that. frem the reading ef È, and à, immediately follow the values of 1’ 


amd Dx af this frequemer m. Measurements of this kind will be discussed in 


à RESF Paper 


TI Srong Gsversien. — By superimposing the observed pulses the empty 
cavia may alse be tuned to frequencies which differ a fixed and accurately 
Emewnm amount — x and + n respectively (of some MHz only) from a fre- 
Quenex »,. If this d done and a gas à admittad into both cavities simulta- 
memesty. them the change in position of the pulses may be used to measure the 
person D— \e'[2Ax at frequency n. 

Tihs may be dome as follows. The triangular modulation voltage is decreased 
se that the frequemery band covered by the frequency modulated klystron is 
è Rick las than 2. If new €, is tuned to a frequency just below this 
Rand and C, te à frequemey just above it no pulses will be seen on the synehro- 
scepe sereem. After this, in addition to the low frequency triangular voltage, 
è Smusendal radio wave of frequency An from a quartz crystal oscillator is 
Td te the repeller af the kiystron. Apart from microwaves of frequency v 
the Rirstroa will then also generate power at the side bands (» + mAw,), in 
which m & am integer. The amplitude of the sine-wave is now adjusted to such 
è value that the larger part ot the klystron power is fed into the first pair of 
Side Bands. Le. for m= 1 at frequencies r-bAr. Each cavity will then re- 
sonale on ome of these side bands and the corresponding pulses will be seen 
om the sereen. If these pulses are superimposed then, the cavities will be tuned 
te frequencies which differ exactly 24» in frequency. If now a gas is ad- 
mitted into beth cavities, having refractive index n; and à, at the resonance 
frequencies vy, — x, — An and à = »+ Ar respectively, then the observed se- 

aration of the pulses will be 


7 de | l À l l 
as ue) 
el \N à, mn, à, 
er puttine n, —R, = In 
don 4 
(it d. = — [n(n — #) + 2Ar (a — 1)]. 
n 


Ti (® — 1) at frequency n is very small or is known from a previous measu- 
rement the dispersion follows immediately from (14). 


23. Weak dispersion. — The empty cavities C, and € are tuned to », and 
mn, respectively, x being integer. The microwave power at a frequency nr 
is obtained by placing an adjusted crystal multiplier between klystron and 
eavity. The measurements at a frequency mi, as listed under 2°1, may now be 


= } CT ii Pt | 
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repeated for the frequency nv, by admitting the gas into C,. One may apply 
the formulae (11) and (12) if the right hand side of (12) is multiplied by n. 

If the gas shows little absorption in the frequency region from », to m 
the relative variation of the refractive index is very small and difficult to 
detect if the refractive index is measured directly at the different frequencies. | 
Such measurements can; however, be carried out easily if the gas, at the same 
pressure and temperature is admitted into both cavities simultancously. The 
quantity w’ in (11) then denotes the change in refractive index over a frequency 
region (nv, — vg). Therefore a set-up like this is perfectly suited for the measu- 
rement of a small dispersion over à large frequency band. 


24. Increased sensitivity for very small absorption. — Up till now in our 
discussion the synchroscope sweep was always initiated by the pulses from 
channel 2, whereas the observed pulses came from channel 1. If alternatively 
for one period of the modulation one would trigger the sweep with the pulses 
from channel 2 and the next period one would 


trigger with the pulses from channel 1, instead qf gg 2e 

of one pulse (Fig. 2) on each of the traces on da E e 

the screen now two pulses would be seen. In era eg 
order to achieve this result an electronic switch TARA | PRES 

should be inserted between the discriminators g-— ———~L---«0 
and the pulser outputs. If both cavities are tuned Figs 3 + Pica dentine 
to the same frequency v,, and k,=k, and Q,=%; screen of the synchroscope 
then dj =d;, d;=d, and A+=A"=0 (see Fig. 3). if an auxiliary electronic 

If we let gas into one of the cavities and then switch is applied. 


tune this cavity to have the pulses on one trace 

cointident again, then it may be easily shown that the distance between 
the pulses on the other trace is four times as large as the distance 6, of 
formula (8). So by adding an electronic switch as indicated the sensitivity for 
absorption may be increased by a factor four. 


3. — Sensitivity of the apparatus. 


The sensitivity of the apparatus is limited by the low frequency noise pro- 
duced by the crystals detectors and multipliers. Because of this noise there 
are fluctuations at the output of the amplifiers, causing irregularities in the 
firing point of each discriminator, #0 that the pulses on the screen show a 
jitter. 

We have for the indetermination Ay in the discriminated frequencies vt, y 
produced by this jitter 
av dP , 


ns 2 ee ri ri 
oily ie Doi 


(15) AV 
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V and N being the signal and noise voltages at the outputs of the amplifiers. 

If the maximum power transmitted by the cavities is of the order of 10-5 W 
then Av can be easily calculated. In this case the noise produced by the 
crystal is essentially the Johnson noise of a resistence equal in value to the 
d.c. impedence of the crystal, and the crystal operates as a quadratic detector. 
We have, from (1) and (15) (?°) 


Av 4V3 VARTE 
jee SOP MO 


B and M being the bandwidth and the overall figure of merit of the video 
receiver, and dP/dv being calculated for frequencies »*,»- at which V/V,,. = 
= 3/4 (k, = 1/V3). 

With P, 10-56 W, M = 100, Q= 30000, we deduce, for the minimum de- 
tectable variation of refractive index: 


Ly 


(= 1)me = (5) Rake? 
Vo /min 


At present, the sensitivity of the apparatus is of an order of magnitude 


smaller, being 
(2) a UO, 
Vo min 


The sensitivity in the measurement of absorption is proportional to that 
of refraction, and depends from frequency, because from frequency depends 
the physical phenomenon from which absorption is deduced. Indeed from (12) 
we deduce. 


Agcy, (Av 
2a" min = . 
a ck, | Vo ti 


* OK OK 
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(19) W. Gorpy, W. V. SmirH and R. F. TRAMBARULO: Microwave Spectroscopy 
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RIASSUNTO 


In questa nota viene descritto un apparecchio per misure del coefficiente di assor- 
bimento e dell’indice di rifrazione di gas a pressioni relativamente elevate, costruito 
nell’Istituto di Fisica dell’Università di Pisa. Confrontando mediante una tecnica ad 
impulsi due cavità risonanti opportunamente accordate, viene ottenuta una sensibi- 
lità di 10-8 nelle misure dell’indice di rifrazione e di 10-5 m-* nelle misure del coeffi- 
ciente di assorbimento. 
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* The Electromagnetic Cascade with Tridents. 


J. W. GARDNER 


Atomic Power Division, English Electric Co., Ltd. 
Whetstone, Leicester, England 


(ricevuto il 16 Luglio 1957) 


Summary. — Recent results claiming anomalously large cross-sections 
for direct pair creation by high energy electrons (mainly above 10 GeV) 
are discussed with particular reference to their significance in the deve- 
lopment of the electron-photon cascade. A qualitative estimate, using 
an admittedly crude model, suggests that the effect of such large « tri- 
dent » cross-sections could drastically change the pattern of development 
in the early stages of a really high energy cascade. The detailed analytical 
formulation of cascade theory is then modified in such a way that nu- 
merical values of the cross-section for trident production, when these 
are satisfactorily established, can be readily inserted in the formulae 
used for computing the moments of the number distribution function. 


1. — Introduction. 


The so-called trident process, i.e. direct creation of electron pairs by inter- 
actions of fast charged particles with Coulomb fields, has been studied both 
experimentally and theoretically by a number of workers (1). Until recently 
it appeared safe to assume from the available evidence that trident production 
by electrons, even at the highest cosmic ray energies, was a much less common 
occurrence than emission of bremsstrahlung or photon materialization. In 
view of this, and the complexity of cascade theory even without tridents, there 
has been an understandable reluctance to include terms representing trident 
production in the diffusion equations describing the multiplication of particles 
in the electromagnetic cascade—notwithstanding an awareness that even a 


(1) Extensive references are given in a paper by M. M. BLocK et al.: Phys. Rev. 
46, 1627 (1954). 
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few tridents might have a non-negligible effect on the initial development of 
a cascade with energy above 10 GeV, say. 

In the past few years, however, the justification for neglecting the trident 
terms in cascade theory has been somewhat shaken by observations (2) and 
calculation (75) which appear to indicate a distinctly larger cross-section for 
trident production, especially above 10 GeV, than would have seemed pos- 
sible from earlier studies. Not all these recent reports have gone unchallenged; 
for example, BLock and KING (*) have criticised KosHIBA and KAPLON'S 
interpretation of their data. [A basic difficulty in the experimental study of 
tridents arises from the tendency to include in the « trident count » a fraction 
of those pairs which are due to conversion of bremsstrahlung photons so close 
to the track of the parent electron that the event cannot be distinguished 
observationally from direct pair creation by the parent electron. There ap- 
pears to be some disagreement among experimentalists as to the allowance 
to be made for such pseudotridents; also for «spurious scattering » and other 
effects leading to under-estimation of energy.] Nevertheless the persistence 
of reports of « anomalously » large trident production, coming from a number 
of independent groups all over the world, does commend a re-assessment of 
high energy cascade theory to include a significant trident component. 

In the present paper we shall first consider qualitatively, in Seet. 2, how 
the electromagnetic cascade would develop if the mean free path for trident 
production were equal to that for bremsstrahlung emission, which at energies 
above 100 GeV, is about the value indicated by references (2-4) (*). 

Even allowing for the crudeness of the model used, the results are so striking 
that it has seemed worth while to extend the quantitative formulation of 
MESSEL and Ports (!°) to include tridents. This is done in Sect. 3, where we 
set up the matrix equations for the number distribution function and its 
moments, and compare them with the corresponding equations when tridents 


(2) J. E. NAuUGzE and P. S. FREIER: Phys. Rev., 92, 1086 (1953). 

(3) J. E. NauGLe and P. $. FREIER: Report of Duke University Cosmic Ray Con- 
ference (1953). 

(4) M. Kosurpa and M. F. KapLon: Phys. Rev., 95, 647 (1954); 97, 193 (1955); 
100, 327 (1955). 

(5) R. WerLL et al.: Helv. Phys. Acta, 29, 437 (1956). 

(8) A. DEBENEDETTI et al.: Nuovo Cimento, 12, 954 (1954). 

(?) T. Murora et al.: Prog. Theor. Phys., 16, 482 (1956). 

(8) J. W. GARDNER: Sydney University Research Report, 1956 (unpublished). 

() M. M. Brock and D. T. Kina: Phys. Rev., 95, 171 (1954). 
(*) Extrapolation of the conventional theory of trident production by electrons 
(see e.g. reference (1)) gives a mean free path in emulsion (7 — 40) of 8 radiation 
lengths at 100 GeV; however, such extrapolation is shown to be unjustified by MUROTA 
et al. (7) who by a more precise calculation obtain a m.f.p. of 3 r.l. at this energy. 

(0) H. Messe and R. B. Ports: Phys. Rev., 86, 847 (1952). 
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are neglected. Numerical solution of these equations is not attempted here; 
while the time and labour involved would of course be greater than for the 
corresponding calculation without tridents it should not be prohitibive, given 
a fast electronic computer, to obtain numerical results for the first few mo- 
ments and use them to construct the distribution function. However, one is 
deterred from embarking on such a computation until firmer evidence is 
available on the numerical values to be inserted for the trident production 
cross-section. It seems likely that further experiments in the course of the 
next year or so will supply this evidence. 


2. — Qualitative considerations. 


A qualitative model for the electromagnetic cascade without tridents has 
been described by HEITLER (11). While making no claim to detailed accuracy, 
the model does give a reasonable picture of the general development of the 
cascade. In the present section we briefly recapitulate this model and then 
discuss how it may be modified to include tridents. Some rough estimates of 
the average numbers of electrons are made, taking the mean free path for tri- 
dent production to be equal to that for bremsstrahlung emission, and using two 
alternate assumptions about energy partition among the three trident electrons. 

In the Heitler model it is assumed that whenever an electron traverses. 
a distance In 2 radiation length (*)—the distance in which its energy has on 
the average dropped by half—it emits just one photon of half its energy and 
retains the other half. Likewise, a photon in distance In 2 is assumed to 
undergo materialization with the photon’s energy shared equally between the 
two electrons. We ignore the fact that the mean free path for photons slightly 
exceeds that for electrons; both are taken at the high energy limit, where 
Compton scattering and collision losses may also be ignored. The model readily 
leads to the following results: 


i) The number of particles in the cascade increases exponentially, 
within limits. 

ii) For not too small depths the number of electrons in the cascade is 
about double the number of photons. 


iii) The depth, ¢,,,., at which the number of particles in the cascade is 
a maximum is given by #4, ~ In (H,/KE.), with E = primary 
energy, E.= critical energy, K = a constant of order unity. 


iv) The number of electrons at the maximum & 2£/3KE.. 


(1) W. HEITLER: Quantum Theory of Radiation (Oxford, 1954), 3rd ed., Section 38. 
(*) Not one radiation length, as erroneously stated in B. Rossrs: High Energy 
Particles (New York, 1952), Section 5'3. 
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To incorporate tridents into the above model let us assume that we have 
not one, but a large number N, of primary particles (electrons and photons) 
all of the same energy incident at t= 0; also in this case it is more con- 
venient to divide the depth into layers of thickness In 1.5 radiation length 
instead of In 2.1 r.l. The solutions for the Heitler model [equations (7) and 
(8) below] predict that in the first In 1.5 r.l. half of the primary electrons will 
emit one photon each, and half of the primary photons will undergo mate- 
rialization. Now our new model postulates that the bremsstrahlung and 
photon-materialization mean free paths are unchanged but that, in addition, 
trident production by electrons occurs with the same mean free path. This, 
incidentally implies the assumption that the mean free path for trident pro- 
duction is independent of energy in the high energy limit. A further assump- 
tion we make is that the probabilities for a given electron to emit a brems- 
strahlung and to make a trident in 1.5 r.l. are additive. This greatly sim- 
plifies the argument and is consistent with the general level of approximation 
of the Heitler model. On the present model, then, half of the primary electrons 
each produce a trident and half each emit a photon in the first In 1.5 r.l.; 
the primary photons behave just as on the previous model, half undergoing 
materialization and half continuing unchanged. Thus in the first In 1.5 r.l. 
all the primary electrons and half the primary photons undergo some mul- 
tiplicative transformation. What happens at subsequent depths is determined 
by our assumptions about the partition of energy among the three electrons 
of a trident. We shall consider two opposite assumptions: a) each electron 
receives about one third of the total energy and is therefore itself capable of 
producing further tridents; b) one electron, the « primary », receives the bulk 
of energy and is the only one capable of producing further tridents, although 
the two secondaries may still participate in the cascade development by emit- 
ting bremsstrahlung. The truth lies somewhere between these two extremes, 
but probably much nearer to b) and a), see, for example, reference (1). 

Let us denote the numbers of electrons and photons at depth t by A(t), 
B(t) respectively and also write x for In 1.5. Then the equations describing 
the multiplication of particles in the cascade are: 


(1) 2B(t + x) = Bt) + A), 
(2a) AGE BOAT, (Assumption a)), 
(2b) A(t +a) = B(t) + A(t) + A(0), (Assumption b)). 


° Equations (1) and (2a) taken together are satisfied by: 


(3) A(t) & 0.78N, exp [2.034], 
(4) B(t) & 0.22 N, exp [2.031]. 
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(1) and (25) taken together yield: 

(5) A(t) & 1.27 N;fexp|[1.42t]— 0.39}, 
(6) B(t) & 0.36 N, {exp [1.42] — 0.391. 


For completeness we also recall the corresponding functions for Heitler’s above 
model which ignores tridents: 

(7) A(t) & 0.67 exp[t], 

(8) B(t) = 0.33 exp [t]. 

To obtain t,., for the cascades represented by the above models one uses. 
the criterion that the average energy per particle (electron or photon) at the 
maximum is of order KH,. One finds then, that on either assumption a) or 
assumption b), the inclusion of the trident process causes the maximum to 
be reached in about half the distance that is required if tridents are ignored. 
To take a specific example, consider a cascade in air initiated by a 100 GeV 
primary: ignoring tridents t,,,, > 6 r.l.; with tridents under assumption a) 
tnx © 3 r.l., and under assumption b) t,.,, & 3.7 rl. 

Equations (3).to (8) further show that tridents, besides giving a much 
faster initial multiplication of particles, also increase the ratio of electrons to: 
photons: without tridents it is 2:1, with tridents (on either assumption) it is 
just over 3.5:1. The effect on the multiplication of electrons of including 
the trident process in the cascade is illustrated numerically in Table I, and 
appears particularly striking, even with due allowance for the crudeness of 
the models used. The qualitative results of this section thus suggest the need 


TABLE I. — Average numbers of electrons from a single incident primary. 
A,(t) = 0.67 exp [t] - tridents ignored, 
A,(t) = 0.78 exp [2.03t] - tridents included, assumption (a), 


A,(t) = 1.27{exp [1.42t] — 0.39} - tridents included, assumption (b). 


ib È radiation ie | A; A; Agen, As/ Ay 
engths 

2 | 4.95 45.3 DIR, 9.15 4.29 

È HS 343.0 89.5 25.4 6.55 
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for having available a more quantitative formulation of cascade theory, into 
which numerical values for trident production cross-sections can be substituted 
when these are more reliably known. Such a formulation is outlined in the 
next section. 


3. — Quantitative formulation. 


In this section we follow closely the method and notation of MESSEL and 
Ports in reference (1), denoted hereinafter by MP. As is done in that paper 
we ignore scattering both elastic and inelastic, i.e. we neglect both ionization 
loss and the lateral spread of the cascade. These approximations are certainly 
good ones at energies where trident production is expected to be significant. 
The only other assumption we make is that the cross-section for trident pro- 
duction is, like those for bremsstrahlung and photon materialization, a homo- 
geneous function of degree zero in the primary and secondary energies. More 
specifically, we assume that the probability for a primary electron of energy Ep. 
to make in distance dx a trident with secondary energies in the intervals dH; 
at H;, dH, at E, and dE, at E, can be written w®(y;, 4%, 91) An: nr Am dx with 
n= E;/Eo, Nr= Ey/Eo, M = EE, and, by energy conservation, 7; + 4.+ M = 1. 
The total cross-section for trident production is then written: 


ileal 


ge) — [Jor Nr» M) AH. AM à 


00 


For convenience we recapitulate briefly the more important definitions 
of MP. The differential cross-sections for bremsstrahlung and for photon 
materialization (Bethe-Heitler full screening approximation) are respectively 


OÙ (xs M) Anr Ande and (7x, 71) Ar di dx ; 


the corresponding total cross-sections are 
1 1 
an = | OMA, mdm and «= | OM Hr) dm - 
0 0 


Phe function po, (mi; +1 Mas Mr Nn ms 2) OM, day 18 Gehined-as) the 

differential probability that after a depth x a primary (j) of unit energy has 
given rise to » electrons with energies in the ranges 7;, 4.+4y,, with k= 1,..., ”, 
in any order, and m photons with energies in the ranges %n41, Mnt + Matty 
with /=1,...,m in any order. Further, the function 


(i) 


drm(M13 DUC; UE) Mn + 19 Ce) Mme x) dy, ss de 
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is defined as the differential probability that after a depth x a primary (j) of 
unit energy has given rise to n electrons and m photons with energies specified 
as above, plus any numbers of electrons and photons with arbitrary energies, 
including if we wish the same energies specified above. A relation between gia 
and pl, is given in MP; hence if either q,, or p,,, is known the other may 
in principle be determined. However, it is found more expedient to derive 
diffusion equations for each of q\’,, and p?,, and solve them directly. 

The usefulness of the quantity q\”,, called by some authors the « product 
density », lies in its very simple relation with the moments of the practically 
useful function pym(7, ©) expressing the probability that a primary (j) of unit 
energy has given rise, after a depth x, to n electrons and m photons with 
energies > 7 and any numbers of electrons and photons with energies < 7. 


This relation is (12): 


(3.1) TOG: x) == fan se fan gio <009 Mn 3 Yn+19 SO} Yntm > x) ’ 
n 


where 7” (74; x) is the (n, m)-th factorial moment defined by: 


© © 1 (; ! 
(3.2) Ds 2) = SEE uma). 
Since the enormous complexity of the analytical expression for y”, makes 
its direct calculation quite prohibitive, even with an electronic computer (1), 
the value of q\”, for obtaining the moments is immediately apparent; for one 
can get certain information about the shape of the distribution from its first 
few moments and may in fact be able to construct a reasonable approximation 
to gp”, from an expansion in the moments. 
We now proceed to set up the matrix recurrence relations giving the p”, 


() 
nm ° 


and q 


3'1. Equation for p;,,. — The diffusion equation satisfied by the p}, is 


n,m * 


obtained by considering all possible last collisions, just as in MP. It is: 


{0/0x + n( (az) eee MA OS Bee 
ae DI DI DIA nn ALP À De Mn +Mns mò ee SCT) Des 20M, Naam) sa 


n m 
CIAO) 


de > De 2 i of) ae Mn +17 RE] Mn + m 3 inn Sig ie Ma M) ae 


or > pe 2,m (Chee Oy) ee ieee + Mi +7, Mn+ 1 ne) (ne O N) ’ 


(2) H. MessEL and R. B. Ports: Proc. Phys. Soc. (London), 65, 854 (1952). 
(5) H. MesseL and J. W. GARDNER: Phys. Rev., 84, 1256 (1951); 86, 808 (1952). 


Ce 
nm 
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with the initial conditions: 


(3.4) DIA TR 0) GE Oo d(1 ane Un) o 

In (3.3) the Cf and C} signify that the summations are taken over all 
possible choices of 7, and 7, ,, n, respectively from the n,, k=1,..., n; the 
C? signifies summation over all possible choices of 7! > TOI eni 
= 1,..., m; O% signifies summation over all possible choices of 9, ,, M, 1 9, from 
the x, k=1,..., n. The only differences between our equation (3.3) and the 
corresponding equation (eqn. (8)) of MP are (i) the coefficient of n on the 
left side is now x +x instead of just «®; and (ii) the right side now includes 
a third term representing trident production as well as the original two terms 
representing bremsstrahlung emission and photon materialization. The initial 
conditions are exactly as in MP; moreover, since any cascade with tridents 
cannot contain less than three electrons, the right side of (3.3) must revert 
to the right side of eqn. (8) of PM for n < 3. 

If we define the Laplace transform of p°,, as: 


(co) 


(3.5) | Pate | exp [— Ax] pin(a) da, 


0 


then (3.3) and (3.4) may be transformed into the following matrix equations: 


(3.6) [AE,+ > a,(k)| Plt.) ---9 My; À) = 


et) 


| 


da ORA UTO DIEU <> Yy—o9 Mn DE Ny A) ci 


0% 
pre yes Mo ae pio 
C3 
and 


(3.7) [AE, ae a (1)]Pi(17,, A) = E, d(i — M) . 


In these equations E, is the unit matrix of order 2". The a,(k) are given 
by the direct product of N matrices: 


To Aa DA) 0 
(3.8) PEC PIO 7 ” CRE CO JA 
k-th place 


The P, is as 2" x 2 matrix the columns of which correspond to P® and P® and 
the rows are ordered by writing 7, ...7y as a binary number with digits 1 and 2, 
standing for an electron and photon respectively. The w,_, is a 2x2" 
matrix in which the only non-zero elements are ©°2(7:,7,) and (73,1) 
ordered according to certain rules. The v, , is a 2"x2" ® matrix in which 


= 2 - Il Nuovo Cimento. 
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the only non-zero elements are OO (M, Mrs M) ordered according to certain 
rules. The rules for ordering the non-zero elements in w, , and v,_, can be 
readily formulated by considering a specific case, say N —3. MP write them 
out expressly for w,_,). 

Comparison of (3.6) with the corresponding equation (eqn. (11)) of MP 
Shows that the introduction of tridents into the cascade replaces the simple 
matrix recurrence relation between P, and P, , by a relation between P,, 
Pee ean Ee os: 

32. Equation for q,,. — The «last collision » diffusion equation satisfied 
by q®,, for the cascade without tridents (eqn. (15) of MP) contains five sum- 
mation terms on the right side, corresponding to the different types of last 
collision involving bremsstrahlung emission or photon materialization. When 
tridents are included three additional terms arise corresponding to the pro- 
babilities that the last collision is a trident production with one, two or three 
of the resulting electron energies falling in the specified intervals. As in the 
case of pl’, the only change to the left side of the diffusion equation is that n 
is multiplied now by «+ instead of just «”. In writing out the new 
diffusion equation below we omit, for brevity, the arguments of all q’s on 
the right side except the three arising from tridents (*), and in these latter 
we omit the 7,41; ---1 Mnim representing photon energies which of course are 
all unchanged by the last collision if this is a trident. In this abridged notation 
the equation becomes: 

(3.9) {ofexe + nla? + a) + ma Mar 203 Mes noid ce. 


> Mn+m) x) i 


— 2 D Qn Wo 10(Nny Tinea) a dn Je 2, miO(Na a, Nn) # 


Ci È; 
1 


ak di Ro Nas Mn) ) du + Za [Eno u— Mn) du + 


0 
1 


(3) ! I 
Ain: Da | Un, m — 1#00(% TRE Yntm ’ Nain) du Sia 
Ci 


(Î) if al ! Ne if 1 I 
ae Sa dala LOL) Yn-3 9 Vn-2 a Mn 3 x) OSS, Mn=1 Nn) ala 


C3 


= (9) ! ! ne ! I ali i 

ate >. | m1 m(Mns ces Mnr2s U} CIO (na) Mn U— Nn—1 — Mn) AU + 
Gh a 

0 
al 


if I ti 
+ 2 | if n (11, M1 U5 DIO (Mn, Vs U — N — 0) AV AY . 


U_U 


(*) The arguments of the other five q’s are exhibited in MP eqn. (15). 
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ì . sz (j . . . . . 
Noting that qy,(r)=1 for all x, our initial conditions are: 


(3.10) die) =. d(1 —7,); D (C—O 


M+-j,2 00% 


nm 


We define the following Mellin-Laplace transform of q\”,: 


a) 5 5 
(3.11) Qnm(819 +++) 8n3 Sat ++) Sntm3 A) = 
1 1 (co) 
= | dj... Jar | de ni... niin exp [— Ar] grim 5 
U 0 0 


together with the quantities: 
4 
; Mai È a a eee 
3.12 1140) So ——— D D (: 5) dn, 
(3.12) es ETC md ane 
0 


foro d'or 2 and 


(3-13) W(s,, $$) = 


Si No Sg 7 Sg | 
= =A || head ( Baile =» | | Een ) AM) Noy Na) An, Ans + 
Na + Me + Ns Nt M + Ns Mi + Ms + Ns 


Operating on (3.9) with (3.11), using (3.10), (3.12) and (3.13), and sup- 
pressing for brevity the arguments of the Q’s on the right side, we have the 
following. transformed equation: 


(3.14) {A+ n(x DL a) + MAP?” (8, #09 Sn3 Sn+19 +9 Spam > À) — Onur, Ont = 
= Da Zim Oum Wu Siem) + Bn Qu amp WP (En as 89) + 
(iin A CA 
ARR ime 0) 2 OMO Coe tee 
0; 
# >> Qua wo, Se 13 2 Oy UAC SE S,) “6 
CA 


ate Di QI of 8,30 = 2. ADS AE CE 7,9 0, 0) © 
Co 


- Again denoting n+m by N we see that each of the Q’s in (3.14) refers to 
either N, N —1, N — 2 particles; the matrix form of (3.12) is then a recur- 
rence relation between Q,, @,_; and Q, ,, which are matrices constructed 
from the Q’s exactly as the P, in (3.6) are constructed from the P’s. The 


a 
ri 
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recurrence formula is: (for N > 1) 
| 
(3.15) [AE, + 25 BIO CL Se 
= 2» Wels 558%) OF SS recy Spa oe oe A) A 
a 2» Vig AOS ee OC eS EE Aloe SA) 


The B,(s,) are given by the direct product of N matrices: 


B;(83) B3(8%) 


(3.16) Bi) = Ba rs 
ae B.(6)  B(s,) “ 
k-th place 
with 
B,(s) = xD + at — Ws, 0) — 3W®(s 0, 0), 
B,(s) = — 2Ws, 0), 
(S217) x 
B;(s) aa te VEO, 8) ’ 
Bia E 


The only respect in which B, differs from the corresponding matrix A, for 
the cascade without tridents (MP eqn. (22)-(23)) is in the inclusion of the x® 
and W® terms in B,(s). 

The W,_, are 2"x2Y7! matrices where the only non-zero elements are 
comprised of W(s,,s,), W(s,,8,) or 3W(s,,8,,0) ordered according to 
certain rules, readily formulated by analogy with the w,_,(#y_1:7») of (3.6); 
the V,_, are 2"x2"-? matrices and are just the v,_, of (3.6) with w®(n, 7,, 7,) 
replaced by W®(s;, 84, 81). 

For N=1 the recurrence relation assumes the quite simple form: 


(3.18) [AE, + B;(s)] Q:(s, 4) + Ey. 


4. — Discussion. 


The relative usefulness of p’,, and gq, is discussed in MP. The solution 
for Q, is the one of direct physical interest because it leads, via the relation (3.1) 
to the general solution for all the moments, including the correlations between 
electrons and photons. As an illustration we write out the solution of N — 1 
explicitly. Equation (3.18) yields, on performing the matrix multiplications, 


a 


THE ELECTROMAGNETIC CASCADE. WITH TRIDENTS 21 


the following four simultaneous equations: 


(4.1) As BOG are On =1, 
(4.2) in (BOVA ABQ; = 05 
(4.3) (AB) 000.0 = 0, 
(4.4) BOO GO BOOT. 


Let us suppose, for example, that we wish to calculate the average number 
of electrons at depth x and above energy y arising from a single primary 
electron. For this we require Q‘\?, which from (4.1) and (4.2) is found to be: 
À + B,(s) 


4.5 Qu — . n TRE 
ca on PPDA EI SOC 


The inverse Mellin-Laplace transform of @Q{? integrated over energy from 7 
to 1 then gives the average number directly: 


Ut? lot io 
de 
Gy = (1) 
(4.6) IZ Cai SE fu exp [Ax]n- Qi, - 
Ug — - io dn to 


To obtain numerical results from (4.5) and (4.6) one must know the values 
of the cross-sections occurring in the equations (3.17) which define the B’s. 
Only B,(s) involves the cross-sections for trident production; the other B’s 
are already known from the Bethe-Heitler full screening cross-sections for 
bremsstrahlung emission and photon materialization. Computations of average 
numbers and correlations using these cross-sections, and ignoring tridents, 
have already been performed on desk machines, using a saddle point ap- 
proximation for the integrations (14). To this approximation it should also 
be feasible to include the effect of tridents in a desk machine calculation of 
average numbers, given reliable cross-sections to substitute in B,(s). Certainly 
with an electronic computer the inclusion of tridents in average number cal- 
culations should present no great difficulty, whether the saddle point ap- 
proximation or direct numerical integration be employed. To obtain the 
higher moments, the electronic computer could be programmed to solve (3.15) 
for Q,, given Q, , and Q, _,, starting with N= 2. In practice it may not be 


(14) L. JANossy and H. MessEL: Proc. Phys. Soc., A 68, 1101 (1950). 
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necessary or desirable to go beyond N —2 or 3 in order to obtain enough 
information about the distribution function for it to be usefully employed 
in the interpretation of experimental results. 

We have not attempted to give a general analytical solution for Q,, cor- 
responding to MP’s equation (24); such a solution may be obtained by analogy 
with their procedure but would be even more unwieldy than that of MP. 
The best hope for obtaining practically useful results appears to be direct 
numerical solution of the recurrence equation along the lines indicated above. 


Most of the work reported here was performed during the tenure of a re- 
search fellowship at the University of Sydney. 


RIASSUNTO. (*) 


Si discutono con particolare riferimento al loro significato nello sviluppo della 
cascata elettronica fotonica alcuni risultati recenti che segnalano sezioni d’urto anor- 
malmente grandi per la creazione diretta di coppie da parte di elettroni di alta energia 
(principalmente al disopra di 10 GeV). Un esame qualitativo, con l’ausilio di un modello 
evidentemente grossolano, suggerisce che l’effetto di tali grandi sezioni d’urto di «tri- 
denti» potrebbe, decisamente cambiare lo schema di sviluppo nei primi stadi di una 
cascata di energia realmente alta. Si modifica quindi la formulazione dettagliata della 
teoria della cascata in modo da poter senz’altro inserire nelle formule usate per il cal- 
colo dei momenti della funzione di distribuzione numerica i valori numerici della sezione 
d’urto per la produzione di tridenti, quando questi siano stati convenientemente ac- 
certati. 


(*) Traduzione a cura della Redazione. 
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Muitiple Scattering Measurements in Nuclear Emulsions. 


S. ROSENDORFF and Y. EISENBERG 


Department of Physics, The Weizmann Institute of Science 
Rehovoth, Israel 


(ricevuto il 18 Agosto 1957) 


Summary. — The practical use of the new multiple scattering function, 
cos ny, introduced by LIPKIN (!) et al., was examined by comparing its 
theory with experiment. At the same time the mean absolute angle ac- 
cording to Moliére’s scattering theory was compared with experiment for 
the second and third difference and the respective dispersions were de- 
termined experimentally. (The theoretical dispersion diverges in Mo- 
lière’s theory.) A good agreement between theory and experiment was 
obtained in the cosine case, both for the mean value and the dispersions, 
as well as in the mean absolute angle case with cut-off. Writing the 
relative error of pf as D,,/4/n (for the second and third differences, 
respectively) we have obtained in the cosine case: D,=0.976 and D,=1.36, 
and in the mean absolute angle case (with cut-off at 4<|p|>): D,=0.97 
and D,=1.37. It is shown that, provided one restricts oneself to a cer- 
tain class of estimates of pf, the above dispersion values are close to 
the minimum possible dispersion of pf. Also, an explanation is suggested 
for the known experimental fact that the ratio of the « scattering con- 
stants » K,/K, is higher by a few percent than the expected value 4/3. 


1. — Introduction. 


Multiple Coulomb scattering measurements of charged particles in nuclear 
emulsions is very often used in order to estimate the product pf of the particle 
momentum and its velocity. Several detailed multiple scattering theories (?) 


(1) H. J. Liexin, S. ROSENDORFF and G. YEKUTIELI: Nuovo Cimento, 2, 1015 (1955). 

(2) G. MoLIÈRE: Zeits. f. Naturf., 2a, 133 (1947); 3a, 78 (1948); 10a, 177 (1955); 
H. Snyper and W. T. Scorr: Phys. Rev., 76, 220 (1949); S. A. GOUDSMIDT and J. L. 
SAUNDERSON: Phys. Rev., 57, 24; 58, 36 (1940); H. W. Lewis: Phys. Rev., 78, 526 
(1950). The present work will refer mostly to Moliére’s calculations. 
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have been developed, by which it was possible to calculate the value of pf 
from the scattering angles of a particle in matter. In nuclear emulsion practice, 
‘ it is common to use the mean absolute value of the projected angle y; between 
successive chords along the particle’s trajectory for determining the value of pf 
(a summary of the theory and of the experimental situation as of 1953 is 
given by GOTTSTEIN in HEISENBERG’s book (*)). The angles y, are usually 
determined by measuring the y-co-ordinates at successive equal intervals 
(equal « cell size »), of a track aligned to be parallel to the x-motion of a micro- 
scope, and by taking the second difference of these co-ordinates (this method 
is known as the « Fowler Co-ordinate Method »). Often, in order to eliminate 
possible distortions in emulsions, third and higher differences are used. 

Recently it has been suggested by LIPKIN et al. (1) that the value of pf 
could be estimated by using the mean of the cosines of angles proportional 
to the scattering angles, namely <cos 7p;> instead of the conventional use 
of <|g;|). The constant multiplier 7 is to be chosen to minimize the dispersion 
obtained for pf. It has been shown in I that by using the cosine function one 
overcomes the main difficulties which are encountered in using the mean ab- 
solute angle parameter; namely, one can calculate analytically, without diffi- 
culties, both the mean value <cos 7y;) and all higher moments. Thus, it is 
possible to obtain the theoretical dispersion of the measured value of pf in 
the cosine case. This is not possible in the mean absolute angle case due to 
the divergence of the second and higher moments. 

The purpose of this work is to check experimentally the theory of the 
cosine function as introduced in I, to compare the cosine method with the 
mean absolute value method and to determine, experimentally, the dispersion 
in the latter case. 

The experimental methods and the range-energy relation used are discussed 
in Sect. 2. A comparison between the experimentally measured and the theo- 
retical quantities is given in Sect. 3 for both the first moments and the dis- 
persions. In Sect. 4 an attempt has been made to find out whether or not 
a better estimation function for pf exists, and the conclusions are given in 
Sect. 5. 


2. — Experimental methods. 


In the experimental part of this work we have used a block of 2 in. «3 in. x 
X 400 um Ilford G5 stripped emulsions, exposed to the pion beam of the Nevis 
Cyclotron at Columbia (*). The average energy of the pions was about 60 MeV, 


(5) Kosmische Strahlung, p. 495. Edited by W. HEISENBERG (Berlin 1953). 
(*) We are grateful to Dr. G. HARRIS of Columbia University for exposing the 
plates for us. 
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and the average range about 53 mm. Thus, most of the pions stopped within 
the emulsion block and their ranges could be determined exactly. From the 
known ranges, the value of pf was determined accurately and compared with 
the values obtained from the multiple scattering measurements using the 
<cos 79;) and the <|w,|> functions. 

The range-energy relations used were derived from the measurements of 
BARKAS et al. (4). We have assumed the usual power low, namely: R—a- M1: E», 
where È is the range, Æ the kinetic energy of the particle and M its mass. 
Using the following values of pion ranges from BARKAS work: 


E,= 29.80 MeV, for k,=(152604 90) um 


E,= 50.64 MeV, for R, = (36610 +180) um 


and correcting them according to the average measured muon range (<R, 
was (602.2 + 2.5) um in BARKAS’ work and (591.0 + 4.4) um in the present ex- 
periment) we have obtained for the range-energy relations in our emulsion 
block the following expression: 


(1) HB (MeV) — 0,186: R99. 9-39 , 


where R is the particle range in microns and M its mass in units of the proton 
mass. Expression (1) is valid for pions in the range interval 9000 um < R < 
< 40000 um. All the pions we have used were in the above range interval. 

A total of 3960 angles between successive chords have been measured. 
The cell-size employed was 50 um. The measurement noise for that cell-size, 
as determined by calibration against a very energetic cosmic ray particle was: 
<dy) = (0.057 + 0.006) um. All errors quoted in this work are Standard De- 
viations. 


8. — Comparison with theory. 


3'1. The scattering constant. — The mean absolute value of the projected 
scattering angle g measured between successive tangents or chords (second 
and third differences) is given by: 


(4) W. H. Barxas, P. H. Barret, P. Cuer, H. H. Heckman, F. M. SmrH and 
H. K. Trono: Phys. Rev., 102, 583 (1956). 
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where K,,, is the «scattering constant» for the tangent, second and third 
difference respectively, and ¢ the cell size. According to the theory of Mo- 
LIBRE (2) (we have used here the notations of Molière), the scattering constant 
is a slowly varying function of 6 and t. Explicitly, it is given by the following 
expression: 


0.982 0.117 
(3) = Xu pf-V Bila (1 + ; }: 


B, B? TUO 


X,; is given for the three mentioned cases in terms of the «unit probability 
angle » per unit length (5) 7. in the following manner: 


I'D Ane4Z? 
(4) re a, ue | sites har i AI A), 


or, for G5 nuclear emulsions: y,= 0.31/pf, if t is 100 um, Z2=1 and pf is 
given in MeV/c. 

The parameters B,, B, and B, are defined by the transcendental equation: 
(5) B;,— In B; = (X;;t/X?,) — 0.154. 

The X,;s are simply related to Molière’s screening angle 7,: 
(6) AN, nX.,,=ny—-2/3: 


a2,3 nu 


The mean number of single scattering processes, 2,, occurring in the thick- 
ness ¢ of the scattering substance is given by: 


(7) 2, = xetlxa * 


The B,’s are therefore, according to (4), (5) and (6), related to 2, by the 
following expressions: 


[B-InB=InQ— 0.154 (B,= B) 
(8) ! B, ln B, = In (42,) + #— 0.154 
| B;—In B, = Q, + 3— 0.154. 


The B,’s as function of log, Q, are shown in Fig. 1. Curves for Q,/t as 
function of 6 can be found in the literature (*). 


(5) H. A. Berne: Phys. Rev., 89, 1256 (1953). 
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B, and B, are approximately equal; with this assumption the ratio K,/ Ks 
is, according to (3) and (4), given by V3. 

The exact ratio of the third to second difference scattering constants 
should be: 


6 = 


K, | /3 Bs 1 + 0.982/B, — ... 
He ET (1980) Be 


ò is, therefore, greater than V by several 
percent. Its exact value depends upon 
the value of B. 

For example, for B= 10, one obtains 
Ô—1.022V3 and for B—5 one obtains 
ò=1.046V3. In the present experiment 
ò=1.028V3. Consequently, this factor 
explains, in part, the known fact (5) that 
the experimental ratio of K;/K, is a few 
percent greater than expected according 
to Molière’s theory (in its approximative 
form), namely Vi. 

In the present experiment the velocity of the pions f varied between 0.51 
and 0.66. The average value of 2, for this velocity range was Q = 630 per 
100 um, and by Fig. 1 the average value of B was B= 7.6, considering the 
fact that the cell size we have used in this work was 50 um. Since the velo- 
city of most of the pions was close to the average value of f and since B changes 
very slowly with the velocity, it is justified to use the above value for B for 
all the cases we had. The value of B, and B, according to (8) is found to be: 


Fig. 1. — The parameters 5, B, 
and B, as function of log, 2). 


(9) | B= 0 (we he, 


This can also be seen from Fig. 1. Let us use, for convenience, the nota- 
tion: (see (4) above): 


= nv BiB, 


= Xs / Li = x.VBi2. 


(5) S. Biswas, E.G. GEORGE, B. Perers and M.S8.Swami: Suppl. Nuovo Cimento, 12, 
369 (1954); M. Dr Corato, D. HirscaBERG and B. LocatELLI: Suppl. Nuovo Cimento, 
4, 448 (1956); S. Brswas, B. PETERS and Rama: Proc. Ind. Acad. Sci., 4, 154 (1955). 


a 
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By (4) and (10) we obtain (for ¢ = 50 um): 


| & = 0.356/p8, 
(10') 
| 03 = 0.448/pp 
3'2. The grouping of the data. — As was mentioned in Sect. 2, we have 


measured a total of 3960 angles. The mean absolute angle and the mean 
of the cosine function were obtained, as we shall see in Sect. 3°3 and 3°5, by 
treating the entire data as one group. For obtaining the dispersions, D/Vn 
however, it was necessary to divide the data in a suitable manner. The di- 
vision we have chosen is the following: 

99 groups, each containing 40 successive angles (namely 2000 um long 
segments of pion track length for each group) for the dispersion in the mean 
absolute angle case, and 396 groups each containing 10 successive angles for 
the cosine case. The reasons for this choice were as follows: 


a) In the mean absolute angle case one knows the relative error of the 
dispersion only for a normal distribution, in which case it is given by: 1/V2N 
(where N is the number of groups). It is reasonable to assume that for n = 40 
the distribution is close to a normal one, and one can use the above value for 
the error of the dispersion. On the other hand, in the cosine case one can 
calculate exactly the errors of the dispersion (see Sect. 3°6) and the requi- 
rement of a normal distribution is not necessary. This enables one to choose 
groups of ten angles (n= 10) in which case N becomes 396 and the statis- 
tical error of the dispersion becomes smaller. 


b) In principle, due to the correlation between the angles, D could depend 
upon n. In the cosine case, for example, it has been shown in I that: 


(11) {D[cos (np)]]?_ Di, Di 
Ü cosnp Sn" n° 


For all values of 7 that we have used, the second term on the right of (11) 
could be neglected for n => 10. Thus, in the cosine case, one can write the 
error as D/Vn ,where D is independent of » for n’s as small as 10. No such 
relation can be derived in the mean absolute angle case. One cannot calculate 
the dispersion in this case, and one does not know to what extent D depends 
upon ». Choosing n= 40 for the mean absolute angle case, makes the dis- 
persion less accurate than in the cosine case, but enables one to write the error 
also in this case as D/Vn, where D is probably independent of n. 
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33. The mean absolute angle. — In order to compare the experimental 


results with the theory, it is convenient to define an angle whose distribution 
is independent of the particles’ energy. Such an angle, for example, is: 


(12) “= g/o, 

where ¢ is the scattering angle and © is defined by (10). Using (10’) we obtain 
La = P2PB]0.356 , 

T3 = Pspp/0.448 . 

The mean of (12’) has to be compared with the theoretical value: 


(13) ass D = V/2/x (1 + 0.982/B, 4 — 0.117/Bi3+...), 


Lo,3 


which can be derived from (2), (3) and (4). 

This comparison is not straightforward, however, since, as was mentioned 
in Sect. 2, (12’) contains the scattering noise. The noise determination was 
carried out in the following way: the calibration track (see Sect. 2) was rotated 
slightly so that it formed a small angle with the X-direction of the scattering 
microscope. By this procedure a condition was reached which was similar 
to that under which the actual measurements of the pions were performed, 
The noise increased by about 20% in this case. The value finally used in the 
noise elimination was the one given in Sect. 2 increased by 18%, since it gave 
best agreement between theory and experiment for the case <|x,|) with cut- 
off. It should be emphasized that <|#|> is quite insensitive to the exact value 
of the noise. A change of 5% in the noise corresponds to a change of merely 
0.3% in <[x|}. 

Assuming a normal noise distribution, one gets (5) for the ratio of the noise 
in the second and third differences: 


In the present experiment we have obtained for the above ratio: 1.76 + 0.25. 
Assuming both the noise and the true scattering to have normal distributions, 
one gets the following relation between the mean absolute values of the ob- 
served angle 4, the true angle x, and the noise x,. 


(14) <a.) ={<|eg}=<e, |}. 


a 
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The final results of the mean absolute value <|x,|) are given in Table I, 
with and without cut-off at 4 <|a,|). The agreement with cut-off is much 
better than without cut-off, indicating too large a contribution to the result 
from large angle scattering in the theory of Molière. 


TaBLe I. — The mean absolute angle. 
fr EAPE NEADE | 
Difference - —— — —_ - —————_ 
Observed | Calculated Observed Calculated | 
pw } pri i Marr PME tri YE | 
No cut-off 0.855 + 0.013 | 0.895 0.864 + 0.017 0.891 
| | 
With cut-off | : 0.820 + 0.011 0.823 0.833 + 0.013 0.822 | 
| à | 


In order to reduce the numerical work, pf was assumed to be constant in 
each 500 um interval. Actually it varied by about 1%. However, the error 
introduced by this assumption was much less than the statistical error and 
could therefore be neglected. 


34. The dispersion of pp. Absolute value case. — In order to know how 
good an estimate of pf is obtained from eq. (2), one has to calculate the dis- 
persion of pf in the theory of Molière. However, such a calculation is impos- 
sible (*) since the second moment of n; diverges in Molière’s theory (*). One can, 
therefore, only determine the experimental value of the dispersion, whereas 
in the cosine case it is also possible to calculate it, as was shown in I. 

Using (14) we obtain for the dispersion of X,: 


as PIED_KaD Ara der {les flesh) 


[æ | «| ao |p? — <| ay |>? 


We have assumed that the noise has a normal distribution, therefore 
D(|a,y|)/|%y | = 0.755/Vn . 


If n is a large number and if the ratio <|#,|>/<|%|> is small, the distri- 
bution of <|X,|> becomes approximately normal and the dispersion of | X,| 


(*) Note added in proof: It should be noted, however, that it is possible to cal- 
culate the dispersion if an appropriate cut-off is applied. For example, HUYBRECHTS 
(Suppl. Nuovo Cimento, 4, 903 (1956)) calculated numerically the dispersion in the mean 
absolute angle with cut-off at 4¢y,> for B = 8, and got D, — 0.93, in agreement 
with the experimental value obtained in this work (see Table II). 
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given by (14) becomes the ordinary standard error. It can be written as 
D./Vn, where D, is the desired coefficient to be determined by experiment 
(as was mentioned above, D, diverges in Moliére’s theory). It should be noted 
that the dispersion of pf is the same as that of |X,|. In Table II are given 
the results of D, determined experimentally, of the second and third diffe- 
rences, for the cases of no cut-off (NCO), overall cut-off (CO) and local cut-off 
(LCO). The overall cut-off (a cut-off by using the average of all 3960 angles) 
is less useful than the local cuto-off (determined separately for each group 
of 40 angles), since the latter is of more practical importance. So, for ordinary 
emulsion work, the value of the dispersion for the LCO case should be used 
(see Table II). 


DARLE SEZ 


dispersion of pB in the mean absolute angle case. 


D,=Vn ED(p6)/P8)} Dig=vn[D(pB)/pBlh | 


CRAN IEEE ee ch 
| 


No cut-off MONTE 1.020 + 0.074 1.435 + 0.106 
With cut-off... 0.878 + 0.065 1.216 + 0.090 
Local cut-off 0.973 + 0.073 | 13722105102 
3'5. The mean of the cosine function. — The comparison between theory 


and experiment in the case of the cosine function has been performed for three 
different values of the multiplier 7 for the second as well as for the third dif- 
ference. The values of 7 used were 7/7, = #, 1, 2, where 7,,, is the optimal 
value of 7 which minimizes the dispersion of pf (see I). a 

In I the same parameter B has been used for both the tangent and second 
difference cases. A simpler and more general expression for the mean value 
of the cosine function can be derived by using the various parameters B, and 
X,; defined in (4) and (8). The result is: 


1 Il 
(16) ‘cos pp) = exp | 5 (70)? (1 sn In + (10?) 5 


where 


This expression is valid for the tangent case as well as for the second, 
third, etc., difference cases (see Fig. 5 of I for the graphical representation 
of (16)). The minimum dispersion of pf in the case of the second difference 
occurs at the point 


(17) (Moot, 02)? = 1.04, 


=) 
in) 
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which yields for the present experiment 

(alg) Mont, = 2-87 pB , 

and in the case of the third difference occurs at the point 


(18) (Mont, 03)” = 1.32, 
which yields 
(18’) Mont, = 2-57 DB , 


where pf is expressed in MeV/c. 

The noise elimination in the case of the cosine function is somewhat dif- 
ferent from the noise elimination in the absolute angle case. In the latter we 
have to assume a normal distribution for the true scattering whereas in the 
cosine case no such assumption is necessary. The relation between the means 
of the cosine functions of the three scattering angles 9,, g, and y, is simply 
given by 


(19) (COS NP;> = <COS YPo)[KCOS NPy> . 


(19) is the analogous expression to (14). We shall assume again a normal 
noise distribution. Therefore we have: 


(20) «cos "Py? = exp [Es 3(70x)?] ’ 


The expression to be compared to the theory is: 


Peo 1 mes 
(21) = © <cos p> = + exp [(ox)*/2] (cos po» 


where N is the number of groups. In the cosine case we had N = 396, each 
group contained n=10 angles. The variation of the exponential function is 
very small compared to the standard deviation of the cosine function. It is 
therefore permitted to take it out of the sum in (21). 

Table III shows the results of the mean value of the cosine function for 
three different values of the multiplier 7. The agreement with theory is very 


good at all points. 
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TaBLE III. — The mean value of the cosine function. 
{cos 7p4) a (COS NP,>s 
Mas | Le dame 
Observed Calculated | Observed | Caleulated 

1 0.957 + 0.002’ 0.955 | 0.948 + 0.003 | 0.945 
1 i 0.572 + 0.009 0.570 0.490 + 0.013 | 0.500 

| | 
2 | 0.137 + 0.015 0.152 | 0.094 + 0.018 0.097 


3'6. The dispersion of the cosine function. — The dispersion of < cos 74) 
is given by 


Zn ie VE > (m;—<m»)*/(m) = DV? , 


where 


| 


Mi = COS 7:55; 


LE) 


= 


is the mean of the cosine function in a group containing n angles, and N is 
the number of groups. 

The theoretical value derived in I cannot be compared directly with the 
measured one, as it does not include the noise. Therefore we have calculated 
the dispersion of the expression 


cos NPo = {exp [in(g.+9y)] + exp [— in(9:+9y)]}. 


The standard deviation 0} of cos 7g is found to be 


(22) of =0; — }(cos 2ny,>(1 — <cos 2npy>) + <cos np.>*(1 — <cos npy)?) . 


6, is the standard deviation of cos 7w,. As the influence of the noise on o, 
is small, it has been decided to neglect the statistical correlation in the deri- 
vation of (22). However, in the numerical calculation of of, according to (22), 
the exact value of of including correlation, as derived in I, eq. (4.9) (for the 
second difference) has been used. 

. Taking into account (19), the dispersion of cos ng, becomes 


MEI ZI 
COS 7Po {cos 7pD (cos NYy> ’ 
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and the dispersion of <cos 79,) becomes 


D(cos ng ) (Pes el | 


208 : 2|3 
(23) J D(cos pu Ì 
COS NP: COS Po 


COS 7Dx o (= 


Because of the assumption of normal noise distribution we have 


D(cos n@y) Linea 
Se E ip sr 


The dispersion of pf is equal to (23) divided by the derivative 
din <cos ng.) /d(In pf) . 


In Table IV the measured and calculated dispersion of pf of the second 
and third difference for three different values of 7 are presented. 


TaBLe IV. — The dispersion of pp obtained by the cosine function. 
| D.,=Vn[D(pB)/pPh Dig=vn[D(pB)/pB}z 
nn, = em ee IIC ea SS 
ra Observed | Calculated Observed | Caleulated 
i 1.349 + 0.120 | 1.614 | 1.488 + 0.144 1.756 
1 0.976 + 0.039 | 1.022 1.362 + 0.061 1.229 
2 | 2.013 +0.222 | 1.741 3.061 + 0.612 3.038 


The theoretical dispersion of pf for the second difference is given in I, 
eq. (4.10). For the third difference the calculation is similar but more tedious. 
We shall give here only the final results, namely 


Deb SA, Q + = exp [22](exp [e, — 1]) + 2 sinh? (9/6) + 
Roe S| 4 
DE exp [2/3](exp [e — 1]) + 2 sinh? (92/3) + 


exp [22/3](exp [es — 1p} /vanco + 38,/2), 


bo] ke 


where 
Q = — 3(on)?[1 — (1/B) In (Go?r?)], 
e, = 270? In 4/B3, &=(yo3)?/3Bz; and ¢& — 8702 In 4/9B,. 


A graphical representation of Vn:(D;(pB))/pB as a function of (703)? is given 


in Fig. 2. 
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The agreement with theory is quite good, but less than in the mean cosine 
case. The statistical errors have been found with the help of known statistical 
methods (?). Explicitly the relative standard 
error of s is given by the expression ap i 


o(s) J 1 A j 
OnE + gq Wale’ — a] /x | | ; | 


where o? is the standard deviation of the 
cosine function and w, its fourth (central) 
moment. The factor (1/n)(u,/o' — 3) measures 
the deviation of the cosine function from 
normality. (Note that for a normal distri- 
bution the error of the dispersion is 1/V2N). 


For relatively small values of » this factor is * Vicasio Sta fa ci 
quite big and therefore cannot be neglected. 

Let us remark here that in the case of the a (703 Ÿ 

absolute angle it is impossible, according 0 dI TA 15 20 
to the theory of Molière, to estimate this Fig. 2. — The dispersion of pf in 
factor. This is the main reason why in the the third difference case. 


absolute angle case the number of angles n 
in each group has been chosen much bigger than in the cosine case (40 in- 
stead of 10; see Sect. 82). 


4. — The minimum possible dispersion of pf. 


In the previous section we have obtained the dispersion of pf by using 
the mean absolute angle and the cosine function of the scattering angles g;. 
In order to check this point we have made use of the following statistical 
theorem (8): Let g;, ...,9, be a set of n independent random variables and 
let A be a parameter which is to be estimated by the g,’s. Let f(y; A) be 
the distribution function of y; and A*(q,,...,~,) be the estimation function 
of A. Then if A* fulfils the following conditions: (1) A* is independent of A, 
(2) A* is a continuous function and has continuous derivatives 0A*/09, at 
all points @= (91; ...: gn); (except perhaps at a finite number of points), (3) 
the expectation value of A* is A: E[A*(q,...,Qn)]=A and (4) if the de- 


(7) H. Cramer: Mathematical Methods of Statistics (Princeton, 1946). 
(REC) 477. 
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rivative of/CA is bound, then, the following inequality holds: 


(24) wpa 471 > fu [fips A) in fear ap} | 


Since we are interested in estimating pf we shall choose A to be: 


Aes Be. 

Not all possible estimates of A fulfil the conditions (1)-(4) above. For 
example, the mean absolute angle with cut-off (see Sect. 3°3) does not fulfil 
conditions (1) and (4), and the cosine function does not fulfil condition (3). 
(In the cosine case this difficulty can easily be removed). Thus the conclusions 
which will be drawn by using the above theorem will not apply to all pos- 
sible estimates of p6 only to a certain class of estimates. 

We shall restrict ourselves to the calculation of the minimum dispersion 
in the case of no correlation between the g,’s and compare the results to the 
calculated dispersion for the cosine function. The actual experimental dis- 
persion, as we have seen in I and in Sect. 3 of this work is about 14% greater. 
(About half of the 14%, increase is due to the fact that the angles y, measured 
by the Fowler Co-ordinate Method are correlated, and the other half is con- 
tributed by the scattering noise). In the case where there is correlation the 
dispersion is increased somewhat, but the same holds for the minimum dis- 
persion. Therefore, the above restriction will not affect very much the re- 
lative position of the dispersions and comparison of the calculated minimum dis- 
persion with the dispersion of the idealized cosine case will indeed show whether 
or not a better estimate of pf exists. For the distribution function we shall 
choose Moliére’s function, neglecting the correction function f®(g) and higher 
order terms, namely: 


(25)  f(p; A) = (1/V24)[(2/Vx):exp [— y?/2A?] + (1/B)fA(p/V2A)], 
fa) = D'(x)-[ (x? —4)(In 4y +20 (x) +1] — 2aD(a) , 


where: «=qg/\/2:A, (x) is the error function and I(x) is defined by the 
integral: 
1 
dZ 
Ia). = ra (exp [1 — Z?]x? — 1). 
0 

The minimum dispersion, using (24) and (25), is obtained by a straight- 

forward but tedious calculation, which involves also numerical integration (°). 


. (9) 8. Rospenporrr: Ph.D. thesis (unpublished), The Hebrew University, Jeru- 
salem (1957). 
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The results for various values of the parameter B are given in Table V. By 
inspecting the table it becomes evident that the dispersion of pB obtained by 
the cosine method is very close to the minimum possible dispersion. Since, 
as we have seen in Sect. 3, the dispersion of pf in the cosine case is very close 


Tape V. — The minimum dispersion of pB and the dispersion in the cosine case with 
no correlation, as a function of 1/B. 
1/B | 0 O05 010025 1020 
Minimum possible dispersion 0.707 0.789 0.845 0.872 0.960 
Dispersion optained by USING | 
the cosine function . . . | 0.707 0.790 | 0.850 | 0.882 | 0.976 
| | | 


to that obtained from the mean absolute angle parameter, it can be concluded 
that among the class of functions considered here, no better estimate of pf 
than the existing ones can be found. 


5. — Conclusion. 


On the basis of the present work it can be concluded that the cosine function 
as an estimate of pf rests on a simple theoretical foundation and its theory 
agrees very well with experiment. In this theory it is possible to calculate 
all necessary moments and therefore to determine exactly the errors of the 
second, third and even higher differences. This is not possible in the case of 
the mean absolute angle because of the divergence of the second and bigher 
moments. On the other hand, it is less convenient to use the cosine function 
than the mean absolute angle for two reasons: 

(a) The optimum value of the multiplier 7 depends upon the value of pf. 
It is therefore necessary to begin with a trial value of 7 and to iterate several 
times. Because, however, of the insensitivity of the dispersion with variation 
in 7, it should be possible to find a good value of 7 by one or two iterations. 
An alternative way of finding the optimal 7 is to first find an estimate of pf 
by the mean absolute angle method. 

(b) The determination of the experimental mean of the cosine function 
involves much numerical work (our analysis was performed on the WEIZAC (*)) 
whereas the experimental mean of the absolute angle is found almost im- 
mediately. 


(*) Weizmann Institute Automatic Computer. 
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We have seen that in the mean absolute angle case one gets good agreement 
with theory by introducing an appropriate cut-off and that the dispersion 
in that case is close to the one obtained using the cosine parameter. This 
shows that both the absolute angle and the cosine function are good estimating 
functions for pf from the statistical point of view. 


We are grateful to Dr. G. YEKUTIELI and Dr. H. J. LIPKIN for many 
clarifying discussions and to Dr. J. GILLIS and Mr. S. BRUDNO for programming 
and actually performing the calculation on the WEIzAc. We are also thankful 
to C. CIVIDALLI and R. PARDO for performing the scanning and part of the 
numerical calculations. 


RIASSUNTO (*) 


Confrontando la teoria con l’esperienza abbiamo esaminato l’uso pratico della nuova 
funzione di scattering multiplo cos np introdotta da Lipkin (+) et al. Contemporanea- 
mente abbiamo confrontato con l’esperienza le differenze seconde e terze dell’angolo 
assoluto medio risultante dalla teoria dello scattering di Molière e determinato speri- 
mentalmente le rispettive dispersioni. (Nella teoria di Molière la dispersione teorica 
diverge). Un buon accordo tra teoria ed esperienza si trova per il coseno, sia per il 
valor medio e le dispersioni sia per l’angolo assoluto medio con cut-off. Esprimendo 
l’errore relativo di pf con D,3/y (per le differenze seconde e terze rispettivamente) 
abbiamo ottenuto per la funzione coseno: D,=0.976 e D3=1.36, e per l’angolo asso- 
luto medio (con cut-off a 4¢\y|>): D,=0.97 e D3=1.37. Si dimostra che, limitandosi 
a una determinata classe di valutazioni di pf, i suddetti valori di dispersione sono 
vicini alla dispersione minima possibile di pf. Si suggerisce anche una spiegazione del 
noto fatto sperimentale che il rapporto delle « costanti di scattering» /3/K, è di 
qualche percento maggiore del valore atteso 4/3. 


(*) Traduzione a cura della Redazione, 


re] 
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Summary. — The contributions made by the u and r°-mesons to the 
frequency of burst production under a thick shield of materials have been 
calculated by making use of the recent results of the soft cascade shower 
theory. The theoretical results thus obtained have been compared with 
the observational data given by various authors. Satisfactory agree- 
ments between the theoretical and observational results are found every- 
where, at high altitude, at sea-level and also underground, using Poisson’s 
distribution for the fluctuation, and taking N(£ = 2me?, t), as the 
average number of particles in a shower. Further it is seen that the 
burst calculations made with the above modifications lead us to decide 
unambiguously in favour of the + spin theory of the u-meson, in con- 
tradistinction to the uncertainties met with previously while trying to 
determine the spin of the u-meson from the observed burst rates. The 
effects of taking N(£=0, t) for the average number of particles in a 
shower or of using a Furry distribution for the fluctuation, have also been 
estimated. 


1. — Introduction. 


Large ionization bursts produced in an ion-chamber under a thick shield of 
materials have been observed by different investigators. The size frequency 
of such bursts has been calculated by CHRISTY and KUSAKA (1) and CHAKRA- 
BARTY (2). Their calculations are based on the assumption that a u-meson 
while traversing the thick layer of materials over the ionization chamber ra- 


(1) R. F. Curisty and S. Kusaka: Phys. Rev., 59, 414 (1941). 
(2) S. K. CHAKRABARTY: Ind. Journ. Phys., 16, 377 (1942). 
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diates a high energy photon which in turn produces the ionization burst through 
cascade multiplication in the material. The knock-on process contributes 
only a very small portion of the large bursts which they considered, and as 
such was not taken into consideration. Apart from the fact that the mass 
of the u-meson taken was smaller and the absolute intensity of the u-meson 
higher than their values now accepted as probable, there are also other sources 
of uncertainty in their results. This arises mainly from lack of exact knowledge 
about N(E, t), the average number of particles in a shower and the cascade 
distribution function. CHRISTY and KUSAKA have based their calculations on 
the Furry distribution formula and Serber and Snyder’s values for N(£, t) 
which except for E= 0, is not very dependable as has been discussed pre- 
viously (*). On the other hand CHAKRABARTY has used values for N(E,t) as 
obtained by BHABHA and CHAKRABARTY which again are rather too low for 
iron. In other respects we have, however, followed CHAKRABARTY in assuming 
that the probability of N particles appearing in a cascade shower instead of 
the average number N is given by Poisson’s distribution. Further we have 
taken the average of the total number of electrons that can be detected expe- 
rimentally in a shower at depth # to be equal to N(E=2mc?, t) rather than 
N(E=0, t) or N;(0, t). A comparison with the results obtained recently by 
some authors (4) while studying the development of cascade showers in mul- 
tiplate cloud chambers justifies this assumption (5). With this modification 
the calculation of the size frequency of burst production gives consistent re- 
sults irrespective of the shielding material over the ionization chamber and 
definitely supports the 4 spin theory of the u-meson which is now accepted 
as correct because of the evidences from the decay spectrum (5). 

At higher altitudes, the frequency of burst production is much larger 
than at sea-level, while the intensity of the u-mesonic radiation does not 
correspondingly increase. This indicates that there must be some additional 
mechanism which will contribute a major part to the frequency of burst pro- 
duction at high altitudes. It is believed that bursts recorded at high altitudes 
under thick layers of material are primarily due to the nuclear interaction of 
the N-rays. The r°-mesons produced within the shield by such interactions 
will, on account of their very short life-time, decay almost instantaneously 
giving two photons which in their turn will produce the ionization bursts 
through cascade multiplication. FAHY (°) has given a rough estimate of the 


K. CHAKRABARTY and M. R. Gupra: Phys. Rev., 101, 813 (1956). 
E. HAZEN: Phys. Rev., 99, 911 (1955); P. A. BENDER: Nuovo Cimento, 2, 


. R. Guerra: Proc. Nat. Inst. Sci. (Ind.) 23 A, 217 (1957). 
SAGANE, W. L. GARDNER and H. W. HuBBARD: Phys. Rev., 82, 557 (1951). 
. F. FAHY: Phys. Rev., 83, 413 (1951). 
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bursts size-frequency based on the neutral meson hypothesis. In the present 
paper, we have made an attempt to determine the contribution to the burst 
frequency from r°-mesons at high altitudes and also at sea-level. 


2. — Bursts produced by u-mesons. 


It has been shown previously by CHAKRABARTY that for large bursts con- 
taining more than 100 particles the contribution from the knock-on electrons 
produced by the u-mesons is negligible. The u-meson may, however, radiate 
a high energy photon in the field of the atomic nucleus. The cross-section for 
this process (assuming spin 4 and magnetic moment e4/2Mec) is given by 
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when the screening of the nucleus by the orbital electrons is incomplete and 
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If N(E,t)is the average number of particles of energy > at a deptht 
in à cascade shower initiated by a primary of energy Æ,, then assuming Pois- 
sonian distribution, the probability of finding more than N particles when 
the average is N is given by 
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Thus the probability that a u-meson after traversing a thickness 7 of the 
shielding material above the ionization chamber will register a burst consisting 
of more than N particles (electrons or positrons) is 
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where / (in em) is the thickness of the material of the shield corresponding to 
one radiation unit and o is the number of atoms per cm? in the materials. 
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The differential spectrum of the u-mesons is given by (5) 


x i (2-10° eV} 
(4) {0 TA (eT 4 i do , 


where 6 = 1.87, eT=1.8-10° eV at sea level and the absolute intensity J, as 
given by PUPPI and DALLAPORTA (°) is- 8.3-10-*s—' (em?)~ sr!  Herce the 
frequency per cm? per s of bursts containing more than N particles is 
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To evaluate the integral on the right hand side we require the values of V 
for different values of Æ, and t. This has been given in closed analytic form 
by CHAKRABARTY and GUPTA (*). Since the exact expressions for N are too 
complicated to be used here, we have, following CHAKRABARTY (?) tried to 
set up empirical relations which give values in close agreement with the nu- 
merical values they have given. We have taken 
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where %, %, %, ete., will depend on the cut-off energy £ and also on the 
material in which the shower develops. 
For E=2me?, we have (for lead) 


(9a) a= OTs ag 205 ey PR = = 0 = 
while for iron we have for the same energy value, 

(90) daro = 20 5 cd IS RARI 
For E=0, we have both for Pb and Fe 
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(8) H. EuLer and W. HeISENBERG: Ærgeb. Exakt. Naturwiss., 17. 1 (1938). 
(*) G. Puprr and N. DALLAPORTA: Progr. Cosm. Rays Phys. Vol. 1, p. 336. 
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For sufficiently thick shields, i.e., for large values of 7, we may now write 
the integral over t in (5) as 
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For N > N,,, the value of the integral (10) is insignificant as pointed out by 
CHAKRABARTY. Since we are interested in large bursts containing more than 
200 particles, this means that an appreciable contribution to B,(N) will come 
only from energies E, > E,, where H,= E, corresponds to N,,= 200. For 
the different cases that we are considering here, y,= In (£,/}) will be bet- 
ween 7.5 and 8.0. Further as the meson energy w > E, > H,, we may ap- 
proximate the meson spectrum S(w) by taking w~°*” instead of (eT+@)7®*P. 

Interchanging the © and €, integration, substituting for N,, and ¢,,, (5) may 
now be written (for incomplete screening) as follows 
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For reasons mentioned in the previous section, we may put (yo — £:)/(%o — ))}R 1 
and take the lower limit e= Be” in the integral over H,. (11) can now be 
easily evaluated and we have finally for the case of incomplete screening, 
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and similarly for complete screening 
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and y is the Euler constant. 
The function f,(N, 6) and f,(N, 6) are in the same form as given by CHAKRA- 
BARTY except that values of the constants 7 and ¢€ are different. 
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B,(N) as given by (13) cannot be directly used for the calculation of the 
frequency of small bursts (50 < N < 100) produced underground. The reason 
for this is that at large depths (say at a depth of 60 m w.e.) the differential 
u-spectrum is given by 
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so that for the burst size under consideration we are no longer justified in 
replacing (12-109 + w)~°*” simply by © +. We have therefore taken 


(2-10° eV)? 


S(c) = "079-61, di , 
with dò = 1.55, which gives a fairly good fit with the actual value for $(©) 
for energies down to — 6 GeV. The contribution to the underground small 
bursts through the bremsstrahlung process can now be evaluated by mul- 
tiplying the right hand side of (13) by the factor .079 and putting dò = 1.55. 

For small bursts the contribution from the knock-on electrons produced 
by the u-mesons by direct collision is also important. The cross-section for 
the knock-on process is 
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where €, is the maximum energy transferable to the electron in a direct col- 
lision, and is given by 
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Substituting for Q(«w, E) in (5) and proceeding as before, the contribution of 
the knock-on electrons can now be determined. The results have been shown 
in Fig. 3 and are discussed in a subsequent section. 


8. — Bursts produced in 7°-mesons. 


For the calculation of the size frequency of bursts produced by 7°-mesons 
it is necessary to have a knowledge of its differential energy spectrum at dif- 
ferent heights above sea-level. But sufficient direct observational data 
regarding this are not available. Following the theory of nucleonic cascades 
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based on the hypothesis of plural production of mesons (1°), we may however 
obtain it from the spectrum of the primary cosmic ray component. The dif- 
ferential spectrum of the primary cosmic-ray component (11) may be approx- 
imated by a power law of the form E7°*”. The value of 6 is however not 
constant but varies from 1’1 (in the low energy region £ < 10 GeV) to the value 
1.67 for sufficiently high energies. In the energy range in which we are inte- 
rested we may put ò= 1.67 and write the primary cosmic ray component in 
the form 


> (10° eV)? AE 
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where a= 1.67. 

It has been shown by MESSEL that for H > Æ, the latitude cut-off energy, 
the primary component is absorbed exponentially in the atmosphere so that 
the spectrum at a depth x g/em? below the top of the atmosphere is 


S(E, x)dE = S(E)exp[— x/X]dE, 


where X g/em? is the absorption mean free path in the atmosphere. 
Assuming the top of the atmosphere to be horizontal the intensity per 
cm? per s at depth x is 
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At a depth ¢ (in radiation units) within the shield over the ionization chamber, 
the intensity of the incident nucleonic component is 


27f(x/X)S(E)dE exp [—t/T], 


where 7 is the absorption mean free path in radiation units for the material 
of the shield. Now, the probability per collision length of a nucleon of 
energy £ colliding with another nucleon and giving rise to a fast nucleon in 
the energy range e, e+de, a meson in the range ©, w+dm and a second 


(19) W. HEITLER and L. JANossy: Proc. Phys. Soc., A 62, 374 (1949); H. MessEL: 
Proc. Phys. Soc., A 64, 726 (1951). 
(1) H. V. NEWER: Progr. Cosm. Ray Phys., Vol. 1, p. 301. 
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nucleon of energy Æ — e—@ is (1) 
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Then the total number of mesons produced within a thickness d(7,— t) (Tp is 
the total thickness of the shield) within the shield is 
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where / and Z are respectively the radiation and collision unit for the ma- 
terial of the shield. Of these mesons, one-third are neutral 7°-mesons which 
decay almost instantaneously into photons. Assuming that the decay photons 
are emitted isotropically in à co-ordinate system moving with the velocity 
of the neutral meson, it may be easily shown that the probability that a decay 
photon will have energy in the range E,, £,4+d4, is (1) 
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and is zero outside this interval, where W is the energy of the meson and © 
is its rest energy. For © <a, we may put 
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(2) H. Messer, R. B. Ports and C. B. A. McCusker: Phil. Mag., 43, 889 (1592). 
(18) A. G. CARLSON, J. E. Hooper and D. T. Kine: Phil. Mag., 41, 701 (1950). 
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Combining (2), (19) and (20) and remembering that two photons are obtained 
from a single decay event, we find that the frequency of bursts produced by 
the neutral mesons is given by 
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For large thicknesses of the shielding material, the integral over { may be 
written as 
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Substituting for N and proceeding as in (10) the expression (23) may be shown 
to be equal to 
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Substituting (24) for the integral over # in (20) we finally get 
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4. — Numerical results. 


The integrals f,(N, 6) and f.(N, 6) have been evaluated numerically by the 


saddle point method and then from 
equations (13a) numerical values for 
B,(N) have been obtained. This has 
been done both for iron and lead shields 
over the ionization chamber and the 
two sets of results thus obtained have 
been compared with the sea-level data 
given by LAPP (1). It should be men- 
tioned here that while Lapp’s burst 
data for iron are obtained with a 35 cm 
homogeneous iron-shield, it is not so for 
lead (10.7 cm thick shield) because of the 
presence of the 1.25 cm thick iron-wall 
of the ionization chamber. The effect of 
the iron which is to reduce the number 
of cascade particles produced in the 
lead shield, may be taken into consider- 
ation by choosing a value of # inter- 
mediate between those for pure lead and 
iron. 
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Fig. 2. — Dashed curve shows Lapp’s sea- 
level data under 10.7 em Pb+1.25 cm Fe 
shield. All the other curves indicate the 
calculated results (Eqn. (13a)) under the 
same shield. O and K- Christy and Kusaka’s 
curve; Oh-> Chakrabarty’s curve. 
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Following the arguments given by 
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Fig. 1. {- Plot of log, B,(N) against 

log,, N. Dashed curve shows Lapp’s 

sea-level data. F indicates Furry distri- 

bution. The screening is taken to be 
incomplete. 


CHRISTY and KUSAKA we find that 
this effective value f,,.,, is 12 MeV 
in this case. For the iron-shield, 
there is, however no such difficulty 
and we have f,, equal to 25.9 MeV 
as usual. 

In Fig. 1 and Fig. 2 we have 
plotted log, B,(N) against log, N 
for different assumptions regarding 
the spin of the u-meson and the ave- 
rage number of particles. Compar- 
ison with Lappîs sea-level data 
(Shown in the same figures) show 
that the 4 spin theory of the u-meson 
is the most plausible one, provided 
we take N(E = 2mc?,t) as the ave- 
rage number of particles. If how- 
ever, we use N(E= 0, t) the theore- 


(14) R. E. Lapp: Phys. Rev., 69, 321 (1946). 
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tical burst frequencies (for spin 3) become larger than the observed ones. In 
this case the spin zero curve gives a better agreement. The reason for this 
discrepancy is that the cascade shower theory in its usual B-approximation 
form gives too large a multiplication near zero energy (as the assumptions 
regarding constant ionization loss etc., are no longer valid.) Further, below a 
certain energy the particles will escape detection and finally it is to be remem- 
bered that the energy of a particle (electron) at minimum ionization is nearer 
to 2mc? rather than to zero energy. 

In Fig. 2 we find that, as expected, the values given by CHAKRABARTY are 
somewhat less than those given by us for H=2mc?. On the other hand, 
contrary to expectation, CHRISTY and KUSAKA’s values (after introducing the 
necessary corrections for the mass of the u-meson ete.) are less than ours 
for H=0. The reason for this is that at a later stage in their analysis they 
have introduced a modification which is equivalent to taking H,/15f as the 
maximum number of particles, while it is nearly H,/96 according to SERBER 
and SNYDER (which corresponds to our values for N,,,, when E= 0). More- 
over the empirical expression for N(H,t) which they have used gives t,,,.=7 
independently of the primary energy while actuallly it varies nearly linearly 
with In (Z,/B). According to their analysis the contribution to the burst fre- 
quencies for any particle of energy 4, is proportional to t., and hence the 
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resulting frequency is thus further under-estimated. | 
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Fig. 3. — A — George and Trent’s under- Fig. 4. — A, B— from Carmichael 
ground bursts data under 15 cm Pb shield; and Steltje’s sea-level burst data under 
B and ©- contributions from brems- 12 em and 27 em Pb respectively; 
strahlung, knock-on electrons respectively ; C, D, E- from yu and r°-mesons at 
Pts. @ > combined contribution from the sea level under 12 em and 27 em of 


two processes. Pb respectively. 
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In Fig. 3 and 4 we have compared our results with the experimental data 
given by LAPP, FAHY and GEROGE and TRENT (!). It is seen that according 
to our analysis nearly 75% of the bursts at sea-level (under 10.7 cm Pb + 
+ 1.25 cm Fe) should be ascribed to u-mesons, while the 7°mesons contri- 
bute (16 to 20)% of the total under the same shield. At an altitude of 3500 m 
above the sea-level, the u-mesons contribute about 10% of the total while 
(70 to 75)% comes from the 7°-mesons. As regards the small bursts observed 
under ground (at a depth 60 m w.e.) itis found that (90 to 95)% of the total 
burst frequency may be accounted for by considering the combined effect of 
the knock-on electrons and the bremsstrahlung photons emitted by the u-mesons. 
As expected the contribution from the r°-mesons at this depth is negligible. 

In Fig. 4 we have also compared our results with those given by CAR- 
MICHAEL and STELJES (!) at sea-level under lead shields of different thicknesses. 
The contribution from the u-mesons amounts to nearly 79% of the bursts ob- 
served under a 12 cm thick Pb shield, while under a thickness of 27 cm of Pb 
nearly 85% of the total may be accounted for in this way. For these two cases 
the contribution from the 7°mesons has also been worked out and is found 
to be & 21% and 12% respectively. Our analysis thus gives results which 
are fairly consistent with the experimental ones at different heights and depths 
and also under different thicknesses of material. 


5. — Dependence on the form of fluctuations. 


So long we have only used the Poissonian distribution in our calculations. 
To have an idea of the effect of fluctuation on the frequency of burst pro- 
duction we shall now calculate the burst frequency by assuming the Furry 
distribution. In this case 


For convenience in numerical calculation we have used a different empi- 
rical form for N(E, t). Following CHRISTY and KUSAKA we have taken 


where #,, is the same as given by (8) and (9) and 


[a=15 for lead (F—2mc!), 


(28) | = 11.5 for iron (E=2me?), 
| = 9.0 for both Pb and Fe (E=0). 


(5) E. P. GeorGe and P. T. Trent: Proc. Phys. Soc., A 64, 733 (1951). 
(15) H. CARMICHAEL and J. F. STELJES: Phys. Rev., 99, 1542 (1955); 105, 1626 (1957). 
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Substituting for P(N+1, N) and N(E, t) from the above equations in (5) 
it may be easily shown that 
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and A and L are the same as given by equations (14 a, b). 

In Fig. 1 we have plotted log, B,(N) as given by (29) for the two dif- 
ferent cases viz., E=2me? and E= 0. It is quite evident from the graphs that 
whatever be the form for fluctuation, N(E= 0, +) cannot be taken as the 
average number. However, for E = 2mc?, the experimental curve lies in bet- 
ween the two extreme cases given by the Poisson and the Furry distribution 
and the former is to be used for this value of the cut-off energy. But the 
Furry distribution may also be successfully used if a higher value of the cut-off 
energy is chosen. 


I am indebted to Dr. S. K. CHAKRABARTY, F.N.I. for help and guidance 
in this work. 

The author is also thankful to the Atomic Energy Department, Govt. 
of India for the financial support during the progress of this work. 


RIASSUNTO (*) 


Si è calcolato il contributo dato dai mesoni u e 7° alla frequenza della produzione 
di burst sotto uno spesso schermo di materia servendosi di recenti risultati della teoria 
della cascata molle. I risultati teorici così ottenuti sono stati confrontati con i dati 
d’osservazione forniti da vari autori. Si è trovato dappertutto un soddisfacente accordo 
tra i risultati delle osservazioni e la teoria, sia a grande altezza, sia al livello del mare, 
sia sotto terra, adottando per la fluttuazione della distribuzione di Poisson e prendendo 
N(E=2me?,t) per il numero medio di particelle in uno sciame. Si vede inoltre che i 
calcoli sui burst fatti con le modificazioni di cui sopra ci portano senza ambiguità a 
preferire la teoria dello spin + del mesone pu a differenza della incertezze incontrate prima 
nel cercare di determinare lo spin del mesone u dalla frequenza dei burst osservati. 
Si sono anche valutati gli effetti del prendere N(£=0, #) per il numero medio di par- 
ticelle in uno sciame o dell’usare una distribuzione di Furry per la fluttuazione. 


(*) Traduzione a cura della Redazione. 
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The Average Multiplicity and Inelasticity 
in r-Meson Production in the Atmosphere. 


G. ISHIKAWA and K. MAEDA (*) 


Meteorological Research Institute - Tokyo 


(ricevuto il 29 Agosto 1957) 


Summary. — The energy spectra of primary and secondary cosmic ra- 
diation are analysed to find the multiplicity-energy relation in the z-meson 
production in the atmosphere at the energy region from several GeV to 
some thousands GeV with the help of simplified considerations on the 
meson production by nucleus-nucleus collisions. The results being sub- 
jected to the assumed values of inelasticity suggest that the multiplicity- 
energy relation agrees with the formula of the type given by Fermi and 
Landau under the condition of a practically constant small inelasticity. 


1. — Introduction. 


Using the data of the integral energy spectrum of the N-component averaged 
over the atmospheric depth J(Æ,), and that of charged x-mesons //(£,, E,) 
the empirical multiplicity of 7-meson production averaged throughout the 
atmosphere m(H,, H,) is given by OLBERT in the simplified expression 


OII | od 
E ’ o) — Am l An ? 
Ronin 

where 7, is the so-called cut-off energy which is well defined by the geomagnetic 


latitude 2 (1). 
For the evaluation of //(£,, £,) the number of charged 7-mesons produced 


(*) Now at University of Nebraska, Lincoln, Nebraska. 
(1) S. OLBERT: Phys. Rev., 96, 1400 (1954). 


54 G. ISHIKAWA and K. MAEDA 


downwards throughout the atmosphere with kinetic energy greater than E 
at the geomagnetic latitude 2, OLBERT made use of his corrected spectrum 
of u-meson production under the condition of H, = 260 MeV. And J(E,) 
is taken from the data compiled by G. PuPPI and N. DALLAPORTA (2). | 

Since his results are derived from the latitude effect of cosmic-rays, the 
obtained multiplicity covers only the low energy region of primary radiation 
i. e. less than about 15 GeV. It is obvious, however, that the estimation of mul- 
tiplicities of meson production in the atmosphere can be extended to higher 
energies as long as both the energy spectra of primary particles. and of the 
produced 7-mesons are known. 

Therefore, we have tried to estimate the average multiplicities of 
meson production in the higher energy region up to 10?GeV, making 
simplified assumptions for the meson production by nuclear collisions to get 
the energy relations between the spectra of primary and of secondary 
COSMIC-TAYS. 


2. — The energy relation of primary particles and produced 7x-mesons. 


We consider the general casea nucleus-nucleus collision in which one 
of the nuclei of mass A,M, has a total energy A,U, in the laboratory system 
and the other of mass 4,M is originally at rest, where M is the mass of the 
nucleon and U, is the total energy of the incident particles per nucleon. As 
we consider such à high energy collision that the binding energy is negligible, 
the total energy carried by both nuclei before collision in the center of mass 
system, U*, is written as 


(2.1) U* ={ Mc?[2A4,4,U, + (Aj + A;)Me]} 


and the velocity of the center of mass of the two nuclei, /., is 


A; 

(2.2) Be = (U2 — Met) [(u+ = Met] | 
\ Sei 

Now let m(U,; A,, 4,) and 7(U,; A,, A.) be the multiplicity and inelasti- 
city respectively for the meson production in the above nuclear interaction, 
then the total energy of a meson produced in the center of mass system, UÈ, 
is by definition 

7 QUE A; A,) 


2.3 pt ieee) U*. 
us m(U,; A1; A.) 


(2) Progress in Cosmic Ray Physics, 1, 317 (1952). 
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It can be assumed that one half of the emitted mesons have the average 
angle of emission, 0*, measured from the line o fflight of the incident particle, 
and the rest are emitted symmetrically backwards, and that the transferable 
energy for meson production is equally distributed to each x-meson in the 
center of mass. 

Hence, using the above equations, the energy of a x-meson, U_, with 
respect to the laboratory system is given as a function of the incident energy 
U, by 

n (Uy; Ay, Ay) - 


(2.4) Ur — mt A 3.) LA UA Mo? tA, cos 0*V U2— Meet] 0<0*<x/2, 
DI HORS 


where -- indicates the emission of the mesons respectively forward and 
backward in respect to the direction of the incident particle (see Appendix). 


3. — The energy spectra of primary particles and of produced 7x-mesons. 


The production spectrum of m-mesons is directly related to the energy 
spectrum of the parent N-rays. The actual relationship of these two quantities 
might be very complicated, but this complexity can be reduced by considering 
the statistical quantities of averaged multiplicity m as follows. 

Let f(U_) be the effective energy spectrum of charged z-mesons at pro- 
duction and P(U,) be the energy spectrum of N-rays, then the relation of 
these two quantities is given by 


(3.1) KU.) = 3m(U,)-P(O,) . 


The factor + in (3.1) arises from the fact that on the average only # 
among the produced mesons are charged and among them only 3 is directed 
forward. As shown in Appendix II, the contribution of the primary particles 
which produces 7-mesons backwards in the center of mass system is in this 
case negligible. 

N-rays include protons, neutrons, z-mesons, «-particles and heavy nuclei 
contribute meson production colliding with the air nuclei. The amount of 
heavy particles in cosmic-rays is very small so that we can neglect its con- 
tribution to meson production. On the other hand, the contribution of «-par- 
ticles to meson production is not to be ignored especially in the topmost 
layers of the earth’s atmosphere where the abundance of «-particles is con- 
siderable. We assume here, for simplicity, that the secondary nucleons and 
r-mesons do not play an important role in the meson production. Then, we 
may write 


(3.2) m(U,) = e,m(U,; 1, Ae) + e,m(U,; 4, Ae) &p+e,= 1 
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where e, and e, are the ratio of protons and «-particles to the total number 
of N-components of cosmic-rays in the atmosphere, respectively. The maximum 
value of e, may not exceed 10 percent (*), then the expression for m(U,) may 
be written in first approximation, : 


(3.3) m(U,) = A-m(U,; 1, As) , 


where 


e, and e, depend not only on the atmospheric depth owing to their different 
attenuation mean free path in the air, but also on the magnetic latitude, 
because the magnetic cut-off energy of these particles is not the same. In 
a similar fashion, the average inelasticity for meson production in a nucleus- 
nucleus collision can be approximated by 


(3.4) n(U,) = A, (Up; 1, A,). 


Hence, the average total energy of the energy of a z-meson produced by the 
primary particle of energy U,, can be written from (2.4) 


an [A,U, + A,Me? + À, cos OV U2 — Me] , 


(3.5) C= 
where U_ is U* in the expression of (2.4). 

The equation (3.5) shows that the average multiplicity in a nucleus-nucleus 
collision can be obtained from the knowledge of the energy spectrum of primary 
particles and the production spectrum of z-mesons, provided that the unknown 
quantity is determined in some other way. 


4. — The production spectrum of z-mesons and w-mesons, and their energy 
relation. 


The effective production spectrum of z-mesons f(U_) in the equation (3.1) 
is not obtained directly from observation. Therefore, instead of the direct 
use of f(U_), we introduce a more convenient quantity, the production spectrum 
of u-mesons. 

The relation between the energy spectrum of z-mesons and that of u-mesons 


(3) Progress in Cosmic Ray Physics, 1 232 (1952). 


56 


a 


THE AVERAGE MULTIPLICITY AND INELASTICITY IN Æ-MESON PRODUCTION ETC. DI 
at production is given with sufficient accuracy in the high energy region by 


B._() 
mu 


T 


(4.1) us(U,, 0) = TU), 


where U,=U_ and r is the ratio of the mass of the u-meson to that of the 
m-meson and is taken as 0.76 (‘). (Because of the two body-decay of the 
m-meson, the energy relation of U, and U, is given definitely, and the above 
expression is obtained with extreme relativistic approximation). 

On the other hand, the energy spectrum of z-mesons is given as a solution 
of the diffusion equation, 


O(U,,#) _ f(Ux) 


(4.2) : NEB orp Pe 


L 


È cui x) , 


SIETE J 
Ox A PROS gir, 


where A and L are the collision and absorption mean thicknesses of the primary 
particles respectively, and A_ is the nuclear absorption mean free path of 
m-mesons. In what follows, we shall assume that À is equal to L, having 
the constant value of 120 g/em?, and that 2. is 60 g/cm? in the air. cf_ is 
the velocity of z-mesons in the laboratory system and B_(x) is the characte- 
ristic energy for the decay of x-mesons at the atmospheric depth x, which is 
given by 


(4.3) Br) = 2, 


where M_ and 7_ are the mass and the mean life of z-mesons and o(x) is the 
air density at the atmospheric depth x g/em?. 

The solution of (4.2) was numerically calculated disregarding the mul- 
tiplicity in meson production and assuming B_(r) to be a constant (*). In the 
actual atmosphere x/o(x) is not necessarily constant and should be carefully 
examined in particular when one considers the temperature effect at low 
altitude. Therefore, we have improved this point as shown in Appendix III 
and obtained the solution in the form 


(4.4) UA) ESE 
where 

3 “i 1 MATE A È it z x az f = 
(4.5) (UR. x) == 110 exp | A == 3) = DO Or LS XY, 


0 


(4) K. MAEDA and M. Wapa: Journ. Sci. Res. Inst., 48, 71 (1954). 
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and 

(4.5) Se ALe 
ol È sala (= 4) az + 
+e Eo ga) ]4 


Then, the relation (4.1) can be written 


À Re 
(4.6) BDs 0) = HU) pg SU 2) 

In Fig. 1 and Fig. 2, the values of x(U,,x) and u;(rU_,%) under the 
above assumption are shown for the atmospheric depth x, with three different 
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values of the energy U_. In these calculations, we made use of the empirical 
formula of the energy spectrum of primary particles given by BARRET et al. (5) 
for P(U), which is expressed by the following formula 


Aer = (0819-1074 Î 


(4.7) 1 


7 


3 V 2.25 + 00.2 ln ((U +3.2)/3.2) 
sim) 


U+3.2 
1.35+0.04 In | o) x 


a 


where U is measured in GeV. 
(5) P. H. Barret, L. M. BoLLINGER, G. Cocconi, Y. EISENBERG and K. GREISEN: 
Rev. Mod. Phys., 24, 133 (1952). 
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From Fig. 2, one can see that the production spectrum of the u-mesons is 
well expressed by the following form in the region of the atmospheric depth 
larger than 10 g/em?, 


(4.3) UTU,, ©) = p(rU_)-exp (— 2/L) . 


The attenuation length L, is found to be nearly equal to that of the primary 
particles, ~ 120 g/cm?. 

The fact, that the depth where the production of u-mesons reaches its maximum 
value is shifted up to around 10 g/cm? whereas that of the 7-mesons is in the 
vicinity of the effective absorption depth of primary particles, can be under- 
stood by considering that the production of y-mesons is proportional not 
only to the intensity, but also to the decay rate of the 7-mesons, and that the 
latter increases exponentially with the altitude. 

Hence, we get the following expression for the production spectrum of 
u-mesons with the help of equations (3.1), (4.6) and (4.8), 


m(U,)-P(U, 
(4.9) ps(t Un) = ps (1U,) SnD . 


where uî(rU_) is the vertical production spectrum of u-mesons calculated with 
the assmption of f(U) — P(U). 


5. — Calculation of the average multiplicity. 


The above considerations show that the multiplicity of 7-meson production 
can be obtained as a function of the energy U,, using the equations (2.4), (4.7), 
(4.9) and the calculated values of w*(rU_) provided that u;(rU_) is known. 

For this purpose the empirical production spectrum of u-mesons given by 
OLBERT expressed in range spectrum (') has been used for numerical calcu- 
lation of u,(rU_) with the help of an energy-range relation. 

In Fig. 3, uy(rU_) and uf(rU_) are shown together with the primary energy 
spectrum given by (4.7). 

Letting 

Lis 


(5.1) De) 41) rx 


AU) 
rU,) | 


The average multiplicity of z-mesons, m(U,) becomes 
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The relation between U, and U, is given by the equation (2.4). The 
actual procedure is rather tedious, namely, first, in order to calculate É(U_) 
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Fig. 3a. 
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for a certain value of U,, we ‘adopt arbitrary values of U, being consi- 
dered eventually as the energy of N-rays which produce forward directed 
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m-mesons of energy U_. 

With these values of U,, we 
get m(U,)’s by the equation (5.2). 
Next, in order to check which of 
these m(U,)s thus obtained is 
correct, we put them back into 
the equation (2.4). If one of the 
values of U_ coincides with the 
one formerly assigned, the corre- 
sponding U,, and consequently 
m(U,), is adopted as the real one. 

The second procedure has been 
carried out actually for several 


different values of 7(U,) under the assumption that 7 is independent of U,. 
The numerical values we used are: Mc? = 0.938 GeV, M_c® = 0.141 GeV. 


and 


A,=1.15, 


À, = 14.4 


cos 0È = 0.786 . 
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The results obtained are shown in Fig. 4. For a comparison, the theoretical 
values computed by Fermis formula of the multiplicity-energy relation is 
shown by a dashed line which is expressed by 


(5.3) . M(U,) = 1.8(0,/Mer)* , 


provided that the ratio of the number of charged mesons to the total is two 
thirds. 


6. — Coneluding remarks. 


The multiplicity calculated by our model depends on the assumed in- 
elasticity. As one can see from Fig. 4, our evaluated multiplicity increases 
with decreasing primary energy below around 20 GeV if 7 is larger than 0.2. 
This apparent increase of multiplicity at lower energies for large values 
oi inelasticity seems absurd, because in the actual case the emission of 7-mesons 
should be suppressed by the threshold energy of multiple meson production. 
As it is plausible to assume that the multiplicity is a monotone increasing 
function of the primary energy, our results indicate either that the inelasticity 
decreases with the U, decrease in this low energy region, or that it has a 
practically constant value less than 0.3. 

It is interesting to notice that the recent analysis by VERNOV and GrI- 
GOROV shows that the inelasticity in the energy region below 20 GeV is about 
0.2 ~ 0.3, increasing with the primary energy (°). 

With respect to the multiplicity in meson production, several theoretical 
estimations have been given especially for the analysis of high energy events 
such as jet phenomena in nuclear emulsions and air showers (71). 

However, recent experiments shows that the energy dependence of mul- 
tiplicity in meson production is favourable to the results given by Fermi’s 
statistical theory and Landau’s relativistic-hydrodynamic theory as long as 
the energy dependence of the multiplicity is concerned. 


S. N. Vernov and N. L. Gricorov: Suppl. Nuovo Cimento, 4, 879 (1956). 
E. FERMI: Prog. Theor. Phys., 5, 570 (1950). 

W. HrIsenBERG: Zetts. f. Phys., 126, 569 (1949). 

H. W. Lewis, J. R. OPPENHEIMER and S. A. WOUTHUYSEN: Phys. Rev., 73, 


948). 
(10) S. Z. BeLEN'KJI and L. D. Lanpau: Suppl. Nuovo Cimento, 3, 15 (1956). 


61 


62 i G. ISHIKAWA and K. MAEDA 


If we express our results in an empirical formula of the type 


(6.1) m(U,) = a(n Uy)”, 


where x and y will be determined for a certain value of 7 by means of a least 
_ squares calculation, we get for the case of 7= 0.2, as an example, x = 1.665 
y= 0.249, and the equation (6.1) is written in this case 


(6.2) m(U,) = 1.665-(0.2-U,)}. 


This is practically consistent with Fermi’s formula (5.3). Although, in 
Fermi’s theory, the inelasticity is not explicitly described, it is obvious that 
the inelasticity is small because among the secondary particles are con- 
sidered not only z-mesons, but nucleons and anti-nucleons. If one of the se- 
condary nucleons is emitted with large energy, this corresponds to the low 
inelasticity because this secondary nucleon is undistinguishable with the in- 
cident primary nucleon. 

On the other hand, Landau’s formula for the multiplicity in a nucleon- 
nucleus collisions is a little different from that of Fermi’s nucleon-nucleon 
collision and is expressed by (1) 


Uae 
(6.3) m(U,) = KA?” (sae | 


Comparing our results with this formula, we get 


We would like to thank to Dr. MARCEL SCHEIN for his kind advice and 
comment on this subject at the Midwest Cosmic Ray Conference at Iowa 
City. We are also very grateful to Prof. ROBERT CHASSON in the University 
of Nebraska and to Dr. JUN NISHIMURA of Tokyo University for their helpful 
discussions. 
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APPENDIX 


I. — The relations (2.1), (2.2) and (2.4) are derived as follows. 
The Lorentz invariance with respect to both systems of center of mass 
(C.M. system) and laboratory (L. system) is 


(I-1) c?p° an U2 — cpt — U*2 e 


The C.M. system is defined so that p* is zero, and if cB be the velocity of the 
incident particle with respect to the L. system, then its momentum is given by 


(I-2) cp = AU, 6 
and this can be also written 

(1-3) B=1—(Me@/U,) 
and 

(1-4) U— AU, + A,Mc. 


Putting above relation into (1-1), the equation (2.1) is obtained. 
Since p* is zero with respect to the C.M. system, the Lorentz transformation 
of the total energy U* to the L. system is given by 


(1-5) U == U* |W 1 = Bee 
where cf. is the velocity of the center of mass in the L. system. This is written 


also, 
(1-6) bid (CFU). 


Putting (2.1) and (I-4) into (I-6), the equation (2.2) is obtained. 

Similarly, the Lorentz transformation of total energy of each x-meson, 
U*, emitted at an angle 0% with the direction of the incident particle in 
the C.M. system, to the total energy in the L. system U_, is written as 


(1-7) U, = (US + Beep* cos 6*)/V1 — B2, 


where cp* = U** and f, is given by the equation (2.2). Using (1-5), this 
can be written further 


U + 
(1-8) U;— U* (Ba dE Pb COS 05] e 


Since we are considering the extremely high energetic case of U,, the relative 
velocity of the emitted x-mesons with respect to the C.M. system fz, can be 


Gol 
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regarded as unity. Then this equation becomes (2.4), where the corresponding 
terms of 6, and U*/U* are substituted by (2.2) and (2.3), respectively. 


II. — The actual production of z-mesons, fa(Ux), can be expressed by two 
terms, i.e. the one due to the mesons emitted forwards and the other due to 
those emitted backwards in the C.M. system of colliding particles. 

Let f'(U_) and f(U_) be the production spectra corresponding to these 
two terms respectively, then we can write 


fa(U,) = f'(U,) + f(U,). 


As the emission of z-mesons in the C.M. system is regarded to. be sym- 
metrical with respect to the direction of the incident particle, we get 


and 


f(U,) = 3m(U,)P(U,) O( Uz) O(U;) <1, 


where U;= is the energy of a primary particle which produces z-mesons of 
energy U, in the L. system originally emitted forward and backward in the 
C.M. system respectively 0(U;) is a correction factor for angular spread in 
the case of originally backward emitted mesons. 

Then from the equation (2.4), we get the following ratio for the extreme 
relativistic case, 


Us As co8:05 


Us _1— cos 0% m(U>) m(Uz) 
UR 00807 


(Un) mU 
n(U}) MU) : 


The inelasticity might be almost independent of energy and the multiplicity 
can be regarded as a slowly increasing function of energy, i.e. 


DI + 
n E, a a (as it is obvious U; < U;). 
D D 


Therefore, we can write 


Uy _1— cos Or 
Ue acon Oa 


According to the Fermi’s statistical theory of multiple meson production, the 
angular distribution of 7-mesons in the C.M. system is a function of the impact 
parameter of colliding nucleons and is shown approximately be const. cos 0%: 
For this case, we get 0%— cos! 0.79 as the average angle of emission. Put 
0*= 0È in the above inequality, then we get 


HO 21 


SEE Va 
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This ratio is an overestimation, because the recent analysis of jet showers in 
nuclear emulsions shows that the seconadry particles are more concentrated 
within the narrow cone around the axis than shown by the above expression, 
consenting with Landau’s theory which gives very sharp forward angular 
distribution in C.M. system. Consequently, the mean angle of emission 0* should 
be far smaller than the above estimation. It follows, therefore, that the ratio 
of U; to U; might be less than of order 0.1. 

Now, the differential energy spectrum of primary ray is approximately 
shown by 


Pa for U, > 20 GeV. 


Then, the ratio of the primary flux which contributes to the backward and 
forward production of z-mesons of energy U,, is 


On the other hand, the ratio m(U;)/m(U;) is not so small compared with 
the above quantity. 

Therefore, the contribution of the parent particles corresponding to the 
backward produced z-mesons can be neglected as compared with that of pro- 
ducing forward primaries. 


III. — The atmospheric depth x (in g/em?) is equal to the barometric pres- 
sure P (in dyne/cm?) divided by the acceleration of earth’s gravity, g, i.e. 


VEDE CI] 


and the pressure P is also given by P — oRT, where 9 =o(x) and T= T(x) 

are the density of air and the temperature at the atmospheric depth x, res- 

pectively. Thus we get 

(III-1) eee 
Therefore, the characteristic energy for the decay of r-mesons (4.3) is 


written by 


(III-2) Bla) = a T (a) 


= Ta), 


where H = 0.521 GeV deg“! in our case. 
In order to get the rigorous solution of the diffusion equation (4.2), we 
made the following assumption of the vertical temperature distribution in 
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the atmosphere, 


Ta)=T, for #< a, 
(III-3) 
=T,+a(a— x) for TT 


Pi 21705 ree 190de gros Oey emis, 


Putting the conditions (III-3) into the equation (4.2), we get the two 
different equations for x < x,, and consequently the two solutions (4.5) 
and (4.5’). 


RIASSUNTO (*) 


Si analizzano gli spettri energetici della radiazione cosmica primaria e secondaria 
per trovare la relazione molteplicita-energia nella produzione di mesoni x nell’atmo- 
sfera nell’intervallo energetico da alcuni GeV ad alcune migliaia di GeV con l’ausilio 
di considerazioni semplificate sulla produzione di mesoni per collisione nucleo-nucleo. 
Dai risultati, confrontati coi valori assunti per l’inelasticità, si rileva che la relazione 
molteplicità-energia si accorda con una formula del tipo di quella data da Fermi e 
Landau a condizione di avere debole inelasticità, praticamente costante. 


(*) Traduzione a cura della Redazione. 
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Shrinkage of Nuclear Research Emulsions. 


W. M. GiBson and J. G. McEwEN 


Physics Department, Queen’s University 
Belfast, North Ireland 


(ricevuto il 31 Agosto 1957) 


Summary. — A method is described for obtaining accurate values for 
the shrinkage factor of nuclear research emulsions, and for testing the 
uniformity of the shrinkage with depth. The results obtained from one 
plate show a shrinkage factor near the surface 6% greater than near the 
glass. The importance of avoiding such variations, when accurate depth 
measurements are required, is stressed. 


Nuclear research emulsions contain a large proportion of silver bromide, 
removal of which in processing reduces their thickness by a factor of 2 or 
more. If accurate measurements in three dimensions are to be made on events 
recorded in the emulsion, this factor must be known. Several methods of 
measuring it have been used, the most accurate being that described by 
DUKE, Lock, MARCH and Munir (1): they used a mechanical gauge to measure 
the difference between the unprocessed and processed thicknesses, while the 
thickness of the processed emulsion was measured with a calibrated microscope. 

The present work was undertaken partly to check the accuracy of measu- 
rements made by the method of Duke et al., and partly to study the uniformity 
of shrinkage throughout the depth of the emulsion. 

A beam of 925 MeV protons from the Birmingham synchrotron was ar- 


(1) P. J. Dure, W. 0. Lock, P. V. Marcu and B. A. Munir: Journ. Sci. Instr., 
32, 365 (1955). 
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ranged to pass through emulsions, nominally 400 um thick, held on a frame 
which could be rotated about an axis at right angles to the beam. The plates 
were set at approximately 45° from the beam direction for one pulse from the 
synchrotron, and then turned through exactly 90° for a second pulse. The 
average angle of dip of the two beams was thus 45°, within about 0.1°, even 
if the plates were exposed in their original wrappings and with 1° uncertainty 
in their orientation. 

The results from three sets of plates are described here: plates A were 
kept in the manufacturers wrapping until a short time before exposure; 
B were unwrapped and kept in an atmosphere at (58 + 5)% relative humi- 
dity for 10 days before exposure; C were kept in this atmosphere for 7 days, 
then sealed with a sheet of thin lacquered copper, and left for a further 7 days 
before exposure. 

The plates were processed in the usual manner (?), and were examined on 
a microscope whose fine-focus scale had been shown to be reliable by an inter- 
ference method (8). 

The examination was carried out in a room with controlled relative humi- 
dity (58 + 5)%, after the plates had been kept there for some days to reach 
equilibrium. The vertical setting was noted at horizontal intervals of 50 um 
along about 50 tracks of each of the two beams, starting 10 um from the end 
nearest to the glass. 

Since the beam was not perfectly parallel, and the emulsions were originally 
about 420 um thick, some tracks gave 7 intervals and some 8: in the latter 
case the middle two intervals were grouped together. The average of the 
angles of dip of the two beams in the processed emulsion was obtained for each 
interval, and the shrinkage factor calculated. Other values were obtained 
by measurements over the whole length of each track, and by the gauge method 
of DUKE et al. The results are shown in Table I, for three sets of plates exposed, 
processed and examined at different times. The statistical uncertainty in the 
figures from the gauge method is 0.02, and in each of the other figures 0.01. 

It will be seen that there is very satisfactory agreement between the means 
of the measurements throughout the depth, the independent measurements 
over the whole length, and the gauge measurement where this was made. But | 
there is strong evidence of a systematic variation of shrinkage factor throughout 
the depth of plate A. In plate B there is a low value for the layer next to 
the glass, but this might possibly result from a systematic error in the method 
of using the microscope, and would need to be confirmed. In plate C, slightly 
lower values are found for the layers closer to the centre of the emulsion. 


(2). A. D. Darnton, A. R. GATTIRER and W. O. Lock: Phil. Mag., 42, 396 (1951). 
(3) J. B. ENGLAND: Journ. Sci. Instr., 34, 208 (1957). 
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The variations observed could in principle be due to any of three causes: 


a) a non-uniform response of the gelatine to the changes involved in 
processing, or 


b) a change in atmospheric humidity some time before exposure, or 


c) a change in atmospheric humidity some time before observation. 


TABLE: I: 
| Mean depth Shrinkage factor 
Part of tracks used of interval ri 
(microns above glass) Plate A | Plate B | Plate C 
LE Ca sE Si e 
First 50 um (near glass) . | 35 2.38 | 2.25 2.34 
SECON desi 85 OS 4-00 TR es 
Th SOU 2 lc oes | 135 944 - |} 9.29 2.30 
Middle 50 um or 100 u AA 200 2.48 2.31 2.30 
Second from last 50 um . 265 2.46 232 DBS 
Next to last 50um .. . 315 | ol 2.82 2.36 
Last 50 um (near surface) 365 | Des 2.35 
| 
NEbVDolellengthi.c. i. ito 208 | 2.46 il 2.34 
Emulsion measurements by gauge and microscope | — DEE Doo 


The procedure used prevented b) and ec) from playing any significant part: 
b) might have contributed slightly in plate A, but the whole effect would have 
required a very large and sudden change in atmospheric humidity (4). But 
in this connection it should be pointed out that the characteristic time for 
the diffusion of water into or out of a 400 ym thick emulsion is 64 hours; this 
represents the time which must elapse after a change in surface conditions 
before the change in mean density of an emulsion reaches 93% of its final 
value, or before the change in density at the bottom reaches 89% of its final 
value (see the analogous theory for conduction of heat, developed for example 
by CARSLAW and JAEGER (°)). It was obtained by measurements in conditions 
where water vapour was removed rapidly from the surface, and a longer 
effective time constant would be observed in a stagnant atmosphere. Careful 
planning is therefore needed before any exposure in which uniformity of the 
shrinkage with depth is required. In some conditions it may be easier to 
sacrifice some of the advantages of thick emulsion for the sake of the smaller 
time constant of a thinner layer (xt? where t = thickness). 


(4) A. J. OLIVER: Rev. Sci. Instr., 25, 326 (1954). 
(5) H. S. Carstaw and J. C. JAEGER: Conduction of Heat in Solids (Oxford, 1947). 
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Once an emulsion bas reached equilibrium, however, the slow rate of dif- 
fusion of water vapour through the surrounding air thereafter helps to keep 
the distribution of water within the emulsion uniform. In experiments ana- 
logous to those described above, but with deliberate extreme alteration of the 
humidity 12-24 hours before exposure, no systematic non-uniformity was 
found. Little would be expected unless conditions allowed unusually rapid 
transport of water vapour. 

If a) is responsible for the observed effect, it would represent a type of 
distortion not previously observed, but by no means unlikely. It could perhaps 
be reduced by the same precautions as are used to minimize distortion in the 
horizontal plane. 

The large effect observed in plates A, and in some others not described 
here, might, however, be due to slight solubility of one of the components of 
the gelatine; this could lead to a reduction of the thickness of the upper layer 
of the processed emulsion, and a greater shrinkage factor for these layers. 
This could be prevented, or compensated, only by specially designed proce- 
dures. 

It may be concluded, therefore, that if depth measurements are required 
to be accurate to better than + 3%, not only the actual shrinkage factor 
but also the uniformity of shrinkage should be checked: care should be taken 
to ensure that no diffusion of moisture is taking place at the time of the expo- 
sure, and a plate from the same batch should be exposed in some apparatus 
similar to that described here and processed with the plates under test. If 
this shows that the shrinkage is not sufficiently uniform the processing pro- 
edure should be modified. 


The authors wish to thank Dr. W. O. Lock and Mr. P. J. DUKE for ar- 
ranging the exposures, and Mrs. S. E. BAKER and Mrs. D. E. ANDERSON for 
examining the plates. 


RUCASISTUENERO) 


Si descrive un metodo per ottenere accurati valori del fattore di contrazione delle 
emulsioni per ricerche nucleari e per saggiare l’uniformità della contrazione con la pro- 
fondità. I risultati ottenuti con una lastra rivelano un fattore di contrazione del 6% 
maggiore presso la superficie che non vicino al vetro. Si mette in rilievo l’importanza 
di evitare tali variazioni quando sono richieste precise misure in profondità. 


(*) Traduzione a cura della Redazione. 
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Summary. — Measurements of ionization loss by relativistic charged 
particles have beeu made in helium, argon and xenon. The experimental 
results agree well with theory in the case of helium. In argon, and par- 
ticularly in xenon, the relativistic rise in ionization found experimentally 
is less than that predicted theoretically. The measurements were made 
by the method of counting drops in a cloud chamber. 


1. — Introduction. 


The theory of energy loss by ionization was first developed by Bour (1) 
in 1913. In the early 1930’s BETHE (?) and BLOCH (*), improved the treat- 
ment by using quantum mechanical methods. The theory they developed 
indicated that the energy loss per unit path length would decrease as 1/f? 
where f=v/c, v is the velocity of the particle and c is the velocity of light, 
until b=.97. At this velocity the energy loss per unit path length would begin 
increasing as log By, where y= 1/(1— #2)’, and would increase indefinitely 
as the velocity increased. The increase is caused by the relativistic lateral 
extension of the electric field of the moving particle. The lateral extension 


(*) Supported in part by the joint program of the U.S. Office of Naval Research 
and the U.S. Atomic Energy Commission. 

(+) Now at E. I. duPont de Nemours and Company, Wilmington, Delaware. 

(*) On leave of absence from the Laboratoire de Physique, Ecole Polytechnique, 
Paris. 

(2) N. Bour: Phil. Mag., 25, 10 (1913); 30, 581 (1915). 

(2) H. A. Berne: Zeits. f. Phys., 76, 293 (1932); Handb. der Phys., 24, 515 (1933). 

(3) F. BLocH: Ann. der Phys., 16, 285 (1933); Zeits. f. Phys., 81, 363 (1933). 
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of the field allows the particle to excite or ionize a larger number of atoms 
as it passes through the medium. 

In 1938, SWANN (‘), suggested that the presence of the medium through 
which the particle was passing would limit the relativistic extension of the 
electric field because of the polarizability of the atoms in the medium. In 1940 
FERMI (5) calculated this polarization effect or density effect, as it is more 
commonly called, for the special case of a medium composed of dispersion 
oscillators of a single frequency. The conclusions of this theory indicated 
that the energy loss per unit path length should increase to some maximum 
value which depended on the density of the medium. 

Fermi’s theory has been extended by several authors (°°) to multifrequency 
models of a dispersive medium, but the general conclusions as given above 
have not been altered. The disagreement between the theories now depends 
on detailed differences in the calculations, differences between dispersion models 
assumed and differences in the interpretation of the role played by Cerenkov 
radiation. 

Accurate experimental results on the relativistic rise in ionization have 
been obtained by GHOSH, JONES and WILSON (1°), who used u-mesons in a 
cloud chamber, and a high resolution magnetic spectrometer. Their measure- 
ments extended over a momentum range from 500 MeV/c to about 30 GeV/e 
and showed a relativistic rise which agreed with theory and also showed the 
existence of the density effect, but did not reach the plateau. PRICE (11) has 
summarized other work on this problem. 

Recently an experiment was performed by BARBER (*) in which he used 
intense beams of electrons having energies from 2 to 35 MeV and measured 
their energy loss by ionization in an ionization chamber. He measured the 
energy loss in hydrogen and helium at two pressures, one and ten atmospheres. 
The experimental points on the plateau, which were observable in both gases 
at ten atmospheres, were a little below the value predicted by STERNHEIMER 
but were not outside the experimental error. 


4 F. C. Swann: Journ. Franklin Inst., 226, 598 (1938). 


(4) W. 
(©) E. Fermi: Phys. Rev., 56, 1242 (1939); 57, 485 (1940). 
(5) O. HALPERN and H. Harri: Phys. Rev., 57, 459 (1940); 78, 477 (1948). 
(7) G. C. Wick: Ric. Scient., 11, 273 (1940); 12, 858 (1941); Nuovo Cimento, 1, 
302 da 

(8) M. SCHÔNBERG: Bull. Cent. Phys. Nucl. Bruxelles, No. 20 (1950); Nuovo Ci- 
mento, 9, 210, 372 (1952). 

(°) R. M. STERNHEIMER: Phys. Rev., 88, 851 (1952); 89, 1148 (1953); 91, 256 (1953); 
93, 1434 (1954); 93, 351 (1954); 103, 511 (1956). 

(10) S. K. GHosH, G. M. D. B. Jones and J. G. WiLson: Proc. Phys. Soc., A 65, 
68 (1952); 67, 331 (1954). 

(11) B. T. PRICE: Report on Prog. in Phys., 18, 52 (1955). 


(2) W. C. BARBER: Phys. Rev., 163, 1281 (1956). 
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The possibility that the relativistic increase in ionization could be used 
in the measurement of the velocity of very high energy particles has led us 
to measure this rise in various gases. 

We used a cloud chamber placed in a magnetic field so that the momentum 
of the particles could be determined by measuring the curvature of the track. 
The ionization was determined by stereoscopically counting drops along the 
track. In order to obtain measurements of ionization over the range from 
By=3 (minimum ionization) to By—1000 both u-mesons and electrons were 
used. The measurements were carried out in helium at about 1.3 atm., in 
argon at about 0.2 atm., in argon plus helium each at about 0.2 atm., and 
in xenon plus helium each at about 0.1 atm. In each case there was in the 
cloud chamber, in addition to the inert gas, alcohol and water vapor at a par- 
tial pressure of about 5 cm Hg. 


2. — Experimental procedure. 


The magnet and cloud chamber used for the helium experiment have been 
described previously (%). In order to measure the ionization in argon and 
xenon it was necessary to construct a new cloud chamber to operate at pres- 
sures below atmospheric pressure. This chamber is slightly larger than the 
previous one but is very similar in construction and operating characteristics 
and fits into the same magnet. 


3. — Photography. 


During the helium experiment the camera was operated at f/16 and power 
for the flash tubes was supplied by about 300 uF charged to 2000 V. A delay 
of 140 ms between the expansion of the cloud chamber and the flashing of 
the lights was used to allow time for the drops to grow to a size large enough 
for photography. A 6 ms delay was added electronically to the unavoidable 
mechanical delay of 12 ms between detection of an event and complete expan- 
sion of the cloud chamber. This delay allowed the ions to diffuse sufficiently 
so that the drops which later formed on them could be counted readily. 

During the argon experiment the camera aperture was opened to f/11 and 
the power for the flash tubes was increased by doubling the number of con- 
densers connected to each lamp. In addition the delay was increased to 250 ms. 
These changes were made because the drops grew much more slowly in the 
argon atmosphere. The rate of drop growth is limited mainly by two things; 


(3) W. B. FrETTER and E. W. FRIESEN: Rev. Sci. Instr., 26, 703 (1955). 
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by diffusion, because as the drop grows it depletes the vapor in the immediate 
vicinity and can grow only as rapidly as vapor diffuses to it, and by thermal 
conductivity, because as the vapor condenses heat is liberated, and the faster 
this heat can get away the faster the drop will grow. 

In the argon experiment it was readily seen that diffusion was not causing 
the difficulty because the ions were diffusing about the same distance as they 
had in helium, because of the lower pressure of argon (0.2 atm.). Therefore, 
we assumed that the drop growth was being limited by thermal conductivity, 
and verified this by adding helium to the gas mixture which caused the drops 
to grow much more rapidly. The thermal conductivity of helium is about 
nine times as large as the thermal conductivity of argon. 

In the argon-helium experiment (0.2 atm each), when the cloud chamber 
was filled with equal pressures of helium and argon, the camera aperture was 
f/16, and the flash lamp delay was 140 ms. 

With the xenon and helium mixture (0.1 atm. each) the camera was set 
at f/16 and the light delay at 240 ms. It was desirable to use f/16 because 
the photographic depth of field was then sufficient to give good drop images 
to within about one centimeter of both the front glass and the back plate. 
With the aperture set at f/11 the depth of field is about 2/3 the total depth 
of the cloud chamber. Therefore, during the pure argon experiment each 
lens of the camera was focused on a slightly different plane in the chamber, 
and good drop images were obtained on either one view or the other over the 
full depth of the chamber. 


4. — Operating procedure. 


Since we desired to obtain the relative value of energy loss by ionization 
over the velocity range which covers from minimum ionization to the plateau 
value, we decided to use both u-mesons and electrons. The u-mesons cover 
the range from minimum ionization (6y=3) to fy=80. The electrons cover 
the range from fy=50 to By=1000. 

In order to compare the electron and u-meson data it was necessary to 
take both types of pictures under exactly the same experimental conditions. 
To accomplish this the pictures were taken alternately, two u-meson pictures 
and then one electron picture and this cycle was repeated automatically. 

In a cloud chamber the condensation efficiency depends on the expansion 
ratio, the temperature, and possibly on the time after last expansion. There- 
fore, great effort was made to Keep these factors unchanged from expansion 
to expansion during each experiment. The expansion ratio controller was 
not touched during each run and the magnet temperature, although automa- 
tically controlled, was checked frequently and found not to vary more than 
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0.2 °C, except at the beginning of the xenon experiment where it was 0.5° 
lower than the average temperature. A 23 minute dead time was used between 
pictures to smooth out the variation in the time between expansions due to 
fluctuations in the coincidence rate. 

The argon and helium used was commercial grade, about 99.7% pure. 
The xenon that was used was considerably purer. The vapor was obtained 
from a mixture of ethyl alcohol and water placed in a cellulose sponge 
pad in the bottom of the chamber. In the helium experiment a mixture of 
2:1 by volume ethyl alcohol to water was used, and in the argon and xenon 
experiments a mixture of 3:2 was used. 

In order to eliminate the air from the cloud chamber prior to filling it for 
the helium experiment, the chamber was flushed ten times to a pressure of 
1.5 atmospheres with helium. For the xenon and argon experiments the cham- 
ber was evacuated with a vacuum pump. 


5. — Data analysis. 


The momenta of the u-mesons were measured by methods previously 
described (78), simplified where possible. For the low energy electrons the only 
difference in the momentum determination procedure was that instead of 
using comparator measurements to determine the curvature, the photographs 
were projected to full size and the track curvature was compared to arcs of 
circles of known radius. 

We estimate that the above procedure for determining momentum was 
accurate to about 5% and since an error of 10% in the momentum corresponds 
to an error of only 1% in ionization, we considered this sufficiently accurate. 


6. — Ionization measurement. 


The ionization was measured by counting drops along the track with a 
stereoscopic microscope which magnified the images 40 times. The drops 
were counted between two easily recognized spots along the track (e.g., where 
the track crossed fiducial lines) and then later, during the projection of the 
event, the distance between these points was measured. 

The drops were counted by pressing a switch connected to a counter for 
every drop seen. When a large cluster of drops was seen, the number of drops 
in the cluster or blob was estimated, and if it exceeded 40, none of the drops 
were included in the data; but the fact that a blob had been seen was noted. 
It was assumed that all of the ion pairs in one of these clusters were secon- 
dary ions produced by a single collision between the primary particle and 
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one of the electrons in the gas. Therefore, the forty drops were a measure 
of the maximum energy transfer in a single collision that was considered in 
the experiment. This maximum energy transfer is calculated below. 

In order to be able to determine the variation of the consideration effi- 
ciency from the top to the bottom of the cloud chamber, the u-meson tracks 
were divided into two sections for drop counting. The top half and the bot- 
tom half of the track were counted separately. It had been shown previously 
that the variation in the condensation efficiency was approximately linear, 
at least nearly enough for the purposes of this experiment (!); therefore, two 
points were enough to determine the curve. We corrected for this variation 
in the case of the electrons which were not centered in the cloud chamber. 

Since the error in ionization was directly dependent on the error in the 
measurement of the length of the track counted, considerable care was taken 
to keep errors in this measurement at a minimum. For the mesons, which 
were full length tracks almost parallel to the back plate (the counter telescope 
limited the angle with respect to the plane of the back plate), the length was 
measured on the full scale reprojection of the event with a ruler. The radius 
of curvature of the lowest momentum pu-meson was large enough so that the 
track length differed from the chord length by less than one per cent. The 
length measured was then corrected for magnification, and a certain length 
was subtracted for each blob detected. 

Since the electron tracks were so curved that the length of the chord was 
not a sufficiently accurate measurement of the length of the arc, and since 
the radius of curvature changed enough along the track (due to the inhomo- 
geneous magnetic field) so that a correction could not be applied to the length 
of the chord, a calibrated wheel that is used in measuring distances on maps 
was rolled along the reprojected image of the track to measure their length. 
In addition to the corrections applied to the u-meson track lengths, the elec- 
tron track lengths had to be corrected for angle with respect to the back plate. 

The number of drops per centimeter was then determined by dividing 
the total number of drops counted by the correct length of track. This number 
is a measure of the average amount of energy lost per centimeter by ioniza- 
tion in energy transfers less than a certain maximum. 


7. — Systematic errors and corrections applied. 


In order to eliminate as many systematic errors as possible, the following 
procedure was used in analyzing the data: 


1) The particles were divided into three or four groups, depending on 
the experiment. Each group had an equal number of each kind of particle 
and an equal number of particles for each momentum range. 
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2) Each of the three or four groups was analyzed by a different indi- 
vidual. 


3) Electrons and u-mesons were intermixed during the drop counting, 
and the pictures were counted in a random order. 


It was expected that there might be systematic differences in the absolute 
value of the drop counts between individuals, but this turned out not to be 
the case. In each experiment the average values of the minimum ionization 
determined for each of the groups were within one standard error of the over- 
all average. Therefore, all of the grups were mixed together. 

In each experiment it was found that the minimum level changed gradually 
during the course of the experiment. In the helium and xenon experiments 
it slowly decreased and in the argon experiments it decreased. The reason 
for these changes is not definitely known, but all effects that have been thought 
of that could produce such a result (i.e., a leak or a slow systematic change 
in the condensation efficiency) would have negligible or no effect on the rela- 
tive ionization. It was assumed that the variations were a linear funetion 
of the event number and a least squares line was fitted to the variation. The 
original data were then corrected. on the basis of the least squares line. 

The time between expansions fluctuated from the 23 minutes dead time to 
over 60 minutes for some electron events. The data were examined to see if 
the different waiting times introduced any systematic difference in condensa- 
tion efficiency between the electron showers and the u-meson events. The 
minimum level determined from pu-mesons was plotted as a function of time 
after last expansion and no systematie variations were observed. 

Another possible source of systematic error is the drop overlap correction. 
As the density of ions produced increases in general, the number of drops 
counted increases but not quite as fast. Previous work done in this labora- 
tory shows that in the present experiments this effect is negligible. The tracks 
used for these experiments were 1.5 times as wide as those used in reference (18) 
which means that the density of drops was down by over a factor of two for 
the same number of drops per centimeter. Therefore, in the argon experi- 
ment where the drops were most dense (34 drops/em) the correction as obtained 
from the previous work is of the order of 0.1 per cent. 

The statistical fluctuations that were found to exist in the determination 
of the ionization of a 40 cm track at minimum ionization were calculated. 
In the helium experiment the standard deviation of a single measurement 
was found to be 8%, in the argon experiments it was 13%, and in the xenon- 
helium experiment it was 11%. The total number of drops in a 40 cm mini- 
mum ionizing track was, in each of the respective experiments, 788, 924 and 640. 
By comparing the standard deviation calculated from all electrons together 
to that calculated by taking the average standard deviation calculated from 
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the single pictures which contained several electrons, it was found that a large 
part of error was due to fluctuations in the minimum level from picture to 
picture. The standard deviation calculated from all electrons was about 
twice that obtained on the single picture basis for each of the experiments. 

When the u-meson data were compared with theory, care had to be taken 
in the high momentum region. In this region the momentum errors are large, 
and the number of mesons per momentum interval decreases rapidly as the 
momentum increases. Therefore, errors in measurement would throw low 
momentum mesons into high momentum groups and thus cause a systematic 
error. 

The procedure used in this experiment to estimate this error was to cal- 
culate the curve of ionization versus momentum that would be expected, 
taking into account the differential spectrum of the mesons, the statistical 
error in a given momentum measurement, and the theoretical curve of ioniza- 
tion versus momentum. This was done as follows. 

Let p, be the experimentally determined momentum for a particular u-me- 
son. In a cloud chamber momentum determination the error in the reciprocal 
of the momentum is a constant for all momenta and the error is gaussian. 
Therefore, in order to make the calculations easier, the inverse of the momen- 
tum is used. The relative probability that a particle actually has the momen- 
tum p when the experimentally determined momentum is p, can be written 


as 
iL 1 1}? 
wo [eft (242), 
pene oe | 
where 1/p' represents a systematic error in the momentum determination 


and k is a constant which determines the width of the guassian distribution 
in the error. 

The quantities k and 1/p’ are measured by periodically turning off the 
magnetic field and photographing a high energy u-meson. The curvature 
of the tracks obtained is measured in the comparator in the same way that 
the curvature of a normal track is measured; and since the particle went 
straight through the chamber, any curvature detected is due to turbu- 
lence in the chamber and measurement errors. In our cloud chamber the 
maximum detectable momentum which is defined as that momentum which 
gives a curvature equal to the probable spurious curvature arising out of 
the turbulence and measurement errors, was 33 GeV for the argon experi- 
ment, 40 GeV for the argon experiment and 30 GeV for the xenon experiment. 
There was no significant systematic error in the momentum measurements 
of the helium experiment, but in the argon experiment p’ was 55 GeV and in 
the xenon experiment p' was 21 GeV. The systematic error means that on 
the average a straight track would be measured to have the momentum p’. 


PT 
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In order to calculate the ionization that would be expected on the average 
for a u-meson with an experimentally determined momentum of p,, one has 
to sum over all possible momenta p, multiplying the probability that a meson 
with a measured momentum p, actually has momentum p times the relative 
number of that momentum and finally times the ionization produced by a 
particle with momentum p. The expression written in terms of the inverse 
momentum is: 


co 


SI (1/p)n(1/p) exp [— k {1/p — (1/p, + 1/p’)}?] 4(1/p) 
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where the denominator is a normalizing factor and n(p) = n(1/p)[ (d(1/»))/dp] 
is the differential momentum spectrum of the u-mesons incident on the cloud 
chamber. In the integral, negative momentum means that the charge of the 
particle with this momentum is opposite to the charge of a particle which 
has positive momentum of the same absolute value. 

This expression was integrated numerically for each experiment using the 
experimental n(p) as given by CARO et al (1). In general it was found that 
the theoretical and experimental curves began to differ at about 8 GeV/c; 
therefore, no u-mesons with momentum greater than this were included in 
the results. 


8. — Comparison of results with theory. 


The results of the present experiment have been compared with two theo- 
retical expressions for the energy loss by ionization (see Fig. 1-4), one derived 
by STERNHEIMER (1952) and one by BUDINI (?°). 


The formula obtained by STERNHEIMER for the total energy loss in col- 
lisions with energy transfer less than a certain amount 7 is: 
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where N is the number of electrons per cm’, m is the mass of an electron. 


(4) D. E. Caro, J. K. Parry and H. D. RATHGEBER: Nature, Lond., 165, 689 
(1950). 
(5) P. BupINI: Nuovo Cimento, 10, 236 (1953). 
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e is the electronic charge, and v is the velocity of the ionizing particle. I is 
an experimentally determined quantity called the average ionization potential 
and ò, which is the density effect correction, is given by STERNHEIMER in 
terms of the plasma frequency of the medium, the oscillator strengths and 
transition frequencies of the energy levels of the atoms in the gas. STERNHEI- 
MER adjusts the observed optical frequencies and ionization potentials to 
give the experimentally observed value of the average ionization potential 
for a given material. Since these values are uncertain we used two different 
sets to compare with our data. 

In order to compare the predictions of Sternheimer’s theory with our re- 
sults, a correction has to be made for the energy loss in the form of Cerenkov 
radiation which escapes to distances greater than the distance b from the track 
of the ionizing particle; b is chosen so that energy absorbed inside that distance 
would be counted experimentally. 

Sternheimer’s correction for energy loss in the form of Cerenkov radiation 
for the helium experiment lowers the plateau value less than one per cent 
for the mixture used in the cloud chamber; whereas for pure helium the pla- 
teau value would be lowered by 34 or 4 per cent. This happens because there 
are lower energy transitions possible for the outer electrons in the vapor atoms 
than in helium, and these outer electrons are the only ones which produce 
Cerenkov radiation which is not absorbed. (The absorption coefficient for 
the substance is small for frequencies below the frequency of the first tran- 
sition but becomes large for frequencies above that.) Therefore, since there 
are relatively few weakly-bound eleetrons, little energy is lost in the form 
of Cerenkov radiation. In each of the mixtures used in the present experiment, 
the Cerenkov correction was neglected when drawing the curves in the figures. 

The expression derived by BUDINI for the behavior of the energy loss by 
ionization beyond minimum energy loss is as follows: 


OE ae Ea Le Are 
| E ne Da stro si 


For the function (6) two possible cases were considered by BUDINI: 


1) The absorption bands for ionization may be considered to be ex- 
tremely broadened lines so that w,/g; & 1. w, is the frequency of the ith absorp- 
tion and g, is the half-width of the line. In this case F(B) may be written in 
the form 

1 ha Ne 
[PAR =1— p+ BEA), 


where (4,)} ~1 and a«=(4aNe?/m). 
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2) The absorption bands are equivalent to moderately broadened absorp- 
tion lines so that (w;/g,)? > 1. In this case 


hia? ~ |} 


[4 II ? 


LENG) E(B?) pt 
where 1 

The main difference between case 1 and 2 is the manner in which the energy 
loss by ionization approaches saturation when (1—?) becomes smaller than 
B?(h2a/I?)(A,,)'. In case 1 the approach to saturation is more gradual than 
in case 2. 

The above expressions have been used to calculate the theoretical curves 
for the present experiment. Two sets of values for the y, as given by STERN- 
HEIMER (°) (1952,1956) have been used to calculate the average ionization 
potential for the mixtures and also to calculate the density effect corrections 
used by STERNHEIMER. 

The maximum energy, 7, was determined by a method of successive ap- 
proximations using Sternheimer’s theory. A value for 7 was assumed and 
used to calculate the amout of energy loss per centimeter at minimum ioniza- 
tion. Then, knowing the number of drops prduced per centimeter at minimum 
ionization, the amount of energy loss 
per drop was calculated. This number | 
was multiplied by the number of “| oropy,, SULLA 5 
drops in a blob (40) minus the ave- 28} 
rage number of drops that would be i 
in the blob due to collisions other 


e ELECTRONS 


than the high energy one. This gave a ap MESONS 2 
new estimate of 7 which was used for ai E 
a new estimate of the energy loss per 29, la 
drop. This procedure was repeated un- | 7 fe E, 5 RE; 


Siren wer) obtamieed Ebevwiscn Fig. 1.— Ionization loss in helium. The theo- 
the assumed and the calculated one. retieal curves are normalized to the u-me- 
The vapor pressure of the liquid sons with By < 30, and are calculated for 
used in the cloud chamber was calcu- different values of the average ionization 
lated from tables in the International Potential for the mixture of gases in the 
Da cloud chamber. Standard deviations are 
Critical Tables. The values of the indicated for a few of the experimental 
various numbers used for the diffe- points. 
rent experiments are in Table I. 

Fig. 1 shows the results of the helium experiment. Each point on the 
graph represents an average of ten tracks. The two lower curves were calcu- 
lated from the theories of Budini and Sternheimer using the ionization poten- 
tials as listed by STERNHEIMER in 1952. The upper curves are for the same 
theories using the higher ionization potentials recommended by STERNHEIMER 


i 6 - Il Nuovo Cimento. 


82 R. G. KEPLER, C. A. D'ANDLAU, W. B. FRETTER and L. F. HANSEN 
TaBLe I. — Cloud chamber operating conditions. 
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| (eV) | I (eV) average Operating | 
: Rae ne ‘all Noble gas | 
Exper; i | maximum |ionization potentia SR SZ tempe- 
aa | energy for the | Pp CH ( rature 
| transfer gas mixture | O) (°C) 
rc ib dee = PE He ee = == | 
| | | 
a { Upper curves | 740 49.4 | OT Woe NT) 
1 | Lower curves | 740 35.7 | 
| Argon 710 | 146 17.5 23.5 
Argon + Helium 700 129.5 Argon 16.1 21.9 
Helium 16.7 
Xenon + Helium | 960 277 | Xenon 5.8 22.9 
| | Helium 6.8 | 
| | 
| i “ 


in 1956. Budini’s approximation I for \/F(f) was used. The theoretical curves 
were normalized to the u-mesons with fy < 30. It is seen that the agreement 
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Fig. 2. — Ionization loss in argon. 
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Fig. 3. — Ionization loss in argon-helium 


mixture. 


between theory and experiment is ex- 
cellent when the higher ionization po- 
tentials are used. 

Fig. 2 shows the results of the argon 
experiment. The two shown curves 
were calculated using the old ioniza- 
tion potentials and were normalized to 
the u-mesons with fy< 30. When 
the new (higher) ionization potentials 
are used, the curves are raised so that 
at fy=10000, the theoretical ioniza- 
tion is changed form 1.69 to about 
1.72. Each point on the graph is again 
the average ionization of ten tracks. 
It is seen that there is an indication 
that the results disagree with theory, 
but the experimental errors are quite 
large because of the poor quality of the 
pictures obtained with argon without 
helium. 

The plus results 
shown in Fig. 3 have lower experi- 
mental errors because of the better 
quality of the pictures. Normalizat- 


argon helium 
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ion to the u-mesons with fy below 30 has the effect of causing almost all of the 
other experimental points to fall below the curve. If the theoretical curve is 
lowered by one standard deviation of an individual point the points above 
By = 30 fit quite well, but all four points below 8y=30 are above the curve. 
The best fit, as with pure argon, and with xenon plus helium, is to a curve 
with a slope smaller than that predicted by the theory. 

The xenon experiment results are shown in Fig. 4. On this graph each 
point represents the average of 20 tracks. Again the curves plotted are those 
calculated using the ionization potentials and normalized to the u-mesons 
with By < 30. 

In the comparison between theory and experiment so far, the theoretical 
expressions for energy loss have been used only to determine the relative va- 
riation of ionization with velocity, and no consideration has been given to 
the absolute energy loss. In order to 


estimate the absolute energy loss, one 30 

has to know the average energy loss slm BUDINI'S ine 
per ion pair and assume that it is inde- A ue THEORY ERO 
pendent of velocity and one has to | ip 
know the condensation efficiency in the ™f do 
cloud chamber during the experiment. 22 ie 
The energy loss per ion pair is inde- 20} 65 + ELECTRONS 
pendent of velocity over quite wide | Lx NE, 1 
velocity ranges in many gases (1216), so Al i + 
this assumption may be justified. Diffi- | + E i 


culty is encountered, however, in esti- Fig. 4. — Ionization loss in xenon-helium 
mating the energy loss per ion pair for mixture. 

a particular gas mixture in the cloud 

chamber. Even though this quantity has been measured for many gases, 
it is well known that the gas has a marked effect on the value obtained (!). 
For the present experiments it was assumed that the energy loss per ion pair 
in the vapor and in helium when mixed with the cloud chamber vapors was 
30 eV. The values used for argon and xenon were 25.8 and 22 eV respec- 
tively (**). 

The condensation efficiency was estimated by counting tracks in which 
the positive and negative ions had been separated by the clearing field. The 
values obtained for the ratio of the number of drops on negative ions to the 
number of drops on positive ions were 0.4, 0.4 and 0.2 for the helium, argon 


J. Weiss and W. BERNSTEIN: Phys. Rev., 103, 1253 (1956). 
W. P. Jesse and J. SADAUSKIS: Phys. Rev., 97, 1668 (1955); 90, 1120 (1953). 
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and xenon experiments respectively. It was assumed that for the above ratios 
the condensation efficiency on the positive ions was essentially 100 per cent (1°). 

The energy loss at minimum determined experimentally for the helium 
experiment was 415 eV/cm compared to the theoretical value 392 eV/cm. 
For the argon experiment the experimental and theoretical values were 430 
and 433 eV/em respectively, and for the xenon experiment they were 320 
and 425 eV/cm. 

The agreement between theory and experiment for the helium and argon 
experiments is excellent, but the discrepancy in xenon led us to perform ano- 
ther experiment in xenon under better controlled conditions. A special counter 
triggering system was used to select u-mesons in the range of fy from 2.4 
to 4.9, so that u-mesons in the range near minimum were observed. The mag- 
netic field was cut off, and grids were placed at the sides of the chamber 
to provide an electric field to separate the positive and negative ion columns. 
Under these conditions it was possible to count the droplets formed in the 
positive and negative column independently, and to determine the ratio of 
these as a function of expansion ratio. 

We found no significant variation of the number of drops in the positive 
column when the ratio of the number in the negative column to the number 
in the positive column varied from 0.15 to 0.4. We, therefore, assume that 
even when only 15 per cent of the ions in the negative column form drops, 
practically all the positive ions form drops. From this we can calculate the 
number of ions formed, and assuming the above values for energy loss per 
ion pair, can compare our experimental value with the theoretical value. For 
our particular conditions, the experimental value was (28.9-+-0.6) drops cm. 
Theoretically the figures are (31.8 + 0.4) drops em for BAKKER-SEGRE (2°) 
ionization potentials (STERNHEIMER,1952) and (29.6 + 0.4) drops/cm using the 
CALDWELL (2!) (STERNHEIMER, 1956) ionization potentials. The errors came 
from uncertainties in the energy loss per ion pair. We believe that the experi- 
mental and theoretical results are in agreement within the errors. The discre- 
pancy in xenon was probably due to a bias in selection of doubled tracks where 
the conditions were not good for counting the negative column. 


9. — Discussion of results. 


The experimental results demonstrate clearly that the slope of the rela- 
tivistic rise in the heavier gases is considerably smaller than that predicted 
theoretically, and seems to be approximately independent of the atomic num- 
ber. The plateau in helium is, however, substantially lower than the plateau 


(20) C. J. BAKKER and E. SEGRÈ: Phys. Rev., 81, 489 (1951). 
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in argon or xenon. Previous experimental results have not been accurate 
enough or have not reached high enough velocities to have found this effect. 
We have explored various reasons for the difference between theory and expe- 
riment which may be summarized as follows: 


1) The value of 7, the.maximum allowable energy transfer, may depend 
on the atomic number. Thus 7, instead of being a fixed value depending only 
on the number of drops in a blob just countable, should perhaps be corrected 
for the binding energy of the struck electron. The actual energy transfer in 
the case of an electron in the K-shell of xenon, when a blob of 40 drops is for- 
med, is approximately 40 KeV, rather than the assumed value 7—960 eV. 
A calculation was made, using Sternheimer’s theory, but averaging over the 
various values of 7 for the different electron shells, to determine the effect 
on the slope of the curve. The slope was indeed reduced, but not sufficiently. 
For example in xenon, the value of the rise of energy loss above the minimum 
at fy=100 was lowered from 1.48 to 1.45, compared with the experimental 
value of 1.38. It seems clear, however, that this effect is valid and important 
and should be included when calculations of energy loss are made. 


2) There has been discussion of the possibility that J, the average ioni- 
zation potential, decreases with increasing velocity (21). Measurements are 
available, however, only as far as fy=0.5, and there are no experiments on 
dependence on fy at the very high velocities with which we are concerned. 
However, if J decreases as Py increases, the slope of the curve would be greater 
in contrast to the experimental result. 


If 7 and J are independent of energy, the value of 7/1? determines the 
slope of the relativistic rise, and together with the energy loss per ion pair, 
the absolute number of ion pairs per centimeter at minimum. Since the theory 
and experiment agree well up to values of Sy near the minimum, it is difficult 
to reconcile theory with experiment beyond the minimum simply by changing 
the value of 7/12. At minimum 7/1*=(960/(277)?)=0.013. In order to fit 
the slope of the rising part of the curve, 7/1? must be 2.73. We conclude that 
this ratio has no theoretical significance in the relativistic region. 

Other possibilities have been considered but no detailed theoretical calcu- 
lations have been made. Various possibilities are: 


1) The ratio of energy loss to excitation and to ionization may not be 
independent of velocity and this dependence may vary with atomic number. 


2) Shielding effects of the electron shells in the heavy elements may 
be important. The impact parameters in the region of the relativistic rise 
are orders of magnitude greater than those below the minimum of ionization, 


(21) D. O. CALDWELL: Phys. Rev., 100, 291 (1955). 
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where the ionizing particle passes directly through the atom. Complete shiel- 
ding of the K and L shells in xenon in the relativistic region produces a large 
effect (in the right direction) on the slope of the relativistic rise. 


It seems clear that further calculations of ionization loss in heavy gases 
in the relativistic range should be made (*). Although the smaller slope will 
reduce the precision of velocity determination in the high momentum region, 
it seems probable that particle identification in this region may be made under 
sufficiently controlled conditions. 


The authors are indebted to T. A. AGGSON and R. W. BIERCE who helped 
with the drop counting and data analysis. We also thank Dr. R. M. STERN- 
HEIMER who sent us the results of his calculations for our helium experiment 
based on his theory and Prof. P. BupINI for valuable discussions of the re- 
sults. Prof. M. A. RUDERMAN participated in informative discussions on 
the theoretical implications of the results. 


(*) Progress has been made recently in the theoretical calculations of TIDMAN (22) 
applying to oxygen, and Bupini and TAFFARA (?) have also published the elements 
of a theory which should be applicable. 

(22) D. A. TIDMAN: Nuclear Physics, 2, 289 (1956) and in press. 

(73) P. Bupini and L. TAFFARA: Nuovo Cimento, 4, 23 (1956). 


RIASSUNTO (*) 


In elio, argon e xenon si sono eseguite misure della ionizzazione prodotta da particelle 
relativistiche cariche. Nel caso dell’elio i risultati sperimentali si accordano bene con 
la teoria. Nell’argon e particolarmente nello xenon l’aumento relativistico della ioniz- 
zazione trovato sperimentalmente è inferiore a quanto la teoria lascia prevedere. Le 
misure furono eseguite col metodo del conteggio delle gocce in una camera a nebbia. 


(*) Traduzione a cura della Redazione. 
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Absorption and Refraction of Some Polar Gases 
as a Function of Pressure at Microwave Frequencies (*). 


A. BATTAGLIA, F. BRUIN (+) and A. GozzINI 


Istituto di Fisica dell Universita - Pisa 


(ricevuto il 9 Settembre 1957) 


Summary. — Absorption and refraction of several polar gases has been 
measured, as a function of the pressure, at 32 mm wavelength and 0 °C. 
For ammonia, trimetilamine, ammonia-argon and ammonia-hydrogen 
mixtures a quadratic dependence has been found. For OCS and several 
amines the absorption increases less than quadratically. At the same 
pressure and concentration, we found the same absorption for ammonia- 
argon and ammonia-hydrogen mixtures. For absorbing gases, at low 
pressures, the refractive index is not proportional to the pressure. 


1. — Theory. 


The spectral absorption and refraction for a single transition in gases at 
microwave frequencies is usually described by the so called Van Vleck and 
Weisskopf formula (') for the absorption coefficient « and refractive index n 
given by (1): 


82? Nf f 1/27T 1/27Tt | 
(1) = | Mi; (Pv? Fei 2 Se er ou (0 
3ckT |(@ = v0)? + (nt) © (v + mn) + (1/27) f 
SIN eee 
(2) n—1 » 2(n—1) = SET bes | ara 


Vv | tesi Ma Lee ns Be fe ) Ls 
So Gv + enr) (+ ¥9)* + (1/200)8)] À? 
(*) The research reported in this paper has been made possibile by support ex- 
tended by C.N.R. and Centro Microonde. 
(+) Now at the Natuurkundig Laboratorium, Universiteit van Amsterdam. 
(1) J. H. VAN VLECcK and V. F. Weisskorr: Rev. Mod. Phys., 17, 227 (1945). 
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in which « = absorption coefficient, in cm, 

N — number of absorbing molecules, per unit volume, 

f — fraction of molecules in the lower state 1, 

Li; = dipole moment matrix element of states 7, 7, 

yy — resonance frequency of the transition between states à, j, 

T — mean time between two molecular collisions, 

e — 3:10 cm/s, &—1.38:10-16 erg/K° and T absolute temperature, 
and 7 represents the contribution to the refractive index of the (1— f)N 
molecules that are not absorbing, being not in the lower state 2. 


In the far wing of a line, at moderate pressures, the quantity (1/277)*< 
< (Vv), and we may bring formula (1) into a more practical form: 


1+ (v/v)? (Ay | 
= 2max Re ra i 
(3) Xing Xmax (Vo/v = v/v0)? | Vo Fi 


a. being the absorption coefficient at the resonance frequency 


max 


(~ maximum absorption). 


Ay the line half width constant Av=1/277p, and p the gas pressure. The 
central absorption «,,. is constant, if 1/t is proportional to the gas pressure. 
Therefore, at a fixed frequency, the absorption in the far wing of a spectral 
line, as a function of pressure, should be a quadratic parabola. 

At the same pressure and frequency conditions, formula (2) can be put 
in the form 


(4) 2(n ga 18) Ri max aioe 
IT 


CP PA Î JÉ nil rbt Mano 

Va — v? \( 

x being proportional to the pressure, formula (4) predicts, that at low pressures, 
(n — 1)/p increases or decreases with pressure according as v2. 

At higher pressures, formula (2) shows that also the contribution of the 
resonant molecules tends to become proportional to the pressure. Therefore 
at higher pressures the Clausius-Mossotti law (n —1)/p = const holds. 

In the following both formulae (3) and (4) have been applied to measu- 
rements on a number of dipole gases. 
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2. — Measurements. 


The apparatus used in the measurements discussed below was a slightly 
simplified version of the one previously described (2). The resonance frequencies 
of two equal adjustable cavities, each receiving half the klystron power, were 
compared by means of a sensitive pulse technique. Most of the measurements 
were performed at 0 °C. For this reason the cavities have been put in a bath 
of melting ice. 


21. Absorption of NH, and NH, mixtures. — As was shown in the first 
section, the absorption in the low wing of a spectral line, for a fixed fre- 
quency », as a function of pres- 
sure, according to the simple 
theory should be a parabola. NH3 

À 


In Fig. 1 two parabolae of the CEE 
type represented by formula eu ie Es if 
(3) are given. The upper curve SELE, 5 rd 
shows the absorption of NH, TAO 
as a function of pressure at a ie 
wavelength of 32mm and at 
a temperature of 0°C. The sol 
lower curve is obtained from a 
similar measurement at 15 °C. 
As the centre of the absorption 
band of NH, is situated at a 
wavelength of 12.5 mm and 
measurements were performed 0 100 
at a wavelength of 32 mm, the Fig. 1. - Absorption of NH; versus pressure at 
quantity 2(1+(v /vo)?) I(vo/v—v/vo)? temperatuers of respectively 0 and 15 °C. 

in (3) for this case is equal to 

0.528. The maximum absorption at a pressure of 100 mm and at room tem- 
perature is according to BLEANEY and PENROSE (*) equal to 49-10-4 em, 
so that the lower curve of Fig. 1 should be 


a-10° cm-! 


p (mmHg) 
1 | 


i 


A \2 
(5) Kana = 2582 10 (~») em. 


0 


At pressures below 150 mm Hg, the curves fitted to the observed points are 
given by 
(6) ca = Oboe for the measurements at 0°C, 


(2) A. Barracria, F. BRUIN and A. Gozzini: Nuovo Cimento, 
(3) B. BLEANEY and R. P. PENROSE: Proc. Phys. Soc., 59, 418 (1947). 
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and 
(7) x=0.433-10-8p?, for the measurements at 15 °C. 

It follows from (5) and (7) that the line width constant should be 
Ay = 31 MHz per mm Hg. 

In order to obtain the theoretical value of Av one should average over all 
rotational quantum numbers J and Æ of ammonia. If this is done one 
finds (*°) for Av about 26 MHz per mm Hg. This means that in the case 

considered here the observed ab- 
107° sorption is 40% higher than 
expected. This discrepancy has 
been noted before (5). 

A second fact of interest which 
may be deduced from Fig. 1 is 
the dependence of the absorption 
on the temperature. From (6) and 
(7) we deduce a 7-* dependence. 
This agrees with the earlier and 
more detailed measurements per- 
formed by the group of N.B.S. 
Washinton (°’) whereas other 
authors report a T° depen- 
dence (°°). No great value should 
be attached however, to this re- 
sult, because the temperature of 
15 °C has not been determined 
with great accuracy. 

In order to check if the curves 
of Fig. 1 are indeed quadratic 
parabolae, in Fig. 2 we have plot- 
ted the absorption as a function 

20 10° 10° of pressure on a twofold loga- 
Fig. 2. — Absorption of pure NH, gas and of ANNE Si EG UE O 
mixtures of NH; and argon as a function of Pendence, the result should be a 

pressure on a twofold logarithmic scale. straight line having a slope of 2. 


p(mmHg) 
SSS 
10 ASA RER MERE Ra ere iI) 


(4) B. BLEANEY and J. H. N. Lousser: Proc. Phys. Soc., 63, 483 (1950). 

(5) C. H. Townes and A. L. ScnawLow: Microwave Spectroscopy (New York, 
1955), p. 347. 

(5) G. BIRNBAUM and A. A. Maryorr: Journ. Chem. Phys., 21, 1774 (1953). 

(?) G. BIRNBAUM and A. A. MARYOTT: Journ. Chem. Phys., 22, 1457 (1954). 

(8) H. KrIsHAJI and P. Swarup: Journ. Chem. Phys., 22, 1456 (1954). 

(*) P. Swarup: Journ. Scient. Ind. Res., 13 B, 389 (1954). 
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A number of diagrams of this kind is shown in Fig. 2. The upper curve in 
Fig. 2 represents the absorption of pure ammonia at 0 °C and at a wavelength 
of 32 mm. Up to a pressure of 200 mm Hg the curve is a straight line having 
a slope of 1.95, whereas at pressures below 100 mm Hg the slope is 2 within 
the limits of experimental error. At higher pressures the slope decreases, but 
this is not surprising as for increasing pressure the half width does not re- 
main proportional to pressure (1°). 

In Fig. 2 is also shown the absorption as a function of pressure of mixtures 
of ammonia and argon for different fractions of NH,. The curves all have 
the same slope, which simply means that we are each time in the far wing 
of an absorption line. A set of curves like in Fig. 2 may be used to determine 
the absorption as a function of concentration of NH; for some fixed pressure, 
after which the collision cross- 
section may be calculated. 1 
Such derived curves are shown 


© MIXTURE NH,-A 


in Fig. 3 for mixtures of am- + MIXTURE NH,~-H> 
| À n — UE Goa 
monia and argon and of am- NES CODES ORNE er Site 
‘ ; SZ A= 32mm 
monia and hydrogen. It is seen È 


TAZIONE 


that the results for both mixtu- P = 200 mm Hg 


res coincide. 

The analytic expression for 05 
the curve of Fig. 3 may be 
easily derived. Having a mix- 
ture of Ne molecules of am- 
monia and N(1 —c) molecules 
of argon, if v, and v, are the 


monia-ammonia and ammonia- 


mean relative velocities of am- a 
Le 


0 50 100 
argon molecules, and o,, o, the 
cross-sections for ammonia-am- Fig. 3. — Relative absorption of NH,-A and 
monia and ammonia-argon col- NH,-H, mixtures as a function of concentration 


7 | i Le h SS r 2 û Le 
lisions, the collision rate for for a pressure of 200 mm Hg 


an ammonia molecule in the 
mixture is 1/7 = Nevo, + N(1— c)v;o,. Putting a = v,0,/v,0, the absorption 
relative to pure ammonia, as function of the concentration c may be easily 


shown to be: 


y el a) ca. 


The curve of Fig. 3 is well represented, at concentrations below 80% by 
y = 0.89¢c?-+-0.11¢, therefore we obtain for a the value of 0.11. 


(9) G. BirxBAUM and A. A. Maryorr: Phys. Rev., 92, 270 (1953). 
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This is the same than that found for H,, but higher than that found for 
A at resonance (11). 


22. Absorption of some amines and of carbonil sulfide. — In Fig. 4 the ab- 
sorption versus pressure curve of OCS and of NH,(CH,); NH(CHs;).; N(CHs)3; 
NH,(CH;); NH(C,H;), are plotted in a twofold logarithmic scale. All measure- 
ments of figure are at 0°C and 
32mm wavelength. The microwave 
spectrum of amines is little known. 
At low pressures, only methyl- and 
ethylamine have been investigated 
and many lines resolved in the fre- 
quency region between 15000 and 
30000 MHz (*"). At high pres- 
sures, HERSHBERGER (15) studied 
the absorption of amines, at 12 mm 
wavelength and at atmospheric 
© pressure. ; 

Our investigation gave the 
following results. The absorption 
of the dimethyl-ethyl-, and diethyl- 
amine, at 32mm wavelength is 
SI found to be proportionalto pressure; 


È - De the value of the coefficient of pro- 
gid: = 1:75 
tg y =1.5 portionality is 5-10’ cm-!/mm Hg, 


for the dimethylamine and the 
L ethylamine; 1.7-10-? cem-!/mm Hg 
i ao for the diethylamine. The pres- 
sure interval for which this linear 


p (mm Hg) 


Fig. 4. — Absorption versus pressure for 
several gases, at 32 mm wavelength and law holds, is (30-300) mm Hg for 
at 0°C. the first two substances, and 


(30 = 150) mm Hg for the diethyl- 
amine. Above these pressures, the substances (at 0 °C) begin to condense. 
At 12mm wavelength, and at atmospheric pressure, HERSHBERGER (15) found 
that the ratio between the values of the absorption coefficients of the di- 
methil- and ethylamine is 7 to 11. According to our measurements, the two 


(41) C. H. Townes and A. L. ScHawLoWw: Microwave Spectroscopy (New York, 
1955), p. 364: 


M. MATRICON and J. BONNET: Journ. Phys. et Rad., 15, 647 (1954). 


) 
18) K. Sarmopa and J. NiscHiLawa: Journ. Phys. Jap., 9, 974 1(954). 
) 
) W. D. HERSHBERGER: Journ. Appl. Phys., 17, 495 (1946). 
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coefficients have equal values at 32 mm wavelength, over the entire pressure 
interval investigated. The dimethylamine was studied at 12.5 mm wavelength 
also. We found that at this frequency the absorption coefficient depends 
on pressures less markedly than at 32 mm, while its value is appreciably hi- 
gher (Fig. 5). 


| TE RC NH-(CH 3) 
o le0c¢ 
o) 
100 3 Ocs 


Usa | usu | 


p(mm Hg) 


D e eee el tI se Che arene) 
1 5 10 100 500 1000 


Fig. 5. — Absorption versus pressure for OCS and NH(CH;), at 12.5 mm wavelength 
and 0 °C. 


The absorption of the methylamine, at 3.2 cm wavelength, is roughly pro- 
portional to the pressure, for pressures ranging from 400 to 700 mm Hg; at 
lower pressures, the increase is more steep. 

Although the absorption of N(CH;); is seven times weaker than that of 
ammonia the dependence on pressure is about the same. This suggests that 
here also we are in the wing of an absorption band. 

The absorption of carbonil sulfide, at 32 mm wavelength and 0 °C, is 
plotted against pressure in Fig. 4. Fig. 5 shows the results of similar measu- 
rements in the 12 mm region, at a frequency slightly below the frequency of 
the J= 1, J— 2 rotational transition. The first rotational line of this linear 
molecule is predicted to fall at 24.6 mm wavelength for the most abundant 
isotopic species 1012C%#S, so that, at 32 mm, we are in the far low wing of 
the resonance. From 300 mm Hg to atmospheric pressure, absorption is found 
to be proportional to the 3 power of the pressure. Below 300 mm, the de- 
pendence is more pronounced and it approaches the quadratic law, but the 
absorption is too small to be measured with sufficient accuracy by our appa- 
ratus, at the time when the measurements were performed. The absorption 
measured at a frequency very close to that of the J=1, J=2 transition 
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- 1 raises quadratically with the pres- 

Na sure at very low pressures. Accord- 
Ca =) i ing to the theory of collision broad- 
ge ening at higher pressures one ex- 

DUT i pects the absorption to be indepen- 
E dent of the pressure. However we 
à i observe a slight but definite increase 
of absorption. Both results of mea- 
surements at 32 and 12.5 mm in- 
dicate that the line width is less 


Pa than proportional to the pressure. 
——0——_o- sE = GS 


(n-1) 10° 


2°3. Refraction of NH;, OCS, CO, 
and A. — The Van Vleck and Weiss- 
i  kopf formula predicts that the re- 
fractive index of a gas in the wing of 
UE a line should not be, at low pressures, 
p(mmHg) = proportional to the pressure. This has 
0 100 00 0 200 I been confirmed. The results are re- 
Fig. 6. — (n—1)/p versus pressure for presented in Fig. 6 where (n — 1)/p 
several gases. is plotted versus pressure for several 
gases. For non absorbing gases as 
argon ad CO, the refractive index is strictly proportional to the pressure, whe- 
reas for ammonia and COS the curves show the behaviour predicted at low 
pressures by formula (5). Unfortunately with the klystrons at our disposal it 
was not possible to make measurements at frequencies in the high wing of 
the resonances, where formula (5) predicts that at low pressures, (n — 1)/p 
should increase with increasing pressure. 


RIASSUNTO 


Si riferiscono i risultati di misure sul coefficiente di assorbimento e sull’indice di 
rifrazione di diverse sostanze, alle lunghezze d’onda di 32 e 12 mm. Tali grandezze 
sono state misurate in funzione della pressione a frequenza e temperatura costante. 
A 32 mm di lunghezza d’onda il coefficiente di assorbimento dell’ammoniaca, di miscugli 
di ammoniaca ed argon ed ammoniaca ed idrogeno e della trimetilammina è propor- 
zionale al quadrato della pressione. Alla stessa pressione e concentrazione di ammo- 
niaca, il coefficiente di assorbimento di miscugli di ammoniaca con argon e di ammo- 
niaca con argon e di ammoniaca con idrogeno è lo stesso. Per NH; e OCS, alle basse 
pressioni, l’indice di rifrazione non è proporzionale alla pressione. 
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IL NUOVO CIMENTO VO VITAE NET 1° Gennaio 1958 


The Interaction between Equilibrium Defects 
in the Alkali Halides: the « Ground State» Binding Energy 
of the Vacancy Pair (*). 


NÉE ADO STE TRE C METAL 


Istituto di Fisica dell Università - Palermo 


(ricevuto il 10 Settembre 1957) 


Summary. — The Born-Mayer model of ionic solids is used to calculate 
the binding energy of a positive-ion vacancy and a negative-ion vacancy 
at the distance of nearest neighbours in NaCl and KCl crystals: the 
results are 0.6 eV and 0.7 eV respectively. This leads one to expect in 
NaCl an equilibrium concentration of vacancy pairs ranging from 10% 
to 20% of the equilibrium concentration of Schottky defects between 
500 °C and 700°C. The available experimental data are compatible with 
the existence in the alkali halide crystals of an equilibrium concentration 
of vacancy pairs of this order of magnitude in the intrinsic temperature 
region. 


1. — Introduction. 


Morr and GuRNEY (:) were apparently the first to point out that positive- 
ion vacancies and negative-ion vacancies in an alkali halide crystal will tend to 
combine in pairs and SEITZ (?) emphasized shortly afterwards that these pairs, 
and the larger vacancy agglomerates, will play an important role in the physics 
of these crystals. To a rough approximation, the binding energy of two oppo- 
sitely charged vacancies at the closest lattice distance of approach, the binding 
energy of the ground state of the vacancy pair, may be taken equal to the 


(*) The first paper of this series was published in Nuovo Cimento, 11, 274 (1954). 
(1) N. F. Morr and R. W. Gurney: Electronic Processes in Ionic Crystals (Oxford, 
940), p. 41. 

(2) F. SEITZ: Rev. Mod. Phys., 18, 384 (1946). 
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electrostatic interaction — e?/Kr, of the effective charges of the two vacancies 
(K, static dielectric constant of the crystal; x, nearest-neighbour distance). 
It is obvious that an estimate of this sort of the binding energy cannot give 
an accurate value for the equilibrium concentration of vacaney pairs at a 
given temperature. In effect this equilibrium concentration is extremely 
sensitive to the value of the binding energy: for a NaCl erystal at 800 °K the 
ratio of the density m of vacaney pairs to the density n of Schottky defects 
(ref. (2), eqn. (26)) with (*)n=1.2-10?* exp [— 11700/7]cm-* is 7.6 if the binding 
energy is Coulombic (0.91 eV), and only 0.085 if the binding energy is 0.6 eV. 
This estimate indicates clearly the necessity of having an accurate value for 
the binding energy of vacancy pairs to draw any conclusion on their equili- 
brium concentration at different temperatures. 

We have computed the binding energy of the ground state of the vacaney 
pair in NaCl and KCI crystals using the Born-Mayer model as first applied 
to the study of point imperfections in the alkali halides by Morr and LITTLE- 
TON (*). We find that the binding energy is smaller than the Coulombic value, 
as it is to be expected (5), and by several tenths of an electron volt. A simi- 
lar, but less accurate, calculation had previously been performed for NaCl 
crystals by REITZ and GAMMEL (°). 


2. — The binding energy of the ground state of the vacancy pair. 


The binding energy of the ground state of the vacancy pair is equal to 
the difference between the work necessary to create a positive-ion vacancy 
in the perfect lattice, and the work necessary to create it in the position of 
nearest neighbour of a negative-ion vacancy. To take into account the elec- 
tronic and displacement polarization of the lattice around the vacancy created, 
one writes the work of removal of the positive ion as the negative of the ave- 
rage of the potential energies in the position of the ion before and after it has 
been removed. 

The important contribution to the potential energy of an ion in an alkali 
halide crystal are the electrostatic and repulsive energies, and these are the 
only terms that we shall consider. For the repulsive energy we adopt the 
exponential form proposed by Born and MAYER (75). 


(3) H. W. ErzeL and R. J. MAURER: Journ. Chem. Phys., 18, 1003 (1950). 
(4) N. F. Morr and M. J. LirtLETON: Trans. Faraday Soc., 34, 485 (1938). 
(3) F. G. Fumi and M. P. Tosi: Faraday Soc. Discussions, 23, 92 (1957). 
(6) J. R. Rerrz and J. L. GAMMEL: Journ. Ohem. Phys., 19, 894 (1951). 
) M. Born and J. MAyYER: Zeits. f. Phys., 75, 1 (1932); J. MAYER and H. HELM- 
HOLZ: Zetts. 7. Phys., 75,19 (1932): 
(8) F. Bassant amd F. G. Fumi: Nuovo Cimento, 11, 274 (1954). 
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3. — The work to create a positive-ion vacancy in a perfect crystal. 


The MOTT and LITTLETON type of procedure to calculate the work to create 
a positive-ion vacancy in a perfect NaCl or KCl crystal has been described 
in detail in a previous paper (*°). For the constants entering the calculation 
we use the same numerical values quoted in ref. (8), Table I except that we 
adopt for the electronic polarizabilities «, and «_ of cations and anions the 
values reported by TESSMAN, KAHN and 


SHOCKLEY (TKS) (#°) for ions in crystal 210 Ta 
lattices, instead of the values given by | RS KCl 
PAULING (11) for free ions. This entails some E) 

î 


numerical changes in the calculation pre- 
viously described (*), and the new values for 
the work to create a positive-ion vacancy 
(4.78 eV in NaCl; 4.44 eV in KCl) are found 
to be somewhat larger, as the TKS polariza- 10-4L 
bility for Cl” ions is smaller than the Pauling 
polarizability. The net increase of the energy 
of formation of a Schottky defect caused by 


0.076 0.104 


the change in polarizabilities is fairly small 007 009 -E+ 01 
(0.09 eV in NaCl, 0.04 eV in KCl) as the ——— NaCl —:—--— KCl 
work to create a negative-ion vacancy is (Nat) 
NaCl) = > n 263° 10-* = 0.008 
somewhat reduced. DO N one 
(Kt) n” 029 
a E 90- = 9 
m(KCT) = er (EC) 1.190-1 0.922 


4. — The work to create a positive-ion vacancy 
near a negative-ion vacancy. Fig. 1. — Displacement and elec- 
tronic polarization of the imme- 
To calculate the potential energy in the diate neighbours of a negative-ion 
position of a positive ion which is nearest vacancy in NaCl and KCl. 
neighbour of a negative-ion vacancy we need 
to determine the polarization of the lattice due to the presence of the nega- 
tive-ion vacancy. Of course the procedure is entirely analogous to that described 
for the positive-ion vacancy in a previous paper (*): the determination of the 
displacement ér, and of the electronic dipole mer, of the immediate neighbours 
of the vacancy is illustrated in Fig. 1 (??). The electrostatic potential energy 
(°) Equation (7) of ref. (8) contains a misprint. The first term inside round brak- 
kets should read 6£/(1+£), and not 6/(14-£). 
(0) J. R. TESsMAN, A. H. Kann and W. SHocKkLEY: Phys. Rev., 92, 890 (1953). 
(11) L. PAULING: Proc. Roy. Soc., A 114, 181 (1927). 
(12) The equations at the bottom of Fig. 1, ref. (8) contain a misprint. A factor 1074 
was omitted from the numerical values of F, introduced in the equations: however, 
the final values quoted for m(NaCl) and m(KCl) are correct. 


Pa 7 - Il Nuovo Cimento. 
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In the formula for AE° the first group of terms gives the change in electro- 
static interaction energy amongst the ten nearest neighbours of the two vacan- 
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cies, while the second (and third) group of terms give the change in electrostatic 
interaction energy between the ions 1, 10 (and the ions from 2 to 9) and the 
rest of the lattice. In the formula for AE” the first group of terms refers to 
the ions 1, 10 and the second group of terms refers to the ions from 2 to 9: 
the dipole-dipole interaction has been neglected. The numerical results are 
illustrated in Fig. 4: for AH” 


and AE” there is a direct | PS PEN | 
(e 
28) - AE | 
— —-—-— KCI 
| | 


correspondence between the 
values for NaCl and KCl as 
7 AE® and (ri/(a, +a«_)) AH” | 
are functions of 7 which de- ti Wee 2% 
pend only on the type of : 
structure and on the charge 
of ions. For NaCl, our nu- 
merical results agree with 
the results given by REITZ 
and GAMMEL (ref. (5), Fig. 4) 
but for KCl they show sub- 
stantial deviations from the 
results reported by DIENES 
(tet. (©), Fig. 6)2-for AE? 
this is to be imputed to the 
inclusion of the next-nearest 
repulsion which had been 
neglected by DIENES, but 
for AH’ the discrepancy is 
unexplained. The values that 
we find for the displacement 
no of the immediate neigh- 
bours of the vacancy pair 
appear entirely reasonable. 
Indeed, if one considers that 
the nearest neighbours of AIA ne dl Se ea 
each of the two vacancies are 708 005 0067) 0083) 010 
pushed out by the adjacent pis. 4. — Determination of the displacement mr, of 
vacaney and pulled in by the immediate neighbours of a vacancy pair. 
the other vacancy, it is clear 

that the uniform outward displacement of the ten nearest neighbourg of the 
vacancy pair cannot be as large as the larger of the displacements around 
the isolated vacancies. This condition is verified by our values for the displa- 
cement nr, in NaCl and KCl crystals, and by the value for 7” reported by 
Reitz and GAMMEL for NaCl (6), but it is not verified by the value for 77% 
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reported by DIENES for KCl (1%) (Table I). The fact that Dienes’ distortion 
around the vacancy pair appears too large had already been noted by PRATT 
(>), who followed, however, an unreliable line of reasoning. 


TABLE I. — Displacements of the immediate neighbours of a vacancy pair and of isolated 
vacancies (“). 


2 Le ia EN pil Saia 
| | 
NaCl PAULING 0.068 (?) | 0.104 | 0.083 (°) | 
l= TKS 0.070 | 0.104 0.083 
= | f ws] 
KOS «| PAULING 0.065 (@) 0.077 (0.10) (4) 
TKS 0.065 | 0.076 | 0.067 


(a) The quantity En (É r) gives the outward displacement of the nearest neighbours of 
a positive-ion (negative-ion) vacancy and the quantity nr the uniform outward displacement 
of the ten nearest neighbours of the vacancy pair. The values reported were all computed with | 
the Born-Mayer repulsive potential: they change slightly if one replaces the Pauling polariz- | 
abilities by the TKS polarizabilities. 
(b) These values were reported by BASSANI and FUMI [ref. (*), Fig. 1]. 
(c) This value was reported by REITZ and GAMMEL (ref. (5), Fig. 4]. 
(d) This value was reported by DIENES [ref. (%), Fig. 6], but it appears questionable. | 


The electrostatic and repulsive potential energies #, and ©, at a positive-ion 
vacancy near a negative-ion vacancy are computed in the same sort of way 
as the energies Æ, and ®,. The relevant formulas are: 


(11) = ’ 


Ga f JE 1 4 
12 DEEE = PI TRES - ; ne 
ta RR | È 


+ 1 


Je “ae ee ee eee 
POET GEE ERNEST RD ES Ae a) 


A+ ml. 
ol 
The polarization term of Æ, must be computed directly as the potential 
energy in the positive ion in question due to the dipoles induced on the ten 
ions which are immediate neighbours of the vacancy pair, and on the second 


uk | Oy 273 70 I 7 
owl [A+ mt + EN À) + exp 


(13) P. L. Pratt: Report of the Conference on Defects in Crystalline Solids held at 
Bristol, in July 1954 (London, 1955), p. 402. 
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and third neighbours of the positive-ion vacancy by the field generated by 
the charges of the two vacancies, and by the displacements of the ten imme- 
diate neighbours of the vacancy pair. The polarization of the distant ions 
is neglected since the vacancy pair is neutral. 

The various contributions to the work to create a positive-ion vacancy 
near a negative-ion vacancy in NaCl and KCl crystals are listed in Table II: 
in NaCl this work is equal to 4.18 eV, while in KCl it is equal to 3.72 eV. The 
numerical results for NaCl show substantial deviations from the results quoted 
by ReITZ and GAMMEL (ref. (°), Table II): for the initial configuration, the main 
source of discrepancy is the difference in the value for the displacement of the 
nearest neighbours of a negative-ion vacancy, while for the final configuration 
there is also an important effect of the replacement of the Pauling polarizabilities 
by the TKS polarizabilities. These effects concur to raise the work to create a 
positive-ion vacancy near a negative-ion vacancy in NaCl from the value 3.86 eV 
quoted by REITZ and GAMMEL to our value. The displacement 0.078 7, that REITZ 
and GAMMEL attribute to the nearest neighbours of a negative-ion vacancy in 
NaCl had been computed by Mort and LITTLETON (*), using a repulsive potential 
rather different from the BORN-MAYER potential that REITZ and GAMMEL use in 
their calculations: the value that REITZ and GAMMEL should have used for 
this displacement is 0.104 7) (Table I). 


TasieE II. — Contributions to the work to create a positive-ion vacancy near a negative-ion 
vacancy (eV). 


ae |a - BE) DRIVE re | Sige Vere (eo 
| NaCl 4,31 De BO A Sr oo 
KCl | 3.75 | 0.51 3:07 294 TETE EN — 193° hoes 


5. — Discussion of the results. 


The values calculated in Sects. 3 and 4 for the work to extract a positive 
ion from a perfect lattice and from the displaced equilibrium position near 
a negative-ion vacancy give at once for the ground-state binding energy of 
a vacancy pair the values 0.60 eV in NaCl and 0.72 eV in KCl, which, of course, 
should not be taken literally, owing to the approximation inherent to this type 
of calculation (ref. (5), sect. 4). The values that we compute are probably some- 
what too large as the approximation involved in adopting with DIENES (7%) an 
equal outward displacement for the nearest neighbours of a vacancy pair reduces 
the work to extract a positive ion near a negative-ion vacancy. On the other 
hand, the slight uncertainty in the value of 77, due to the flat minimum of 
AE in Fig. 4, appears of little consequence for the binding energy of the 
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vacancy pair which is not affected by small changes in the value of 7. Simi- 
larly, the replacement of the TKS polarizabilities by the Pauling polarizabili- 
ties has little effect on the binding energy of the vacancy pair, as the work 
to extract a positive ion from a perfect lattice and from the position ot 
nearest neighbour of a negative-ion vacancy are reduced by comparable 
amounts (1): for NaCl, the binding energy becomes 0.66 eV, as previously 
reported (17), while for KCl it becomes about 0.76 eV, a 

value sizeably different from the one reported previously 

03 (0.85 eV)(17), which was based on Dienes’ distortion 
around the vacancy pair. It should finally be stressed that 

the somewhat larger value of the binding energy in KCI ap- 

01 pears quite reasonable as the outward distortion around a 
THC negative-ion vacancy is appreciably smaller than in NaCl. 

300 600 70 The ratios of the density m of vacancy pairs to the 

Fig. 5. — The «theore- density n of Schottky defects to be expected in NaCl ery- 
tical» ratio of the stals in the intrinsic temperature region on the basis of a 
density m of vacancy simple association theory (ref. (2), eqn. (26)) with a binding 
pairs to the density » È : ; = i 
IS defects in energy of 0.6eV taking n =1.2-102%exp[—11700/T] cm 
NaCl in the intrinsic are reported in Fig. 5. The corresponding values for KCl, 
temperature region. | Which range from 0.19 at 500 °C to 0.82 at 700 °C, are not 

reported in the figure also because the available expression 

for the equilibrium concentration of Schottky defects in KCl, n = 6.7:10?:- 
‘exp [— 13900/7] cm-? (18), appears somewhat doubtful (!*). The correction 
for the long-range Coulomb interaction between the two vacancies (ref. (2°), 
pp. 300, 307-309) would reduce somewhat the ratio m/n both in NaCl and 
KCl, but not greatly owing to the small concentration of vacancies involved. 
The available experimental evidence on the process of matter transport 

in alkali halide crystals is compatible with the existence of an equilibrium 
concentration of vacancy pairs in the intrinsic region of the order of the con- 
centration that we compute for NaCl, but unfortunately the precision of the 
data is insufficient to allow firm conclusions to be reached (!°). Some evidence 
on the role of vacancy pairs is provided by the studies of MAPOTHER, CROOKS 


0.4 


02 


(15) The replacement of the Pauling polarizabilities by the TKS polarizabilities 
has practically no effect also on the ground-state binding energy of the impurity com- 
plexes, and the values computed by Bassant and Fumr (8) with the Pauling polariz- 
abilities remain valid. For Sr++ in NaCl and KCl, the new values are 0.43 eV and 
0.38 eV respectively, while the values reported by Bassani and Fumi are 0.45 eV 
and 0.39 eV. 

(1?) F. Bassani and F. G. Fumi: Suppl. Nuovo Cimento, 1, 114 (1955). 

(18) H. KELTING and H. Wirt: Zeits. f. Phys., 126, 697 (1949); F. SEITZ: Rev. 
Mod. Phys., 26, 7 (1954). 

(19) R. J. MAURER: private communication. 

(2°) A. B. Liprarp: Handb. d. Phys., 20, 246 (1957). 
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and MAURER (?1), and ASCHNER and MAURER (??), on fhe validity of the Ein- 
stein relation between the diffusion coefficient and the ionic conductivity in 
NaCl, NaBr and KCl crystals, and by the measurements of ScHAMP and KATZ (23) 
on the diffusion of Br- ions in NaBr. MAURER and his coworkers (21:22) report 
that the measured tracer diffusion coefficient for the cation in NaCl, NaBr 
and KCl in the intrinsic temperature region equals within the experimental 
accuracy the diffusion coefficient for the cation D, calculated from the mea- 
sured conductivity o by the Histein relation 


kTo 
(14) Di =a tes 
where t, is the transport number for the cation, and N is the number of molecules 
per em°. Actually, owing to correlation effects, the form of the Einstein rela- 
tion to be used to calculate the tracer diffusion coefficient for the cation D* 
is (25) 
kTo 


(15) D = 0.18 te 


Thus the measured tracer diffusion coefficient for the cation in NaCl, NaBr 
and KCl (2122) appears to be somewhat larger than it should be, suggesting 
that the vacancy pairs play some role. Similarly, if one keeps account of the 
correlation factor, the measurements of ScHAMP and KATZ (2°) show that 
the tracer diffusion coefficient of Br- ions in the intrinsic temperature region 
is of the order of one tenth of the diffusion coefficient deduced from the ionic 
conductivity, which refers essentially to the Na+ ions. This is quite compa- 
tible with the existence of an equilibrium concentration of vacancy pairs of 
the order of 10%, especially if one considers that the mobility of a vacaney 
pair in an alkali halide crystal may be smaller than the mobility of an isolated 
cation vacancy (2°22). The result of DIENES (1°) that the theoretical migration 
barrier for a vacancy pair in KCI is lower than the theoretical migration bar- 
rier for an isolated cation vacancy in NaCl, as calculated by Mort and LITTLE- 


(21) D. MapoTHER, H. N. Crooks and R. J. MAURER: Journ. Chem. Phys., 18, 
1231 (1950). 

(22) J. F. ASCHNER: Self-Diffusion in Sodium and Potassium Chloride, Ph. D. Thesis 
(University of Illinois, 1954). 

(23) H. W. ScHamp and E. Katz: Phys. Rev., 94, 828 (1954). 

(24) J. BARDEEN and C. HERRING: Imperfections in Nearly Perfect Crystals, edited 
by W. SHocKkLEYy (New York, 1952), p. 261; K. Compaan and Y. Haven: Trans. 
Faraday Soc., 52, 786 (1956); W. D. Compron and R. J. MAURER: Journ. Phys. Chem. 
Solids, 1, 191 (1956). 

(25) Y. HAVEN: Report of the Conference on Defects in Crystalline Solids held at 
Bristol in July 1954 (London, 1955), p. 261. 
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TON (4), cannot be considered as really probatory, also because the calcu- 
lation of migration barriers, especially for negative ions, is yet the less settled 
side of lattice calculations on point imperfections in the alkali halides (ref. (5), 
Sect. 3). 


It is a real pleasure to acknowledge the friendly discussions that the 
Authors have had on the subject-matter of this paper with Prof. SEITZ, and 
with Dr. HAVEX and Dr. LIDIARD, during the recent Varenna Summer School 
on Solid State Physics. The Authors are also indebted to Prof. MAURER for 
informative correspondence. 


RIASSUNTO 


Il modello di Born e Mayer dei solidi ionici è utilizzato per calcolare l’energia di 
legame di un posto vacante da ione positivo e di un posto vacante da ione negativo 
alla loro minima distanza reticolare in cristalli di NaCl e di KCl: i risultati sono 
rispettivamente 0.6 eV e 0.7 eV. Per il NaCl, questo porta ad una concentrazione di 
equilibrio di coppie di posti vacanti variabile fra il 10% ed il 20% della concentra- 
zione di equilibrio di difetti Schottky nell’intervallo fra 500 °C e 700 °C. I risultati 
sperimentali disponibili sono compatibili con l’esistenza negli alogenuri alcalini di una 
concentrazione di equilibrio di coppie di posti vacanti di questo ordine di grandezza 
nella regione intrinseca di temperature. 
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LETTERE 


ALLA REDAZIONE 


(La responsabilità scientifica degli scritti inseriti in questa rubrica è completamente lasciata 


dalla Direzione del periodico ai singoli autori) 


On the Motion of Charged Particles 
in the Complex-Symmetric Unified Field Theory. 


J. MOFFAT 


Department of Mathematics, Imperial College - London 


(ricevuto il 25 


In the complex-svmmetric unified 
theory (*?) the equations of motion of 
charged particles were derived from the 
field equations by an application of the 
EINSTEIN, INFELD and HOFFMANN ap- 
proximation method (*:4), and the tech- 
nique for deriving equations of motion 
introduced by INFELD (5). These equa- 
tions of motion contained the Lorentz 
force. 

However, at each stage of approxi- 
mation in the general theory, the equa- 
tions of motion proved to be incom- 
patible. We wish to pomt out that by 
taking into account the conditions of 
integrability of the field equations this 
incompatibility is removed. 

In the first treatment of this problem 
we used a source tensor expressed linearly 
in terms of the Dirac 6-function. The 


J. Morrar: Proc. Camb. Phil. Soc.,-53, 

473 (1957). 

(3) J. Morrar: Proc. Camb. Phil. Soc., 53, 
489 (1957). 

(3) A. EINSTEIN, L. INFELD and B. Horr- 
MANN: Ann. Math., Princeton, 39, 66 (1938). 

(*) A. EINSTEIN an? L. INFELD: Can. Journ. 
Math., 1, 209 (1949). 

(6) L. INFELD: Acia Phys. Pol., 18, 187 (1954). 
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Agosto 1957) 


conditions of integrability of these field 
equations can be taken into account—and 
the incompatibility removed—by apply- 
ing the «dipole procedure » introduced 
recently by INFELD and PLEBANSKI (°). 
The new treatment is based on the com- 
plex-symmetric field equations of empty 
space (*) 


(1) R,, =0. 
These field equations can be written 
(2) Diy + 2A, = 0. 


The ®,, are certain linear expressions 


in Y, Where y,, is defined by 
‘ TRE Ly n° 
(3) Vuv — Euv T 299) Exp . 
and 
f Mmn = —Omns Nom= 93 No 1; 
(4) | > 
Env — Guy a Muy : 


g,yis the complex-symmetric fundamental 


(6) L. INFELD and J. PLEBANSKI: Bull. Acad. 
Pole Sci. 14,768 (1956). 

(*) Greek indices run from 0 to 3 and Latin 
indices from 1 to 3. 
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tensor (*) (7) 


A 
(5) Iuv = "da + Iuv ti 
The 4, are defined by certain linear and 


uv 
non-linear expressions in y,,. 

Let us take a two-dimensional closed 
surface surrounding, say, the s-th sin- 
gularity. We take a surface integral over 


this surface and obtain 


s 


(6) Je + 2A,,)n,dS = 0, 


where n, are the components of the 
normal unit vector to the surface. Owing 
to the mathematical structure of @,, 
we know that [®,n,dS is identically 


MT: 
zero (4). Therefore, we are left with 


(7) fran as = (Ne 


These equations tell us that the surface 
integral is independent of the shape of 
the surface. 

Separating the field equations (2) 
according to the E.I.H. approximation 
method we obtain 

[ Poo + 2400 = 0, 


2l—2 21-2 


D ont 2A ma = 0 > 
(8) ep 


Pin ne vie =0. 
21 21 
The integrability conditions of these field 


equations for the known ,,s are (be- 
cause of the structure of ®,,): 


Ss 


‘A PS == (() 


| 

| 21—1 21—1 
Sa. 

| A * A : 

| f= | 2 Ayn tm AS = 0, 


(*) The * and /\ denote real and imaginary 
quantities respectively. 

(7) J. Morrat: Proc. Camb. Phil. Soc., 52, 
623 (1956). 


and 

s 

s 

ob = 3 , = 

mu 2A Gn Ny ds = 0 > 
21 


s 


s 

A A 

Y > / y 

CARRA nds =. Oe 
24 


As pointed out by EINSTEIN and 
INFELD (4), we can only insure the in- 
tegrability of the field equations (8) by 
using the « dipole procedure ». The first 
set of equations (8) is of the form 


21-2 1-2 
(11) A A 
Y00,rr 2 Ago - 
21-2 21-2 


To force the integrability conditions (9) 
and (11) we must add to *»,, gravitational 


21-2 


pole and dipole solutions and to D 00 


21-2 
electric pole and dipole solutions. The 
additional single pole and dipole expres- 
sions contained respectively in *»,, and 
SR 


s 


| Saratoga). 


| gad) ele 20-2 
(12) 4 Fg s i 
SET 
I \s=1 20-2 2-2 


After the surface integrals are cal- 
culated we find that the integrability 
conditions are identically satisfied if we 
choose 


| 4m = *C,, 

| 2-1 21-1 
(13) oe 

| DR 

21—1 20-1 

and 

| RD, = MDS 

| 21 21 
(14) 


Let us now remove all additional 
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gravitational dipoles by putting 


(15) Di Soe On 
lA 


21-2 


On differentiation with respect to the 
« auxiliary time » 7 we get 


(16) Yas, = > 2*0,, = 0, 
l 26 n 21 


which are the 39 equations of motion 
of the charged particles. In addition we 
have 


a 5s a $ 
ARSA 
1 s=1 2! 1 s=120 

and these equations determine the elec- 
tric dipole moment of the system. By 
calculating *C, and O, we find the m’s 
and e’s from (13). 

Thus in the case of the complex- 
symmetric field equations gravitational 
and electromagnetic dipoles are .added 
at each stage of approximation to insure 
integrability. After the field equations 
have been solved the gravitational di- 
poles are equated to zero and this re- 
stricts the motion of the charged par- 
ticles in the form of 39 equations of 
motion. The additional electromagnetic 
dipoles cannot be removed, because this 
would lead to inconsistent equations of 
motion, and in contrast to the gravi- 
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tational dipoles these solutions are phy- 
sically meaningful (8). At each stage of 
approximation the field equations deter- 
mine the motion of the charged particles 
(by exclusion of gravitational dipoles), 
and specify in addition the type of sin- 
gularity which is consistent with the 
motion; these equations of motion are 
compatible. 

If we neglect non-linear gravitational 
effects and cross products in mass and 
charge we obtain in the sixth order the 
equations of motion previously deduced 
(see reference (!) Eq. (12.5)) together 
with additional electromagnetic dipole 
forces, and these equations of motion 
contain the Lorentz force. 

Fuller details about this work and 
the calculation of the equations of mo- 
tion will be published shortly. 


I am grateful to Professor W. H. 
McCrea, Dr. W. B. Bonnor and Dr. D. 
J. CANDLIN for valuable criticisms and 
to my colleague Mr. R. KERR for many 
stimulating discussions and suggestions. 


(8) L. LANDAU and E. LIFSHITZ: The Clas- 
sical Theory of Fields (Cambridge, Mass., 1951). 
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Flux of Slow u-Mesons and Protons Near the Geomagnetic Equator. 


A. SUBRAMANIAN, S. NARANAN, P. V. RAMANAMURTHY, 
A. B. SAHIAR and SIDDHESHWAR LAL 


Tata Institute of Fundamental Research - Bombay 


(ricevuto il 29 Agosto 1957) 


In the course of an experiment on magnetic equator have been compared 
cosmic-ray particles stopping in a mul- with those obtained by others (1%) at 
tiplate cloud chamber, it was possible to higher latitudes near 50° N, to indicate 
evaluate the vertical flux of slow u-me- the latitude effects. 

The experimental arrangement as 

e. 000 shown in Fig. 1, consists of a circular 
cloud chamber, 45 em in diameter, with 

a 10cm illuminated depth, fitted with 

eleven quarter-inch thick brass plates, 

each equivalent to 4.2 g em? of air. The 

chamber was triggered by a coincidence 

in the Geiger-Müller counter telescope 

ABC above the chamber, accompanied 

by an anticoincidence from Geiger-Miiller 

trays D and E. Counters in tray E were 

= 8 000 used to reduce spurious triggerings due 


E : 
©) È 3 È to side-showers. The half-angles of the 
counter telescope were 11° in one plane 
(ue) 5 . 
€ = and 4° in the perpendicular plane. 
The particles stopping in the plates 
0 of the chamber were classified as electrons, 
pr u-mesons (*) or protons, from the visual 
<ocm estimates of ionizations of their tracks 
0 LR in the various compartments before stop 
Fig. 1. — Fxperimental arrangement. 


() B. Rossi: Rev. Mod. Phys., 20, 537 (1948). 

(2) M. CoNversi: Phys. Rev., 79, 749 (1950). 
sons and slow protons in air at Oota- (8) K, W. OGILVIE: Can. Journ. Phys., 33, 
camund (11.5° N, 76.7° E), 7500 ft above 746 (1955). 


i i RIC oo Na nati (*) The number of locally produced 7-mesons 
S id : Je o ae ; ; i E i 
sea level, an ADOU AN geomagnenle in our sample of stopping particles was estimated 


latitude. These flux values near the geo- to be negligible. 
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ping, and their range in brass. Protons 
were easily identified since they ionized 
>31;n in at least three compartments 
before stopping. For u-mesons, if they 
traversed the plate vertically and stopped 
at the bottom of it, the ionization in the 


TALI 


in brackets in Table I. We have not 
been able to separate from this sample 
of protons, the possible very small num- 
ber of K-mesons which on stopping did 
not give rise to any visible secondaries. 

Table I gives the details of the data 


TABLE I. 


x 


Flux values of slow u-mesons and slow protons near the geomagnetie equator (2° N). 


| Number | 
Particles of particles | 
| observed | 
En si - STO Ca 
u-mesons 272 + (49) 
| Protons 94 + (18) | 


and has been included in the range limits. 


compartment above this quarter-inch 
plate could be as low as 21,4,, while 
that for an electron would be minimum. 
Due to the lack of comparison tracks of 
known ionization in the neighbourhood 
of the track of a stopping particle, it was 
difficult to distinguish visually between 
ionizations of 2Zni, and Inin in all con- 
ditions of chamber photography and 
throughout the illuminated region of the 
chamber. Therefore the electrons were 
separated from the u-mesons by con- 
fining ourselves to particles stopping in 
the last four plates. The electrons could 
be identified by their scattering in the 
top seven plates, assuming that they did 
not loose energy by an undetected radia- 
tive process before stopping. For the 
determination of the proton flux, we have 
restricted ourselves again to the last four 
plates in order to separate from our 
sample, particles heavier than protons, 
such as deuterons, tritons etc. The num- 
ber of particles which could not be unam- 
biguously identified as protons or par- 
ticles heavier (*) than protons is given 


(*) Seven. definite cases of particles heavier 
than protons have been identified and so have 
not been included in the bracketted figure. 


Leni 
tesi 


Effective time 
of observation | 
| 


| 
Flux (*) 105 
(particles g-! air s-! sr!) 


(hours) | 
1072.2 | (3.7-=-4.4) + 0.2 
1072.2 | (0.9-=1.1) + 0.1 


(*) This flux refers to a range interval 47+ 64 g em? air equivalent. The contribution 
of the matter in the roof, counter walls, chamber walls etc., is 17.4 g em ? air equivalent 


obtained, and the resulting flux values. 

The uncertain cases of mesons bra- 
cketted in Table I, are due to the diffi- 
culty in identifying them as y-mesons 
or electrons in an unambiguous manner. 
In computing the area x solid angle factor 
for the acceptance of particles stopping in 
the bottom four plates, we have taken 
into consideration the loss of particles 
through scattering out of the illuminated 
region of the chamber, gain of particles 
through scattering into the illuminated 
region of the chamber, and the reduction 
in the solid angle of acceptance due to the 
dead space between adjacent counters in 
the trays A, B and C. We have also cor- 
rected for the loss of u-mesons due to 
some of the positive mesons stopping in 
the plates and decaying into positrons 
which discharge the anticoincidence tray 
below. The corrections due to scattering 
and the decay (for mesons) decrease the 
area x solid angle factor by about 35% 
for u-mesons, and about 9% for protons. 

Comparing the above limits of the 
slow u-meson flux with the value 
1.2-10-5 particles g7! airs sr! calcu- 
lated for 800 9 em? (+) residual pressure at 


(+) At Ootacamund 7500 ft a.s.l.) the re- 


sidual pressure is 800 g cia ?. 


ME? A. SUBRAMANIAN, S. 


50° N in the range interval 47 64 g em? 
air equivalent, we are led to a latitude 
effect ranging from 2.6 to 3.4. In the 
calculation of the above flux value 
at 800g em? at 50° N, we have used 
the differential range spectrum of slow 
u-mesons given by Rossi (1) and the 
data on the rate of variation of slow 
u-meson intensities with altitude given 
by Convers! (?). In Table II we have 
collected the latitude effects of slow 
u-mesons observed by different workers 


NARANAN, P. 


V. RAMANAMURTHY, A. B. SAHIAR ETC. 


effect suggested by OLBERT (7) and the 
geomagnetic latitude effect due to dif- 
ferent energy cut-off’s in the primary 
energy spectrum at the top of the atmo- 
sphere. The atmospheric latitude effect 
between 50° N and 0° as calculated by 
OLBERT (8) is of the order of only 10%. 
Hence the major portion of the latitude 
effects observed by various workers must 
be due to the geomagnetic latitude effect. 

We have given in Table III, the lati- 
tude effect, between 2° N and 50° N for 


TABLE II. - Latitude effects of slow u-mesons. 
| Residual | Residual 
| Experi- | Fane es | mos- : ’ 
| er ‘of mesons) - Latitudes Latitude 
Observers mental | b a pheric a fect 
technique °?S°FV®C | pressure SRO vari 
| (gem? | SH; 
2 | (gem?) 
| alr eq.) | 

| | | | 
| CONVERSI (? i[Delayed | | 
| ‘ C) |roincidence| 100117 306 | 59°N,  9°N | 1.87-40.12 
| DA pes eee | | 40 °N 29° N | 1.45-+0.1 
| DAVILA-APONTE : fGM | 3 i 

vILA-APONTE and |] of G 100 | 1030 |51°N, 29°N| ~1.35 

DeL Rosario (1) | [counters | 
| Kanexo et al. (5) . || 1545 | 1030 |50°N*, 24°N | 1.9:2.2 
| Tsar-Cuu and Max || Nuclear op ) 881 602 Net Ono NG 355.9 

Moranp (5) . ] emulsions | 
| » » il leto Sh GeO eNom aS 3.0--4.4 
| | Multiplate | | 
| Present authors . . | cloud | 47-64 | 800 |50°N*, 2°N | 2.6-34 
| | chamber | | 
| 
| (*) The flux values at 50° N latitude were calculated from the literature (1). 
| (4) L. DeL Rosario and J. DAVILA-APONTE: Phys. Rev., 88. 998 (1952). 

(5) S. KANEKO, T. KUBOZOE and M. TAKAHATA: Journ. Phys. Soc. Japan, 10, 915 (1955). 

(€) Tsar-Cau and M. MoRAND: Compt. Rend. Acad. Sci., 239, 415 (1954); 235, 1502 (1952). 


and those deduced by us for latitudes 
where only flux measurements have been 
reported. We have deduced the latitude 
effects by comparing these flux measure- 
ments with the flux values calculated at 
50° N for the appropriate atmospheric 
height and the relevant residual range 
interval. 

The latitude effect of slow u-mesons 
arises due to, the atmospheric latitude 


slow protons at 800 g cm-?, by comparing 
our flux value as given in Table I, with 
the value expected at 800gcm-? at 
latitude 50° N from the summary of the 
data (*) available on the slow proton 
intensity at sea level. We have used 


(7) S. OLBERT: Phys. Rev., 96, 1400 (1954). 
(8) S. OLBERT: M.I.T. Technical Report, 
61 (1954). 
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TaBLe III. — Latitude effects of slow protons. 

| 
Momentum | Residual | an | | 
Observers interval | pressure ANI a AUCUN 

(MeV/c) |(gem-?) compare effect 
iC re 7 | 
| = | | 
| CONVERSI (2) Py, en 984-1 043 306 59° N; 29° N=), 3.2--0.5 “4 
| Tsar-CHU and M. Moranp () | 724: 802 856 BO°N#*, 4°S-| (5-24 
| Present authors . 762 830| 800 | 50°N*, 2° N | 275.5 | 


(*) The flux values at 50° N were calculated from OGrLvir’s (**) data. | 


an absorption mean free path of (122 + 
+ 12)g cm? fer slow protons in air as 
measured by OcILvIE (*). For protons in 
the momentum interval 752 830 MeV/c 
(which corresponds to the residual range 
interval 47-64 g cm-? air equivalent), a 
vertical intensity of (0.7--1.0)- 107° cm? 
s1sr1(MeV/c)! is to be expected (39) 
at a latitude — 50° N, at 800gcm??. 
Taking our flux value per unit momentum 
interval as ((2.0--2.4) +0.2) -10-7em-? s+ 
sr! (MeV/c)-1 we are led to a latitude 
effect ranging from 2.7 to 5.5 for slow 
protons. We have given in Table III, 
the latitude effects of slow protons, ob- 
served by different workers, and those 
deduced by us for latitudes where only 
flux values have been reported. 

Our flux value of ((3.7 + 4.4) + 
+0.2)-10-6 particles g-! air s-! sr! for 
slow u-mesons in the range interval 
(47-64 g) em? air equivalent at 2° N 
obtained with a multiplate cloud cham- 
ber, compares well with the value (*) 
(3.3 + 0.6) - 10-9 g-! air s-! sr! at 4° 5, ob- 
tained in a lower range interval, viz. 
<5gcm- air equivalent by Tsar-CHU 
and M. Moranp (6) using nuclear emul- 
sions. However, for slow protons, in the 
momentum interval 752830 MeV/c and 
at the residual pressure of 800 g cm~?, 
we have a flux of ((2.0+2.4) + 0.2). 


I 


(8) Keo Wie 
555 (1955). 

(*) This is the extrapolated value for our 
pressure of 800 g em ?. 


OGILVIE: Can. Journ. Phys., 33, 


8 - Il Nuovo Cimento. 
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-1077 particles cm? st sr-1 (MeV/c)— at 
2° N, which is about half the value (+) 
(4.5 + 0.4)-10-7 particles cm? s- sr 
(MeV/c)! obtained by Tsar-CHU and 
M. Moranp (8) at 4° $, using nuclear 
emulsions. 


* OK OK 


We have great pleasure in thanking 
Professor H. J. BHABHA for his encoura- 
gement, and keen interest in the work. 
We would also like to thank Messrs. A. R. 
APTE and K. F. Dinsuaw for their 
assistance in running the chamber at 
Ootacamund. 


Note added in proof. 


BACCALIN, Bassi and MANDUCHI 
(Nuovo Cimento, 1, 657 (1955)) have ob- 
tained a vertical flux of (3.22 +0.22)-10-7 
particles cm? st srt (MeV/c)! for slow 
protons in the momentum interval 
(8901200) MeV/c, at a residual atmos- 
pheric pressure of 820 gcm-? and geo- 
magnetic latitude 47° N, using a Ceren- 
kov detector. Comparing our flux value 
of ((2.0--2.4) + 0.2)-10-? particles cm-2- 
-s1sr-1(MeV/c) with the value of 
BACCALIN et al. reduced to our altitude 
and momentum interval, we obtain a 
latitude effect of 3.0--4.9 between 47° N 
and 2° N. 


(*) This is. the value extrapolated for our 
pressure of 800g cm" and corrected for the 
difference in the momentum intervals. 
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Note on the Properties of the u-Meson. 


S. Hirokawa and H. Komori 


Institute for Theoretical Physics, Nagoya University - Nagoya, Japan 


(ricevuto il 29 Agosto 1957) 


In our previous work (*), we conclu- 
ded from the analysis of cosmic ray 
bursts that the anomalous magnetic 
moment (6) of the u-meson does not 
exceed eighty percent of the normal 
Dirac moment. In the study of bursts, 
the cascade theory of showers and the 
problem of fluctuations are important 
points which should be carefully exa- 
mined. As to the former, the Monte 
Carlo method (2) would be most reliable. 
In A, the ordinary cascade theory was 
corrected by the Monte Carlo method, 
though not so explicitly stated. Thus, 
at the footnote of page 741 of A, the 
factor of Christy and Kusaka 13.5 was 
reduced to 11.3. Hence, there would 
be no ambiguity in our cascade theory. 

The theory of fluctuations is discussed 
by Mirra (3). It is indicated that the 
Furry model by which we have dealt 
with fluctuations overestimates the fre- 
quency of bursts by a factor of about 1.5. 
Therefore, our conclusion that 6 is less 
than eighty percent of the normal Dirac 
moment needs not be modified. The 


(1) S. HIROKAWA, H. Komori and S. OGawa: 
Nuovo Cimento, 4, 736 (1956), cited hereafter 
as «A». 

(?) R. R. WILSON: Phys. Rev., 86, 261 (1952). 

(*) A. N. MITRA: Nucl. Phys., 3, 262 (1957). 


possibility of spin — 3 of the u-meson 
is also discussed by the same author (8). 
On the other hand, the bremsstrahlung 
cross section in the case of spin — 3 
obtained by MATTHEWS (4) exceeds the 
corresponding spin — + cross section 
by a factor at least 3 at the photon 
energy 101° eV. At the photon energy 
104 eV, the ratio is well over 10 
and still increases rapidly with increas- 
ing energy. Because the spin — + 
bremsstrahlung in this energy range 
accounts for nearly all of the bursts, it 
is expected that the curve of spin — 3 
bursts frequency in Fig. 2 of A would 
almost coincide with the curve of spin — + 
and ò=1. Since the ambiguity of our 
cascade function is not significant and 
since the uncertainty of the Furry mo- 
del is of the order of 50%, we con- 
clude that the possibility of spin — 3 
is excluded. 

Recently, the angular distribution of 
electrons produced by the decay of the 
polarized u-mesons has been observed (5) 
in order to prove the non-conservation 
of parity. Simultaneously, the variation 


(4 J. MATHEWS: Phys. Rev., 102, 270 (1956). 

(5) T. Corrin, R. L. GARWIN. L. M. LE- 
DERMAN, S. PENMAN and A. M. SACHS: Phys. 
Rev., 106, 1108 (1957), 
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of the counting rate of the decay electrons 
is also investigated by varying the mag- 
netic field. From these experiments, it 
is concluded that the gyromagnetie ratio 
of the u-meson is 2.00+0.01. This 
means (°°) that the spin of the u-meson 
is + and that 6 is very small. 


(5) The anomalous magnetic moment, 
— 0.4 <8<0.2 is obtained by Rawitscher 
making the assumption of spin —4 for the 


In conclusion, the assumption on the 
spin of the u-meson made in A is justified 
and no modification is necessary of our 
previous conclusion. 


u-meson, in which the cross section of the 
u-meson pair creation by photons is estimated, 
G. H. RAWITSCHER, preprint. We wish to thank 
Dr. G. H. RAWITSCHER for sending us the manu- 
script before publication (Phys. Rev., 107, 274 
(1957)). 


ERRATA - CORRIGE 


S. HIROKAWA, H. Komori and S. OGAWA — On the Anomalous Magnetic 
Moment of u-Mesons: Nuovo Cimento, 4, 736 (1956). 


Line 15, page 743 of A: «uncertainty of the cascade function » should read « uncer- 


tainty of the fluctuations ». 


Lines 9-14, page 745 of A: «But his reduction may not be...» should read as follows (1): 
«According to Carmichael, the number of ion pairs in the experiment of Driggers 


is small; 3.0-105—1.2-106 ion pairs. 


So, the energy range covered by Driggers’ 


experiment is different from that of Schein and Gill». 


(1) W. D. WALKER: private communication. We are grateful to Professor W. 


D. WALKER for 


communicating us his valuable opinion on this point. 


8* - Il Nuovo Cimento. 
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A Remark on a Paper by Barasenkov. 


C. HAYASHI 


Department of Nuclear Science, Kyoto University - Kyoto 


(ricevuto il 24 Settembre 1957) 


In a recent paper BARASENKOV (') has made a remark that the integrability 
condition is not satisfied in the author’s (2) Hamiltonian formulation of quantiz- 
ed fields with a non-local interaction. Such a conclusion cannot be true since 
the Hamiltonian in this formulation has originally been derived in accordance with 
a general scheme which contains a prescription completely equivalent to theinte- 
grability condition. It will be shown in the following that the calculations made 
by BARASENKOV are in fact misleading. 

For the sake of comparison, the notation of BARASENKOV will be used. The 
canonical variables A,[x; 0] and wy,[x; 0] obey the equation of motion 


SA gr; 0] _ i[H(2'/0); Al@;0) ]. 


(1) do(x') 


The Hamiltonian has been expressed as space-time integrals of products of the 
operators Ay, y, and factors which depend on o explicitly (2). Then, we have for 
the Hamiltonian 


(2) OS ee Al 
E : do(x') Co(x') 


where 5/30 means a total derivative, and 0/00 is a derivative with respect to o 
which is included explicitly. The last term of (2) has been overlooked by Bara- 
SENKOV. The contribution of this term to [à/50(x"); 3/do0(x')]Ap[x; 0] just cancels 
the other contributions as calculated by BARASENKOV (1) (*). Using eq. (1) and (2) 


(1) V. S. BARASENKOV: Nuovo Cimento, 5, 1469 (1957). 

(2) C. HAYASHI: Prog. Theor. Phys., 10, 533 (1953); 11, 226 (1954). 

(*) His expression (29) itself should be corrected as follows. The last two terms on the right 
side of [07 ,(x'/0)/00(x"); Ale; 61], that is, Am (x, — 21) Vols; O1 Volts; 01+ Ado — 22) vd [X13 ol ve (21; 0] 
should be replaced by A4m(w, — x3)wd Les olpolts: 014 Am(s — x1) vi Lt: olwlts; ol. The expression 
for [S/So(x"): S/So(x')] Alt ; c] should also be corrected correspondingly. 
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we have 


So(a")’ So(x') 


3 8 
(3) | = facts: al 


Sa de TR ei 


Li È 0o(x") 0o(x') 


and an expression of the same form for pol; 0]. Each commutator on the right side 
of (3) has been calculated in the g?-approximation (1). It is easily found that (3) 
really vanishes in this approximation. More generally, à direct calculation gives 


Sc nio 


(4) 5 , È " Wa 
o(x') Oo(x”) 2 


which confirms that the equivalent integrability condition, 


(5) ae a + [A(x'/0); H(x'/o)] = c-number , 


is satisfied. 


In conclusion, the author thanks Prof. S. Hrroisur for valuable discussions. 
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Theory of the Dirac Field as a Continuous Assembly 
of Small Rotating Bodies. 


T. TAKABAYASI (*) 


Institut Henri Poincaré - Paris 


(ricevuto il 24 Settembre 1957) 


The Dirac spinor y (in c-number theory) defines a Lorentz frame at each point 
of space with the set of four 4-vectors all real and mutually orthogonal. They are 
quantities related to y in bilinear manners, i.e. (*) 


A 


(1) Sy ipyly, S,= ipy°yl y, 


and the real and imaginary parts af and «? of the complex vector 


(2) E,= py'y = ay + ie, 
where y° denotes the charge conjugate spinor defined by 


(3) y= 097, 


(4) C't0=1, oT=— OC, (ye)? = — OHO. 


Accordingly £, may also be written as 


LH 


= 
de 


(5) E, = — yt yety . 


If we introduce the complex vector 27 which is imaginary conjugate (-) to en 


(*) On leave of absence from Nagoya University (Japan). 

(+) p=ytyt. In this paper + signifies adjoint, * complex conjugate, and 7’ transpose. A 4-vector 
is real when its space components are real while its time component is pure imaginary. Greek 
suffixes run from 1 to 4. Precisely, Sy is a pseudovector. 


(7) Here «imaginary conjugate » means that Si Eh= for u=123, = — Zi for u=4. 
Such £} can be written &X = ptyty4ytay* = — yyy’, and is shown to be a vector. 
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oO) d x? : : 
x, and a, are written in the form: 
1 1 
4 Qa) = PX (12) DIE = TX 
6) cy = (Ent Zu): ca pe 
Actually, we can verify the identical relations 
@)\2 (2)\2 G2 2 2 2 
(au) = (x, ) = S, = 2 ’ Sto ? 


(7) 


Mib ONE OE, _ Log ANT è 7 
| cu Cy = Sy a, S, = a, Ju = ch, S,—=S,8,=0, 


where o is the scalar quantity defined by 


(8) 0 = (Q?+ 22), 
with 
(9) Q= py, 2 = ipyy 


Incidentally we note that {9, 0} defines the other variable 0 (pseudoscalar) by 


(10) te 0 — 0/2. 


A 


é : 3 2 ; : : 2 
If we introduce unitary vectors af, (£—1, 2, 3, 4) in the directions of x, x?, S, 


and S, by 


È) È) - 
( op = oa, (E = 1, 2) 
11 
ae è (3) g (a) 
| a > 0a, Ar Loa, ’ 
then (al?) satisfies 
€ (ES) N (È) (9) 
(12) ay a, > Onp > an dy = den ? 


and exactly represents a Lorentz frame. 

Furthemore we can show that the Dirac field y is represented equivalently by 
this Lorentz frame assigned to each point of space, together with two more va- 
riables o and 0. 

Also we take into account the fact that the existence of a Lorentz frame (al?) 
is equivalent to the co-existence of a translation with a velocity v, (*) and a 
rotational orientation of spatial axes represented by a rotation matrix (C;,), which 
naturally obeys 


(13) (OPA OP = Ò je . 
(*) vy is a unitary 4-velocity satisfying vi =— 1, and is related to the 3-dimensional velocity 
V through 
v;=vV;le , %=i% , (Fe): 


Latia suffix runs from 1 to 3. 
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os E . 
Explicitly ae are expressed in terms of v, and C,, as follows: 


| Oy = —Wy Gia 100 
(14) Cine, 
ik nee 
1 — iv, 


We have thus obtained the representation of the Dirac field in terms of the 
set of variables 


(15) to (Ca); 0, 0}. 

The above result shows that we can regard the Dirac field as a continuous 
assembly of very small rotating bodies (*) whose rotational orientation and trans- 
lational velocity are specified by (C,,) and v,, respectively, @ means simply the 
distribution of the number density of such bodies in the rest reference frame and 0 
is an additional internal variable possessed by the body, other than the rotational 
degrees of freedom. 


At this point we remark the following. Under a gauge-transformation (+ 
le] t=) =) 


li RE 


ie 
(16) ARTEATRO AE y —> exp a 
ie 


o, 0, a? and a‘? are invariant, while a‘? and af are not so. It is therefore profitable 


to revise the definition of £, slightly so that it becomes gauge-invariant. Namely 
we separate a phase factor such that 


(17) E,= exp (24/1): E, 


where À is a redundant scalar function with the property that it transforms, under 
(16). like 

e DI 
(18) A->A+-4. 


lo 


Henceforth we adopt ni and corresponding do and a (€=1, 2), and write them 
simply as &,,, ae They are all gauge-independent quantities and retain 
the properties in the original definitions. The Dirac field is then represented by 
the set of variables, 

(19) : (a 0 0, Ae. 


HZ 


Finally the time evolution of the Dirac field is expressed as the equations of 


(*) By « rotating body » we only mean a certain body whose rotational orientation of internal 
rotation has a definite meaning, but do not assume further properties concerning the usual concept 
of rigidity. - 

(*) Ay is the poteutial of the electromagnetic field acting on the Dirac particle; 


lo) 


Ory 


en 


ZIE 
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motion for our variables (19). They are derived from the following Lagrangian den- 
sity (*), 


(4) (3) 5 


Lt a, 0a) sr 


2 he 
(20) P= — 9 {metal k, + cos 0) + 5 ( 


(3) : 
ay 0,9 _ VE yyy À 
(9) (6) 
+ PACA ay dra Ô y) 5) 
where the derivations are to be taken with respect to the arguments 


È) 
(21) {ay » ©» 9, AR hg == vus ? 


while k,, in (20) are dependent variables defined by 


Ds (2) 
ay 0,0 


e 
(22) mck,, = Au end 


ml & 


This is also gauge-independent and has the meaning of a particle momentum- 
energy 4-vector. The equations of motion following immediately from (20) contain 
the redundant variables 4,,, other than (aC: 0, 0,k,}, but 2,, can be easily elimi- 
nated therefrom. (+) Anyway we can actually verify that the set of our equations 
of motion resulting from (20) is equivalent to the original Dirac equation. 

The present theory has close relationship with the hydrodynamical represen- 
tation of the Dirac field previously constructed (1). The variables o, 0, v, and k, 
are exactly the same things as employed in the hydrodynamical scheme. Also a 
is nothing but the spin pseudovector w, used there. 


(*) Eyv,a is the signature of the permutation (urx2). 

(+) The equations of motion are considered in a separate note. 

(1) T. TAKABAYASI: Nuovo Cimento, 3, 233 (1956); Phys. Rev., 102, 297 (1956); Suppl. Prog. 
Theor. Phys. (Jepan,) No. 4 (1957). 
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Variation of Electron Mass with Velocity. 


F. C. CHAMPION 


King’s College, University of London 


(ricevuto il 4 Ottobre 1957) 


—FaRAGO and JANOSSY (!) have re- 
cently concluded that the experiments 
of CHAMPION (2) failed to distinguish 
between the formula of Abraham and 
that of Lorentz, on the grounds that the 
experimental error was too large. Tak- 
ing the most unfavorable case selected 
by them, where they represent the error 
as three times the standard error, amount- 
ing to 12% of the rest energy, I would 
point out that Fig. 1 of FARAG6 and 
JANOSSY’s paper indicates a ratio of 
300% in the values predicted by the 
two theories for (m,/m)? at the energy 
considered by CHAMPION. This was ex- 
pressed by him in the form of expected 
collision angles of 58° and 75.2° accord- 
ing to the two theories respectively, the 
experimental angle being found to be 
(75.2 + 1.5)°, with a difference therefore 
of over 25%. Details of thirteen other 
independent and technically favorable 


(1) P. S. FARAGO and L. JANOSSY: Nuovo 
Cimento, 5, 1411 (1957). 

(3) F. C. CHAMPION: Proc. Roy. Soc., A 136, 
630 (1932). 


instances of electron-electron collisions 
were given, all of which supported the 
relativity formula with the same order 
of experimental accuracy. 
While I am entirely in agreement with 
FARAGO and J'ÂNossy that the accuracy 
of collision experiments cannot compare 
with that attainable by spectroscopic 
techniques such as the fine structure 
doublet separation of hydrogen-like spec- 
tra, I must insist that with the most 
unfavorable estimation of the errors in- 
volved in the collision experiment, there 
is no doubt whatever that they show 
that if one wishes to predict experi- 
mentally how momentum and energy 
will be shared during the elastic col- 
lision of two electrons in the energy 
range (0.5—1) MeV, a much more ac- 
curate forecast will result if the relati- 
vistic formulae are used in preference 
to any other existing relations. Whether 
this conclusion « proves » the validity of 
relativity theory, or to what degree of 
accuracy any experiment has to be 
carried to «prove» any theory, is a 
matter of definition and opinion. 
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On the X-Hyperon-Nucleon Interaction Cross-Section - II. 


F. FERRARI and M. TONIN 


Istituto di Fisica dell’ Università - Padova 
Istituto Nazionale di Fisica Nucleare - Sezione di Padova 


(ricevuto il 22 Ottobre 1957) 


As is well known, an apparent discrepancy exists between the data on K_ 
capture in hydrogen (1) and in emulsion nuclei (2) in so far as the interaction rates 
leading to different hyperons are concerned. In a previous paper (3) it was 
supposed that the figures of — 0.4 for the ratio (A,_+%%)/(X£7 +2;) as given 
by ALVAREZ et al. (1) and < 3 for the ratio (Ag+2%)/(Zg+Yà) obtained from 
photographic emulsion experiments, might be due to a large cross-section for the 
process 


(1) ini ep NSA Ne, 


transforming, inside nuclear matter, a large number of Y's into A’s. The fitting of 
the data with such an hypothesis leads to a reaction cross-section oy, of about 
50 mb, while the cross-section osy for the process 


(2) Zz4N> MAN i=+,—) 


results much smaller. 

This value however may be greatly modified if one takes into account that 
the hyperons are produced in nuclear matter and consequently, as suggested by 
various authors (4), may be subjected to a mean nuclear potential: the hyperons 
in the nucleus will be considered as being inside a nuclear potential Vy to which 
the Coulomb potential V, has to be added for the charged hyperons (5). In what 
follows we will consider a nucleus described by an ideal Fermi gas defined by 
R=r,At, where r,—1.4:10-1 em, and having the binding energy of the last nucleon 
equal to 8 MeV. Furthermore we will assume that the K particles are absorbed 


(3) L. W. ALVAREZ, H. BRADNER, P. FALK, J. D. Gow A. H. RosENFELD and R. D. TRIPP, 
F. T. SOLMITZ: Nuovo Cimento, 5, 1026 (1957). 

@) W. F. Fry: Phys. Rev., 100, 950 (1955); D. W. FALLA, M. W. FRIEDLANDER, F. ANDERSON, 
W. D. B. GREENING, S. LIMENTANI, B. SECHI-ZORN, C. CERNIGOI, G. IERNETTI and G. POIANI: 
Nuovo Cimento, 5, 1203 (1957). 

(3) N. DALLAPORTA and F. FERRARI: Nuovo Cimento, 3, 743 (1957). 

(4) Padua-Venice Conference on «Mesons and recently discovered particles» 1957. 

(*) H. Capes: Phys. Rev., 107, 239 (1957). 
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by a single nucleon in the nucleus exclusively by the reaction 
(3) Kees NESY ar 


In fact, as is experimentally known, this reaction is the predominant one; the 
probability that the absorption involves two nucieons, 


4) Erbe 


being smaller than 10%. 

Taking into account the relative isotopic spin wave functions of the initial and 
final states of reactions (3), and the ratios observed in hydrogen by ALVAREZ 
ét alo (2), Le. 

WO sis Spar a A-9-9- 
DB Bec epee Ny RASO 20 


Sp Ap: dp + À 


we have that the ratios in nuclei are: 


(5) REA REL 
or 

SOLA 
(6) EVI 

=D ato 2, 


Let us now consider a hyperon produced at a point P distant @ from the center 
of the nucleus. As in reference (*) the probability that a charged 2° may be con- 
verted to 2° or A is taken as being related to the-strength of the reactions (1) 
and (2). However we have here the additional fact that the hyperon may be re- 
flected from the nuclear surface: the reflection probability (or the transmission 
coefficient ti(E)) is a function of the X-nucleus potential Vs, of the electrostatic 
potential V, and of the kinetic energy of the hyperon itself. 

The probability that a hyperon generated at a point P by a capture event, (*) 
will escape from the nucleus as X° or Y°, will be given by the following relations 
(see also formulae (4) and (5) of reference (°)): 


R+o foo) 
È exp [— t/Apx] 
7 ITs: = [dt W(R, 0,1) = 
(7) 3 = far Poe iui yy + ofl — P(R,R, La — T/T 
R-o ms 
ea i exp [— t/Zrx] | 
= po | LA 
re ass 5 crei 
TV 
and nd * 
R+o EE 
(8) IT de W(R, 0, 1) Is Sp + lx Ha n° 
R-o ; ù 


sin far PI (1 exp [—t/Aqy] | 
1+ 2f1—P(R, R, A,x)]0 — T/T) |f 
1 sio 


al 


(*) In the derivation of expressions (7) and (8) we have supposed that the number of 2° which 
are converted into X° is approximately equal to the number of £° which. are converted into ’. 
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In the above, P'(£) represents the energetic distribution function of a hyperon 
with charge à inside the nucleus; 1 — 7° gives the probability that a hyperon >‘, 
after a collision, does not transform into a A and may have sufficient energy 
to escape from the nucleus. The function P(R, 0, Z,y) which gives the escape pro- 
bability, has been defined by CLEMENTEL and Puppr (*) and 7” is given by: 


fae Pi(E)ti(E) 
Vy-V, 


ME 
fab Pi(E) 


LEA 


The introduction of the potential Vy allows one to fit the experimental data as- 
suming different values for the cross-sections os, oxy and for the hyperon-nucleus 
potential Vy. In order to discriminate between the various possibilities which 
are thus presented, it is convenient to consider other experimental data. For this 
purpose we have also considered the ratio 2 /X* which in hydrogen turns out to 
be X;/2; ~ 2, while from emulsion experiments it turns out to be Xg7/Z2d7 <2. 

The ratios (=87+ Agz)/(Log7-+ Ngz) and Loz/Zdz can easily be calculated from for- 
mulae (7) and (8). 

The results are plotted in Fig. 1 and in Fig. 2, as functions of the total cross- 
section o,» and for several values of the Y-nucleus potential: 


1) The full line has been obtained by assuming that the cross-section of re- 
action (1) is greater than the cross-sections of reactions (2), namely cyy=0.2 67». 


2) The dotted line is obtained if one makes the inverse hypothesis. oy, = 
— 02 CAM 


(9) E. CLEMENTEL and G. Puppr: Nuovo Cimento, 10, 197 (1953). 
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From an inspection of the figures and taking into account that the observed 
ratios are (Za + Agz)/(Zgz+Ugz) S 3 and Z77/Za<2 (these values have been cor- 
rected for the loss of =~), it can be seen that the X-nucleus potential must lie between 
~ 10 MeV and ~ 20 MeV, unless one assumes very small oy, or 0yy cross-sections, 
say smaller than 10mb. These limits may not be considerably modified if one 
considers a more refined nuclear model (see H. Capps, reference (°)). 

It is to be remarked that the hypothesis 1) leads to total cross-sections of the 
order of magnitude ~ 10 mb, while the hypothesis 2) can also lead (see Fig. 2) 
to greater total cross-sections. 

As a consequence, in this last case one can obtain agreement with the experi- 
mental data with a quite small cross-section for the processes (1), say oy, ~ 10 mb, 
while the cross-sections for the processes (2) can be of the same order of magni- 
tude as the geometric nucleon-nucleon cross-section, say cry ~ 40 mb. (We have 
assumed Vy — 10 MeV). It is interesting to compare these values with the cal- 
culated total S-wave hyperon-nucleon scattering cross-sections given by LICHTEN- 
BERG and Ross (’) (see Table I). 


Cross-sections in mb 


in MeV | Ap— Ap X+*p->X*p 
An-An |X n> n 


a È Sarzano 
x p>An|Z°p>Ap|X p>X p|2°p> Lp] x p— 2°n 
Ztn > Ap|Zn+An]|XZ*n> Stn} Xn > X°n|X*n+ 2p 


Tu i 15 
CNE 


| | 


== a 
| | 
| | 
LANDE 
Although a detailed analysis cannot be made on the basis of our simple optical 
model, it is seen that at last qualitatively these results are not in disagreement 
with those derived by LICHTENBERG and Ross from scattering data. 
* * * 


We are indebted to Prof. N. DALLAPORTA for many discussions and advices. 


Note added in proof: 


The most recent experimental values of the ratio Xx/Xy are considerably smaller 
than 2 (8). Such a value can be explained by considering also the processes 
Zn > Up, 


Zip > En, 


(7) D. LICATENBERG and M. Ross: Pion Contribution to Hyperon-Nucleon Forces (in press). 
6) F. C. GILBERT, C. E. VIOLET and R. S. WHITE: Phys. Rev. 107, 228 1957, K= Stack col- 
laboration; Padua-Venice Conference 1957. 
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which in our calculations have not been considered (see note at foot of page 2), 
If we include these processes and assume a cross-section ogey- = oyoy+ —15 mb, 
and take into account the ratios given in formula (5) and the energetic distribution 
function of the hyperon after the collision (the hyperons are mainly of low energy 
owing to the Pauli principle), we have that, with Vs—10 MeV, 80% of the X° 
produced by &° are captured while 90% of the X* can escape from the nucleus. 
With these hypotheses we have 


EX + Ax 


SFILA 


The figure of 1.2 is in good agreement. with the latest experimental data quoted 
above. A detailed discussion on these topics will be published. 


We wish to thank Dr. D. Evans and Dr. D. PROWSE for discussions. 
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A Search for the Existence of Asymmetries in Positive K-Meson Decay. 


A. WATAGHIN (*) 


H. H. Wills Physical Laboratory - Bristol 


G. ALEXANDER and R. H. W. JOHNSTON 


Dublin Institute for Advanced Studies 


(ricevuto il 23 Ottobre 1957) 


1. — Introduction. 


It is a consequence of the two com- 
ponent theory of the neutrino, that an 
asymmetry should exist in the angular 
distribution of the meson secondary to 
the K,, decay mode, if: a) the spin of 
the K-meson is different from zero, and 
b) the beam is polarized, and c) serious 
depolarization does not take place in the 
focussing and analysing magnetic fields. 

In the course of an experiment (1), 
the purpose of which was the determi- 
nation of the relative abundances of the 
K-particle decay modes and the study 
of the 3 body decay-spectra, records 
were kept of angular data relating to 
2697 K, decays. The angular distri- 
bution of the secondaries in this sample 
of 2697, which may be regarded as being 
random, has been analysed in this work. 

For reasons explained at length in (1), 
it is considered that the scanning bias 
is negligible. The number of K-particles 


(*) On leave of absence from Faculdade de 
Filosofia, S. Paulo University, with a fellowship 
from Conselho de Pesquisas of Brazil. 

(1) G. ALEXANDER, R. H. W. JOHNSTON and 
C. O’CEALLAIGH: Nuovo Cimento, 6, 478 (1957). 


coming to rest for which no secondary 
was visible is less than 2%. 
2. — Conditions of exposure. 

The experimental arrangements for 
the selection of the primary K-particles 


were those shown in Fig. 1 of BIRGE 


200 


100! 


a 
0 90 180 270 cru 


Fig. 1. — Angle « between initial secondary di- 

rection at edge of stack, projected on the 

emulsion plane. Distribution of values of « 
from 2697 Kt decays at rest. 


et al. (2). A beam of secondary particles 
from a copper target was selected at mo- 


() R. W. BIRGE, D. H. PERKINS, J. R. PE- 
TERSON, D. H. SToRK and M. N. WHITEHEAD: 
Nuovo Cimento, 4, 834 (1956). 
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mentum 409 MeV/c at 90° to the cir- 
culating proton beam of the Berkeley 
Bevatron by a system of quadrupole 
focussing magnets followed by an ana- 
lysing magnet. The field of the latter 
produced a deflection of the secondary 
beam through an angle of 45°. 

The emulsion block was exposed with 
the pellicles perpendicular to the plane 
containing the 409 MeV/c focussed beam 
and the initial 6.2 GeV proton beam, 
so that the momentum, as defined in 
the analysing magnet, was constant in 
each plate but varied slowly across the 
stack. 


3. — Depolarization of the beam. 


In the event of the K-meson beam 
being polarized in the production process, 
one might expect a certain degree of de- 
polarization to take place in the quadru- 
pole focussing magnets; particles tra- 
velling along different rays within the 
cone of acceptance will encounter mag- 
netic fields of widely differing strengths 
and orientations. However, over a wide 
central area of the magnetic lens-system 
the field is relatively small, so that to 
assume as an order of magnitude 104 G 
for the depolarizing magnetic field would 
appear to be reasonably conservative. 
The time of flight in the quadrupole 
focussing system is about 2-10-° s. Under 
these conditions one would expect pre- 
cessions of the order of 10° to take place, 
so that the beam would not undergo 
serious depolarization. 

The time of flight in the analysing 
magnet is of the same order of magnitude. 
The field in this case however being uni- 
form, one expects precession of all the 
particles to take place by the same 
amount, so that the percentage po- 
larization of the beam should not de- 
crease. 

The stray magnetic field in the region 
of the plate was found be ~ 50 G. After 
the K-meson has come to rest, the pre- 
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cession that one can expect due to this 
field in 10-§s (mean lifetime of the 
K-meson) is of the order of 0.5°. Thus, 
the depolarization effect arising from 
this cause may be neglected. 


4. — Possible spin effects. 


If the K-mesons have spin different 
from zero and the beam were polarized, 
some anisotropy in the angular distrib- 
ution of the particles emitted should 
exist unless the secondary particles from 
the K_, mode (abundance 23%) have 
possessed such asymmetry as to com- 
pensate the supposed asymmetry in the 
K,,, (61% abundance). Our results, given 
in the next section, on the identified K,., 
K,, and the total sample, do not lend 
support to such an assumption. 

It should also be remarked that, if 
the angular distribution of the x in the 
K_, mode is supposed isotropic (there is 
no reason at the moment to suppose 
that it is not), an asymmetry of ~ 15% 
in the K,,, part of the total sample would 
be clearly seen as a ~ 10% asymmetry 
in our total sample. 

It is of interest to note also that the 
observation of an asymmetry in the 
decay of u-mesons in emulsion (*), in 
spite of possible triplet muonium for- 
mation, suggests that K-meson effects, 
if they exist, ought not to be obscured 
by the analogous process of K-mesonium 
formation. 


5. — Experimental results. 


Our results, given in Tables I, II, III 
and in Fig. 1, are consistent with iso- 
tropy in space of the direction of emis- 
sion of the charged emitted particle. 
For the K, sample, a peak to valley 


(ERI GARWIN, L. M. LEDERMAN and 
M. WEINRICH: Phys. Rev., 105, 1415 (1957). 
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excess ratio > 10% in any direction in 
space should be seen either in the pro- 
jected « distribution (Table I, Fig. 1) or 
in the data relating to dip given in 
Table III. The angle « is the projection 
on the emulsion plane of the angle bet- 


A.. WATAGHIN, G. 


ALEXANDER and R. H. W. JOHNSTON 


ween the direction of the K-beam when 
entering the stack and the direction of 
the charged particle emitted in the decay. 
The « distributions for identified cases 
of K,, and K,, are given in Table II. 
The distributions of these samples, though 
smaller, are also consistent with iso- 


TABLE I. — Number of K,. A 
0° a 45° ‘ 348 TABLE II. 
45°— 90° 318 0 = Bp, Ka, 
90° — 135° 348 0°— 60° | 33 11 
135° —180° 316 | 60° — 120° | 28 16 
180° —225° 336 | 120° —180° 47 | 14 
225° —270° | 334 180° — 240° | 23 12 
270° —315° 359 240° —300° 28 12 
315° —360° 338 300° —360° 36. | 12 
Total 2 697 Total 195 77 
TABLE III. 
ID YA | | N umber Percentage | Dm ran 
(Processed emulsion).| (Experimental) (Experimental) | of isotropy 
0<d< 6 408 13.4 + 0.5 12.9 
6< 0 < 20 699 23.0 + 0.9 25.0 
20 < 0 | 1936 63.6 + 0.9 62.8 
Total 3043 100.0 100.7 
Toward glass 1169 
Toward surface 1153 
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6. — Discussion and conclusion. 


An investigation of the angular dis- 
tribution in space the secondary of 
436 cases of K,-decay has been carried 
out by M.I.T.-Johns Hopkins Group, (4). 
Further, F. and M. BRuUIN (5) have ana- 
lysed 500 K,, (the projected angle dis- 
tribution only). No evidence for ani- 
sotropy was found in either case. 

The above analysis based on 2697 


(4) Reported on behalf of M.I.T. - Johns 
Hopkins Group C. O’CEALLAIGH: 7th Rochester 
Conference Report. 

6) F. and M. BRUIN: Physica, 28, 553 (1975). 
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events enables an upper limit of about 
10% to be placed on the anisotropy with 
1% confidence on unresolved modes, or 
a — 15% limit on the K,, part of the 
total sample, if the distribution of the 
other modes is supposed isotropic. There 
is no indication of any asymmetry 
in the identified K,, and K_, beyond 
that expected through statistical fluc- 
tuations. 
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